General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




USE OF MAGSAT ANOMALY DATA FOR CRUSTAL STRUCTURE AND 
MINERAL RESOURCES IN THE U.S. MIDCONTINENT 


(E84- 10112) USE CF MAGSAT fi KC K ALJf DATA FOB N8^-2 

CEUSTAL 31BUCTUBE AND MIKE5AL EESCUBCES IN 
THE US MIDCONTINENT Final Frcject Beport, 

Dec. ^80 - Sep. 1983 (Iowa Uciv. ) 177 p Oncla 

HC A 39/M? A3 1 CS CL 38G G3/U3 33112 


Robert S. Carmichael 
Department of Geology 
The University of Iowa 
Iowa City, Iowa 52242 


October 1983 

Final Project Report for period Dec. 1980-Sept. 1983 


Prepared for 

GODDARD SPACE FLIGHT CENTER 
- Greenbelt, Maryland 20771 




USE OF MAGSAT ANOMALY DATA FOR CRUSTAL STRUCTURE AND MINERAL RESOURCES 
IN THE U.S. MIDCONTINENT 


Robert S . Carmichael 
Department of Geology 
University of Iowa 
Iowa City, Iowa 52242 


October 1983 

Final Project Report for period Dec. 1980 - Sept. 1983 


Prepared for 

GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 20771 










NASA 




]. Report N*. 


TECHNICAL REPORT STANDARD TITLE PACE 


12* Gov«rnrn*nt Accriiion Na. 


3* Recipient*! Catalog No. 


XTIirr3nT3uEti33“ 

USE OF MAGSAT ANOMALY DATA FOR CRUSTAL STRUC- 
TURE AND MINERAL RESOURCES IN THE U.S. MIDCON- 
TINENT 


5- ft « part Dot* 

October 1983 


6* Performing Organisation Code 


7, Aulhor(s) 


Robert S* Carmichael 


8- Performing Organisation Report No. 


9. Performing Organ! ration Name and Addrctt 

Department of Geology 
University of Iowa 
Iowa City, Iowa 52242 


70. Work Unit No. 


IT* Contract or Grant No. 

NAS 5-26425 


! 


12. Sponsoring Agency Nome and Address 

Goddard Space Flight Center 

National Aeronautics and Space Administration 
Washington s D.C* 20546 

Techa! cal, Monitor. ;.^Har,o.l.d-Qa^ r . gf 


13. Type of Repart and Period Covered 

Final Project Report j 
Dec* 1980-Sept. 1983 j 


Id. Sponsoring Agency Code 


75. Supplementary Hotel 


16. Abstract 

Magnetic field data acquired by NASA’s Magsat satellite is used to 
construct a long-wavelength magnetic anomaly map for the U*S. mid- 
continent* This aids in interpretation of gross crustal geology 
(structure, lithologic composition, resource potential) of the re- 
gion. Magnetic properties of minerals and rocks are investigated 
and assessed, to help in evaluation and modelling of crustal magne- 
tization sources and depth to the Curie— temperature isotherm. 


17. Kay Words {S^Ktd by AutW(*» 

U. S* Midcontinent 
satellite magnetic anomaly 
magnetite 

crustal magnetization 


18. Distribution Statement 


19. Security Closstf. (of this report) 

20<* Security Classif. (aftfiti page) 

21. No. af Pag«* 


Unclassified 

Unclassified 

176 

1 22. Pric* 


Far sale by die Clearinghouse for Federal Scientific and Technical Irlorasauan, Springfield, Virgin** 22151- 



PREFACE 


The objective of the work was to investigate the use and applicabi- 
lity of satellite magnetic anomaly data, acquired by NASA's Magsat mission 
of 1979-80, in reconnaissance interpretation of crustal geologic character. 
This would include major structure, tectonic development and evolution, pet- 
rologic (rock) composition, and thermal conditions in the lithosphere. The 
Magsat data provides long-wavelength magnetic information with a coverage, 
accuracy, and resolution not previously available. As such, the satellite 
data and maps constitute an important new aid for interpreting the crust 
of the Earth. 

The scope of the project has been to 

i) process the satellite magnetic data to yield an anomaly map of the 
central U.S. This has involved data selection to obtain an optimum data 
set, and then analysis including reduction to a l°xl° latitude/ longitude 
grid at common altitude datum (400 km), wavelength filtering, and reduction- 
to— the-pole to aid intepretations . 

ii) relate the magnetic anomaly map to prospective crustal sources, in 
terms of crustal thickness variation, depth to Curie-isotherm, and petro- 
logic variation between different geologic provinces. This can be done in 
conjunction with correlative geological and geophysical information such 
as regional gravity, aeromagnetics , seismic refraction to define crustal 
thickness, and available geological/structural data. 

iii) investigate the prospective relation between satellite magnetic ano- 
maly mapping, deep-seated geologic conditions and properties, and reconnais- 
sance exploration for and assessment of supplies of such resources as min- 
eral deposit provinces, energy fuels (oil, gas) in basins, and geothermal 
energy. 

iv) assemble and evaluate data on magnetic minerals (primarily titanomag- 
netite) whose presence and properties are the source of the magnetization 
of crustal bodies. This will help in modelling and interpreting magnetic 
anomalies as witnessed at satellite altitudes, in terms of crustal compo- 
sition, thickness, structure, and thermal regime. 

One may conclude, at this early stage of analysis and interpretation 
of Magsat satellite magnetic data, that the regional and global anomaly maps 





provide useful long-wavelength information on gross features of the Earth's 
crust. The satellite anomalies reflect, primarily, broad and deep-seated 
mid- to lower-crustal properties and conditions (magnetic mineral content 
of rocks, metamorphic state of basement rocks, crustal thickness, geother- 
mal gradient). As such, they portray major geologic structures and terranes, 
such as continental blocks sutured together, rather than the much shallower 
and more localized sources such as intrusions in the upper crust which yield 
intense but much shorter-wavelength anomalies on aeromagnetic or surface 
maps. 

In the U.S. central Midcontinent, the magnetic anomaly map reflects 
major geologic terranes and deep crustal petrology. It does not, at least 
at this stage of the investigation, indicate a magnetic signature for the 
MGA (Midcontinent Geophysical/Gravity Anomaly) , which is the gravity and 
aeromagnetic signature of the late-Precambrian-age CNARS (Central North 
American Rift System)— a major failed paleorift. This is for several rea- 
sons: the CNARS/MGA major portion which trends from Lake Superior down to 
central Kansas, 

i) is oriented along the north-south azimuth of the satellite orbit 
track, and is less likely to show an anomaly than an east-west trending 
structure or boundary. 

ii) the CNARS, although a long and deep structure, is relatively nar- 
row (60-150 km) for the resolution of satellite magnetic surveying, unless 
it is very anomalously magnetic; which leads to the following, 

iii) although the CNARS/MGA major portion is an immense seam of volcanic 
(basaltic) extrusives and intrusives in the upper- to middle-crust, and 
could be expected to have a large induced magnetization (due to suscepti- 
bility contrast) and an even larger remanent magnetization, its magnetic 
anomaly signature is modest in intensity at aeromagnetic heights and not 
apparent at satellite elevations. From geophysical modelling, this is be- 
cause there is the presence of both normally and reversely polarized re- 
manence in the rock complex due to a geomagnetic field reversal in Keweena- 
wan time. The resulting net magnetization is much less than it would be 
for, say, induced magnetization alone. 

For modelling of deep crustal sources, considerable progress can 
be made with the operating assumptions that, 

i) the main magnetization derives from the lower crustal rocks, es- 
pecially if having mafic composition and intermediate-grade metamorphism. 


iv 


ii) the dominant mineral special of concern is magnetite, with a 
Curie temperature (and corresponding isotherm as a lower bound to magnetiza* 
tion) of 560-57 5°C. 

ill) the magnetization is dominantly induced, such that the vector of 
magnetization is in the direction of the local geomagnetic field. 

It is recommended that work continue on analysing the existing 
Magsat anomaly data, on both a global and regional scale. Development and 
application of an optimum data-processing scheme, as for example including 
reduction-to-the-pole to spatially relate anomalies to sources, and cross- 
correlation with other digital data sets such as gravity and lower-crustal 
thickness, will continue to improve the resolution and use of the satellite 
magnetic data. For example, we have begun two new efforts to refine terres 
trial application of the Magsat data: application of the techniques used 

for our initial study, to interpreting the Indian subcontinent/Himalayan 
region; and study of satellite magnetics to help intepret recent volcanic/ 
geothermal regions such as hot spots (mantle plumes) and incipient rifts. 
The concept of using reduced-to— the— pole continental anomaly patterns (as 
UA5A workers are now doing) to try to reassemble Mesozoic Gondwanaland , or 
disassemble Precambrian plate collisions in cratonic regions, is an intri- 
guing application that holds great interest for the geoscience community. 
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INTRODUCTION 


The Satellite Mission 

NASA f s Magsat satellite was launched October 30 1979 from Vanden- 
berg Air Force Base in California (Langel et al, 1982a; and see Figure 1) . 

It was placed into a sun-synchronous orbit, that is, in the twilight plane 
(dusk meridian) to avoid the highly conductive (Sq) dayside current re- 
gion with its associated transient magnetic field effects* The orbit in- 
clination was 96.76°, initial apogee was 561 km, and initial perigee was 
352 km. As the mission proceeded, the orbit decayed to a more circular 
one (see Figure 2) , untiL after 6% months the apogee was under 400 km and 
the satellite came down on June 11 1980. 

Magnetic field data were collected from early November 1979 until 
the following May. On one complete orbital track, correction of the data 
for altitude effects will vary from, say at the early stage of the mission, 
from about 550 km to 350 km and back to 550 km again. The satellite loca- 
tion was known to about 30-60 meters radially from the Earth, and to about 
300 meters horizontally. 

The satellite had its magnetometers on the end of a 6-meter boom, 
to diminish the effects of spacecraft fields. The measuring apparatus in- 
cluded a cesium-vapor magnetometer for absolute (scalar) field measurements, 
and a triaxial set of flux-gate magnetometers for directional (vector) field 
measurements. The scalar instrument operated sporadically but was used to 
help calibrate the vector instruments. The vector field components, for 

A 

the "attitude-corrected" data, are accurate to better than 6 nanoteslas (nT) 
and the total (vector sum) field is accurate to better than 2 nT. 

The magnetic fields as measured by the satellite included contri- 
butions from the terrestrial core field (of magnitude up to about 50,000 nT) , 
and geomagnetic transient fields from external ionosphere/magnetosphere 
current systems (of magnitudes up to 100's of nT) , and magnetic anomaly fields 
from magnetization originating in the Earth's crust (of magnitudes *20 nT 
from an arbitrary zero level), at Magsat altitudes. 


nanoTesla in S.I. units; equivalent is gamma in cgs units 
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Figure 2 . Magsat orbit dimensions, for 1979-80 (from EOS, v. 61, p. 545, July 1980) 


Altitude Approximate 


Satellite 

Inclination 

idegs) 

Range 

ikon 

Dates 

Instrument 

Accuracy 
( nT) 

Coverage 

Sputnik 3 

65 

440 - 600 

5/5S- 6/58 

Fluxgates 

100 

USSR 

Vanguard 3 

33 

510-3750 

9/59- 12/59 

Proton 

10 

Near ground 
station* 

1963-38C 

Polar 

1100 

9/63- 1/74 

Fluxgate 

(1-axis) 

30-35 

Near ground 
station 

Cosmos 26 

49 

270 - 403 

3/64 

Proton 

Unknown 

Whole orbit 

Cosmos 49 

50 

261-488 

10/64-11/64 

Proton 

22 

Whole orbit 

1964-83C 

90 

1040-1089 

12/64- 6/65 

Rubidium 

22 

Near ground 
station 

OGO-2 

87 

413 -1510 

10/65- 9/67 

Rubidium 

6 

Whole orbit 

OGO-4 

86 

412-908 

7/67- 1/69 

Rubidium 

6 

Whole orbit 

OGO-6 

82 

397-1098 

6/69- 7/71 

Rubidium 

6 

Whole orbit 

Cosmos 321 

72 

270-403 

1/70- 3/70 

Cesium 

Unknown 

Whole orbit 

Azur 

103 

384 - 3145 
& 

11/69- 6/70 

Fluxgate 

(2-axis) 

Unknown 

Near ground 
station 

Triad 

Polar 

750- 832 

9/72 - present 

Fluxgate 

Unknown 

Near ground 
station 

““Near ground station” indicates no on-board 
sight of a station equipped to receive telemetry* 

recorder. Data were 

acquired only when the spacecraft was in 


Table 1 . Satellites that have measured the near-earth geomagnetic 
field (from Langel, 1980) . 



One of the mission objectives, and our particular project interest, 
has been to map and analyse variations in the crus tal-produced magnetic 
anomaly fields. This is to aid in assessing the utility of satellite mag- 
netic data in interpreting crustal geological character and geophysical 
properties, such as major geologic structures, nature and evolution of con- 
tinental tectonic regimes, crustal composition (rock type) , prospective 
localization of economic resources (mineral deposit provinces, oil and gas 
basins, geothermal regions), and variations of crustal thickness. 

Magsat was the first satellite to provide a global survey of vector 
magnetic fields and, because of its lower -altitude orbit, provides the 
most accurate and highest-resolutlon survey of the geomagnetic field and 
anomalies to date. 

Other satellites have measured the near-earth's total (scalar) 
field, but to a lower accuracy. They are listed in Table 1 (Langel, 1979 
and 1980) . None of these missions were designed for solid-earth magnetic 
field studies. Only the TOGO (Polar "Orbiting Geophysical Observatory") 
satellites of 1965-71 provided fairly accurate global surveys, and of the 
earth's total-field only. Those had satellite altitudes varying from 400 
to 1500 km. 

It was of considerable interest that even with the POGO data set, 
one could extract a component of field due to crustal magnetic sources. 

That is, one could begin to map crustal anomalies (Regan et al, 1975). 

Table 2 lists a comparison of crustal anomaly mapping using the 
POGO data, and that from the more accurate, higher-resolution Magsat mis- 


Table 2 . Comparison of crustal magnetic anomaly 
mapping by satellite. 


POGO (1965-71) 

- scalar field measurements 

- most of the data from altitude 
above 500 km. 

- spatial resolution about 500 km 

- observed anomalies up to irlO nT 


MAGSAT (1979-80) 

- both scalar and vector component 
field measurements 

- data from about 300-500 km. alti- 
tude 

- spatial resolution about 200-250 km 

- observed anomalies up to i20 nT 






An aspect of the vector component measurements is that they iden- 
tify fields not parallel to the dominant main core field, so one has the 
prospect for interpreting the possible presence of remanent magnetization 
in large cirustal sources. That is, one could separate the induced and any 
remanent component of magnetization. If would be very instructive to be 
able to identify and interpret the presence of crustal remanence from sat- 
ellite data, if this proves to be possible (Galliher and Mayhew, 1982; 
Bhattacharyya, 1977) . 

Figure 3 shows how, in theory, the magnetic anomalies from two 
localized crustal sources should look at heights of 300, 450, and 600 km, 
assuming induced magnetization in a horizontal Earth's field. There should 
be a substantial increase in both anomaly amplitude and resolution, as one 
descends to the lower (satellite) altitude. 

The Data 

The satellite measured vector fields sixteen times a second and 
scalar fields (when operable) eight times a second. The data were pre-pro- 
cessed by NASA/Goddard Space Flight Center, and distributed on 

CHRONFIN magnetic tapes - calibrated data in a North, East, and Z (down' 
ward) system 
or 

INV-B ("Investigator") magnetic tapes - data decimated to one sample 
each 5 seconds of time (i.e. one-fortieth of scalar data, one- 
eightieth of vector data) ; this corresponds to a point about every 
36 km of orbit track. 

The Investigator tapes are the preferred product for workers interested in 
assembling magnetic anomaly maps for regional crustal interpretation. Those 
tapes also include accessory information such as data/sample location in 
latitude, longitude, and radius, predicted geomagnetic field from a spheri- 
cal harmonic core field model, local time, dip latitude, and an estimate of 
the external field. For the time-invariant crustal anomaly field, one would 
like data from magnetic "quiet" times, when there were low temporal magnetic 

variations. Parameters indicating this are the geomagnetic indices K and 
* P 

D gt » ant * Figure 4 shows how these varied over the 7-month lifetime of Magsat 

* 

e.g. D st index measures temporal variation of magnetic fields produced by 
the equatorial ring currents (about 10^ amperes, at about 5-8 earth radii 
and near the geomagnetic equator), relative to values on "quiet" days. 
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ORIGINAL PAGE iW 
OF POOR QUALITY. 


Days recommended by 
Japanese Team 


GEOMAGNETIC INDICES FOR MAGSAT MISSION 


▼ Very Quiet Days 


(DAILY AVERAGES) V Quiet Days 



NOV 

1979 


DEC 


JAN 





On Figure 4 are indicated the best ’’quiet-time" days, as selected by the 
Japanese project team. 

Another time-varying component contributing to measured magnetic 
fields could be that due to currents induced in the earth by the ionos- 
pheric currents. Lateral variations in electrical conductivity in the 
crust and upper mantle will give current distributions which could produce 
magnetic field patterns which might be fairly persistent in time. These 
are due to transient electromagnetic induction, and not true crustal mag- 
netization. Their presence would be minimized by appropriate time-averaging 
and comparison of the data, particularly for magnetic "quiet- time" data 
being used. 

The data were to be pre-processed by NASA at two levels of accuracy: 
"intermediate-a.ttitude"-corrected, to 20 arc-minute accuracy, and "fine- 
attitude’’-corrected, to 20 arc-second accuracy. 

Spherical harmonic modelling of the magnetic field seen by Magsat, 
by Langel and others, yields a power spectrum of components that suggests 
that the internal core field has terms up to about order/degree 15, and 
that the crustal (anomaly) fields have terms from order/degree 13 and up . 
That is, the field originating from the crust dominates at order ("n") 15 
and above. In processing Magsat data to get anomaly data, many workers 
typically thus take out a field model of order/degree 13. 

After considerable data processing, NASA/Goddard provided global 
magnetic anomaly maps (Figures 5 to 8 , Langel et al, 1982b). These were 
generally compiled from quiet-time data, and plotted on a global grid of 
2°x2° latitude by longitude (i.e. about 220 x 220 km). 

These maps were compiled as of November 1981. The zero levels are 
arbitrary, particularly for a given local geographic region, and in fact 
may vary in a relative way across the global map. The scalar anomaly map 
(Figure 8) is for *50° latitude, has an average data altitude of 404 km, 
and has about 25 points per 2°x2° block used for the map assembly. For 
this map, the data range is from +18 nT to -22 nT. There is some apparent 
east-west banding, which can be non-crustal in origin. This could be intro- 
duced by, 

i) processing of data, e.g. by the along-track (north-south) filter- 
ing of data used to reduce the north-south striping originating 
from satellite passes having different altitudes and the result- 
ing data having different arbitrary zero-levels 
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Ii) the satellite tracks are north-south, which will tend to em- 
phasize east-west features that are crossed in transit 
iii) incomplete removal of the field effects of the equatorial ring 
current, near the magnetic equator. 

Thus east-west anomalies are emphasized and the north-south anomalies are 
relatively suppressed. This has the effect of reducing the detectability 
of north-south crustal features. 

The vector-component anomaly maps have field variations of up to 
110-15 nT. The magnitudes are, 

NASA anomaly map field variations 

AX (positive northward) +11 to -14 nT 

AY (positive eastward) +10 to -10 

AZ (positive downward) +15 to -14 


For these maps, the latitude range given is 160°. There were about 

eight points per 2°x2° block, and there was an average altitude of about 

430 km. In averaging the data, no adjustment was made for the different 

altitudes of the data, which varied from 340 to 552 km. It can be noted 

* 

that the delta-X anomaly map has more apparent east-west anomalies be- 
cause of the along-track filtering in the north-south direction. The delta- 
Y map is noisier. This is presumably due to the influence of fields from 
ionospheric currents; for example, the equatorial ring current ll5° from 
the magnetic equator, and field-aligned sheet currents in the auroral zones 
(about 10-15° from the geomagnetic poles) . 


This Project 

Our objective has been to use magnetic field data acquired by the 
Magsat satellite to construct a magnetic anomaly map of the II. S. central 
midcontinent. The data set will be the best that is practicable; for 
example, selected for reliability, and averaged on a l°xl° latitude/ longi- 
tude grid spacing rather than the customary 2°x2° spacing used for POGO 
studies and NASA's Magsat maps. 

The magnetic anomaly map can then be used to help interpret (see p.l3bj 
excerpt from Langel 1980) gross crustal character (geologic and tectonic 
terranes, crustal thickness and petrologic composition), and assessed as 
an aid in understanding the nature of the large-scale structural feature 


There can be an overlap of effects between east-west anomalies due to 
crustal sources being relatively emphasized , and apparent east-west 
anomalies being introduced as an artifact of the data processing 
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in the area— -the Central North American 
Rift System (CNARS) , especially the 
major portion of this which is expres- 
sed as the Midcontinent Geophysical/ 
Gravity Anomaly (MGA) . 

The MGA/ CNARS is a major crustal 
feature that extends over 1000 km from 
Lake Superior through Minnesota, Iowa, 
and to E’. -central Kansas. It repre- 
sents an extrusion and intrusion of ma- 
fic rocks in late-Precambrian time (Ke- 
weenawan age, about 1100 million years 


In addition to the continental-scale studies, 
several investigators will study more limited areas 
or particular tectonic features. Godivier 
(ORSTOM) will extend the work of Regan and 
Marsh 24 in the region around the Central African 
Republic where ORSTOM has a large amount of 
correlative data. Gasparini (Osservatorio Vesu- 
viano) will investigate the Curie depth and 
volcanism in the Mediterranean area. Won (North 
Carolina State University) will study a combination 
of Magsat, aeromagnetic, and regional gravity data 
in the eastern Piedmont of the United States. Car- 
michael and associates (University of Iowa) will 
study the central midcominent of the United States 
with particular attention to the known midcomi- 
nent geophysical anomaly. 

(from Langel, 1980) 


ago), with crustal downwarping and then reactivation to form a central horst 
bounded by faults which have up to 10,000 meters of vertical displacement 
on them. The most prominent central portion of the MGA is in Iowa, where 


the central basaltic horst and flanking sedimentary troughs produce the most 
pronounced gravity anomaly in North America — over 160 milligals, in a lateral 
distance of less than 40 km. The MGA trend is the major part of the CNARS 


system, and represents a failed paleorift zone. Although the CNARS and its 


Precambrian basement rocks are buried over much of its extent in the central 


midcontinent, where the analogous and probably genetically-related terrain 
is exposed to the north end (northern Minnesota and Wisconsin) there are 
important economic mineral deposits — iron, zinc, lead, copper. 

The regional, long-wavelength satellite magnetic data can be corre- 
lated with other geophysical and geologic data sets to help interpret the 
deeper and broader features and geologic setting. This thus adds new in- 
formation for the geologic, tectonic, and structural framework of the region. 

The data analysis techniques used should help in determining and re- 
fining optimum methods for processing the satellite magnetic data, for fur- 
ther application in interpreting magnetic anomalies for other geographic 
regions and geologic/crustal terranes. 

Further on this project, an objective has been to compile and eval- 
uate magnetic property data of crustal rocks of importance in modelling mag- 
netic anomalies. This would include consideration of geothermal gradients 
and depth to the Curie-temperature (T c ) isotherm, at which magnetization 
is lost. 


Activity Schedule 


1980 

June 11 - Magsat satellite ends mission 

Dec. 3 - Project contract starts 

3-6- Magsat meeting at Goddard Space Flight Genter 

1981 

Mar. 30 - we received sample data tape (limited profile data) 

July 10 - received first intermediate-attitude-corrected tape with 

profile data (for first 80 days of mission); i.e. first tape 
of our Investigator data set 

July 30-Aug. 14 - NASA Magsat meeting and IAGA Scientific Assembly, 
Edinburgh; and Internal. Conference on Tectonics, Oslo 

Aug. 19 - received second tape of Investigator data set 

Sept. 17-18 - Ann. Midwest Amer. Geophys. Un. meeting, Minneapolis 

Nov. 12 - received Investigator "final output" data tape, for quiet 
time data for period Jan. 19-May 19 1980 

1982 

Oct. 17-21 - Ann. Soc. of Explor. Geophysicists meeting and Magsat 
session, Dallas 

Dec- - we requested a "no-cost" extension of contract, since Principal 
Investigator would be absent on a Leave assignment overseas 
for 5 months (as geoscientist on a college ship) 


1983 


Jan: -May - P.I. on Leave overseas 
Sept. 3 - contract expires 


DATA PROCESSING 


In order to investigate the utility of satellite magnetic data 
for geologic interpretation, and resolution of causative crustal sources, 
we need an optimum data set. This would incorporate 

i) quiet- time data, from time intervals when the transient (external) 
field variations were a minimum 

ii) data "cleaned" of spurious information introduced in the data- 

acquisition or pre-processing operations, i.e. that is non-terres- 
trially based 

f iii) comparison and correlation of data from spatially nearly-coincident 

t orbit tracks, to best assure a reliable and consistent data set 

l iv) data "reduced" to a common elevation and grid, so that data are 

equalized and distributed and maps can be contoured to reveal mag- 
netic anomalies. 

The objective is to remove contributions to the data set of time-varying 
fields (due to geomagnetic transients) and spatially-varying field values 
(due to differences in satellite altitude) , and artifacts of the pre-pro- 
cessing of the raw data. This will yield the invariant field "morphology" 
which is due to crustal magnetization, and hence to the best resolution 
that can be achieved for the Magsat satellite orbital altitudes. 

Correction of Data 

A 

The magnetic field profiles on individual orbital tracks have 
been examined to identify "bad" data points or profiles. The isolated 
spikes of the former can have magnitudes of tens of 1000' s of nanoteslas. 
Data problems can include 

i) individual spikes of erroneous data 
ii) raised or lowered effective "zero level" of a profile, due to 

— hardware problems 

— pre-processing, before Investigator tape sets are produced 

— long-period (and hence long-wavelength) diurnal or geomagnetic- 
transient field effects on a particular orbit pass 


the "Investigator" tape data are decimated, so data points are about 
36 km apart 

data processing here was performed on a Perkin-Elmer 3220 (Iowa Geolo- 
gical Survey) and IBM 370 (University of Iowa) 


Changing the zero level of (part of) a profile, due to long-period field 
variations, would lead to striping of the anomaly map along the azimuth 
of the satellite tracks. This is roughly north-south. Such striping 
would obscure and confuse the identification and interpretation of crustal 
sources (structure, and features such as geologic province boundaries) that 
strike in a north-south direction. 

Figure 9(a) shows an example of a pass (no. 85) over the U.S. mid- 
continent, from 25-50° N. latitude. The individual erroneous data spikes 
have magnitudes up to 57,000 nT, superimposed on the true anomalies with 
magnitudes of a couple of tens of nanoteslas. The true anomaly is thus 
rendered indistinguishable by comparison. When these spikes are removed, 
as seen in Figure 9(b) , the proper anomaly field is revealed. It has a 

fc 

total amplitude here of about 24 nT . 

Calculation of scalar (total-field) anomaly data 

Since the total-field magnetometer did not perform continuous data 

gathering, the scalar magnetic field (H^) magnitude can be calculated from 

the vector summation of the X (or E. , for northward) , Y (or H , for east- 

x y 

ward), and Z (or H z , for downward) component fields as measured by the vec- 
tor magnetometers. 

For example. Figure 10 shows the "raw” (NASA pre-processed) inves- 
tigator data on pass 85, having a track profile over the U.S. midcontinent. 

Figure 10 (a) , (b) , and (c) show the H , H , and H components respectively 

y z 

of the "anomaly" field, with the calculated scalar total-field being 
shown in Figure 10(d). At these relatively northerly magnetic latitudes, 
the H z component and total-field H^, are similar; near the magnetic equator, 
the H x component and would be similar. 

As is generally the case, at least at these magnetic mid-latitudes, 
the component-field profile data are "noisier" than the scalar magnitude 
profile. 


The zero-level datum on the data tape, and for contoured anomaly maps 
(for NASA's, and here) is arbitrary. In fact, because of processing 
corrections and techniques, the relative zero level can vary across a 
geographic region which is sufficiently large. In analyzing and inter- 
preting profile or map anomalies, the extrema terms "high and low", 
"maximum and minimum", and "positive and negative" are fairly synonymous. 
A magnetic "negative" could be, with a lowered arbitrary zero level for 
the local region, actually be a relatively low-amplitude "positive". 
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Magnitude (in 10 2 /nT) 
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The short-period ’’noise" variations, recorded as a function of 
time, are displayed as ' ' shcr t-wavel eng th" variations when plotted as a 
function of profile (track) distance. They have amplitudes of several 
nano teslas. They are clearly temporal variations, since 

i) they are different for each pass over the same or near-neighbor 
tracks 

ii) they are too short to be produced at satellite elevation by 
true sources in the earth’s crust. 

These variations are thus not terrestrially-based, and need to be removed 
from the desired time-invariant crustal anomaly data set. This can be 
done by along-track filtering, or smoothing by a running average, of either 
the noisy component-field data or (preferably) of the calculated scalar 
magnitude data. 

Correlation of tracks and passes 

To develop an optimally reliable data set for the mapping of mag- 
netic anomalies, we need to compare profile data for the same satellite 
track but at different times (i.e. different passes), for nearly-coinci- 
dent tracks, and for tracks that cross. This comparison will allow recog- 
nition of time-varying magnetic-field effects that could masquerade as an 
"anomaly" on an individual profile. These will include the long-period 
field changes which would change the apparent "zero" base level for pro- 
files which would contribute to striping in the north-south direction when 
profile data are averaged. 

To illustrate the time-varying magnetic effects left in the "ano- 
maly" profile data that are supplied. Figure 11 shows H component data 

X 

for passes 254 (Figure 11(a)) and 531 (Figure 11(b)). These can be com- 
pared with Figure 10(a) for pass 85 over the same track. In terms of the 
geomagnetic-transients and base-level changes, the comparison shows, 


pass 

approx. 

"base 

relative 

level" 

approx * short-period 
deviations from mean curve 

85 

- 6 

nT 

± 10 nT 

254 

0 

(arbitrary 

+14 to -10 

531 

- 4 

reference) 

± 6 


Similarly, Figure 12 shows H component data for passes 254 and 
531, and can be compared to Figure 10(b) for pass 85. 
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Figure 11 . component field data at different times over 

same track. Field (nT) vs. latitude (degrees North) 
(a) Pass 254 
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A comparison of as recorded and pre-processed prior to dis- 
tribution, for the three passes over the same track (see Figs. 10(b) and 
12(a) and (b)), gives. 


approx, relative 
"base level" 


approx, short-period 
deviations from mean curve 


+3 to —1 3 
+3 to -6 


Smoothing of track (profile) data 

Data along individual profiles can be statistically smoothed, to 
reduce the effect of geomagnetic transient disturbances. Figure 13 shows 
the effect of a running-average smoothing of the previous pass 85 (H^ 
scalar data, as seen in Figure 10(d)). 

Reduction of data to grid and datum level 



:o; 









As seen, in Figure 14, each data prism is l°xl in latitude/ longi- 
tude (i.e. about 110 km x 110 km), and varies in height from 350 up to 
560 km . The data will be interpolated and weighted-averaged vertically 
to a common datum altitude of 400 km, and laterally to the center of the 
prism, to form the net average data value (point) for that prism. 

The variation in measured anomaly field values, and the averaging 
that will have to be done to produce a net data point, is illustrated 
by a track over the central midcontinent that has the three nearly-coinci- 
dent passes 85, 254, and 531. This satellite track is shown in Figure 15 
(for pass 85) , and the corresponding orbit (height) is shown in Figure 16 . 
Again, the latitude ranges from 25°N. (north-central Mexico) to 50°N. 
(north of North Dakota). The radius for this varies from about 6765 to 
6735 km; subtracting a mean Earth's radius of 6371 km, the altitude varies 
from 394 down to 364 km. Comparison of the scalar (H^,) field for passes 
254 and 531 is shown in Figure 17 . This can be compared to Figure 9(b) 
for pass 85. 

On this track, the large positive magnetic high is located over 
North Texas /Oklahoma. The anomaly amplitude (maximum relief) varies with 
the altitude of the satellite on a particular pass. This is summarized 


below. 


anomaly amplitude 


mean altitude of pass 



The anomaly amplitude falls off with increased altitude, as expected 
Plotting magnetic field anomaly maps 


The processed selected data can now be gridded, plotted, and con- 
toured to display the magnetic anomaly field. In the maps to follow, the 
grid is l°xl° in latitude/ longitude, the data, are interpolated to a 400 
km altitude datum, and the map projection is Albers equal-area. The ac- 
tual area of processed data is about 50% larger than the study area shown 
here, to include a beltway surrounding the data area and thus diminish 
"edge effects" in the subsequent data processing to follow. 


Recall that on one-half a satellite orbit, the altitude can vary as 
much as (e.g. near the start of the mission) from about 550 km down 
to about 350 km. 
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The study area is the U.S. midcontinent — see Figure 18(a) — extend- 
ing from 80° to 105° W. longitude, and 25° to 50° N. latitude (i.e. about 
2800 km) . The contour interval is 2 nT. 

The and component magnetic-field anomaly data, mapped in 
Figures 18(b) and (c) , are fairly irregular and exhibit apparent striping 
of the anomalies along the azimuth of track directions. This is due to 
incomplete equalization of the zero-levels for multiple passes. Simple 
wavelength filtering did not improve the appearance significantly. The 
and (calculated) scalar H^, data are more regular, and are also similar 
in appearance because of the fairly high magnetic latitudes of the study 
area. The H anomaly map is shown in Figure 18(d), and the EL, map in 
Figure 19 . The geographic latitudes of 35-50° translate, because of the 
proximity to the geomagnetic pole (about 79° N. latitude) to about 45-60° 
magnetic latitude ♦ For the simple equation for the magnetic-field incli- 
nation I at geomagnetic latitude 
tan 1=2 tan'O- 

the inclination is thus quite high, 64-74°. For a crustal induced magne- 
tization, the vertical-component and total-field anomalies should be simi- 
lar in appearance. 

Wavelength filtering 

The and maps were less severely striped along the orbit track 
direction, and wavelength filtering (two-dimensional) removes the apparent 

icie 

effect of this. 

Figure 20 shows the scalar (H^, or total-field magnitude) anomaly 
map, after applying a 4° high-pass wavelength filter (i.e. passing wave- 
lengths greater than 4°, or roughly 400 km). The anomaly values range 
from highs of 16 nT (over Teimessee/Kentucky) and 10 nT (over north Texas) 
to a modest high extending from the latter up to the northeast toward Lake 
Michigan, to lows of -4 nT (over S. Dakota/Nebraska) and -16 nT in the 
Gulf of Mexico with a diminishing low extending up the Mississippi Embay- 
ment. The anomaly values are relative to the arbitrary zero-level chosen 

A 

applying a high-pass wavelength filter (i.e. a low-pass wavenumber filter) 

C * 

We should note that for our purposes — interpreting crustal geologic struc- 
ture and composition— we would like to do as little high-pass filtering 
as possible. This is so as to retain the short-wavelength anomaly fea- 
tures associated with crustal sources, and achieve the best resolution 
of these that is possible. 
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Figure 18 . Magsat vector component anomaly maps for U.S. mid 
continent. 1 :cl latitude/longitude grid, data 
weighted-averaged to altitude of 400 km. Albers 
projection. (Black, 1901) 

(a) Base map. Scale: 5° of latitude is about 550 km. 




















ORIGINAL PAGE J3 
OF POOR QUALITY 



Figure 19 . Magsat scalar (total-field) anomaly map, for 
U.S. midcontinent. l°xl° latitude/longitude 
grid, data weighted-averaged to altitude of 
400 km. Albers projection, contour interval 
of 2 nT. (Black, 1981) 









for the area. However, the total anomaly range over the midcontinent 
(excluding the Gulf of Mexico) is about 24 nT. 

Reduction to the magnetic pole 

To relate magnetic anomalies (highs and lows, and gradients) to 
their causative crustal sources, it is helpful to reduce the map "to the 
pole". This compensates for the variation of magnetic inclination with 
latitude (from zero at the magnetic equator, to ^90° at the magnetic 
poles) , by effectively transposing the area to the north magnetic pole 
where the magnetic field is vertically downward. There would thus be 
radial magnetic polarization, and the anomalies are rendered more symme- 
trical. Implicit in this process is the assumption that the magnetiza- 
tion direction is in the direction of the local Earth's field; that is, 
that it is all induced magnetization. 

Figure 21 shows the scalar anomaly map of Figure 20, reduced to 
the pole. This has shifted positive anomaly features to the north by one 
to two degrees of latitude, and reduced anomaly asymmetry. The altera- 
tion of Figure 20 is not dramatic, however, because of the relatively 
high magnetic latitude of the study area. The major magnetic highs have 
been shifted northward: the north Texas one up to the Oklahoma border, 

and the Tennessee one up into Kentucky. 

The agreement between the anomaly map and the reduced-to-the- 
pole map supports the general assumption that, at least on a large (i.e. 
long-wavelength) scale and as seen at satellite elevations, it is induced 
crustal magnetization which is responsible for the major anomalies. 

It is of interest to consider how data treatment and processing 
as done here could affect, or improve, the resolution of the satellite 
magnetic mapping. Figure 22 shows a comparison for the midcontinent of 
a Magsat scalar anomaly map, as excerpted from a preliminary global data 
set (NASA, March 1981) with our filtered, reduced to the pole anomaly map. 

There is more detail in our anomaly map (Figure 22(b)), and it is 
useful to try to assess how far the resolution of the satellite data can 
be "pushed" by processing to reveal crustal structure and properties. 


* 


Or, possibly, induced magnetization plus a viscous remanent magnetiza- 
tion (VRM) acquired in the direction of the present earth’s field, under 
the influence of elevated temperature at depth, and long time. 


Figure 21 . Scalar anomaly map, l°xl° data at 400 km altitude, high-pass 
filtered and reduced to the magnetic pole. Contour interval 
2 nT. (Black, 1981) 













The greater detail here (see Figure 22 again) is due to 

i) the care taken in removing "bad" data (isolated points, and 
tracks), and comparing profiles (e.g. of near-neighbor tracks) 
to optimize the consistency and reliability of the final data 
set 

ii) treating the data on a finer l°xl° (110 x 110 km) grid, rather 
than the customary but coarser 2°x2° (220 x 220 km) grid spacing 
used for many global-scale maps 

iii) for a given data prism, doing a weighted-averaging of the data 
laterally to the central point and vertically to the datum alti- 
tude of 400 km. The data closest to the averaged point are thus 
favored over more extreme values. In contrast, some other compi- 
lations of data do a simple averaging of all data in a prism ele- 
ment to arrive at an average value, without regard to the place- 
ment of individual points in the prism or the distribution of the 
whole set of points used. 

Figure 23 is a foldout map, repeating Figures 19 (scalar H^) and 20 (fil- 
tered H^,) together for comparison. 
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INTERPRETATION OF CENTRAL MIDCONTINENT REGION 


Geologic Setting 

Figure 24 shows some major tectonic and geophysical features in 
the larger U.S. midcontinent region. Figure 25 shows our particular study 
area, the central midcontinent from the Dakotas and Minnesota in the north 
to Oklahoma and Arkansas in the south. Figure 26 portrays the middle re- 
gion of the study area, featuring the major arm of the arrested CNARS 
(Central North American Rift System). 

The basement rock of the region consists of a Precambrian craton, 
exposed in northern Minnesota, northern Wisconsin, and northwestern Mich- 
igan and farther north as the Canadian Shield. To the south, it is covered 
by a veneer of Paleozoic and younger sediments. Major tectonic and geo- 
physical features include the following, several of which are aulacogen 
(i.e. failed rift arm) structures (e.g. Keller et al, 1983): 

i) the Central North American Rift System, which extends from Lake 
Superior in a major arm over 1000 km southwesterly through Minnesota, Iowa, 
and eastern Kansas, where it is subparallel to the Nemaha uplift. The 
CNARS has another arm extending southeasterly from Lake Superior down 
through Michigan. 

ii) the Great Lakes Tectonic Zone, running west-east through Minne- 
sota from South Dakota, and then along the south side of Lake Superior 
where it is subparallel to the "Great Lakes magnetic lineament" . It is 
a boundary between (Superior province) Archean granite and greenstone to 
the north, of age 2500-3100 million years, and Archean gneiss in central 
and southern Minnesota to the south, of age 2500-3800 million years. 

iii) the Southern Oklahoma aulacogen, to the south. This is repre- 
sented by the Anadarko basin (western Oklahoma) and Ardmore basin (SE 
Oklahoma) , these being flanked on the southwest by the Amarillo/Wichita 


Presently still an informal term, but used for example in Hilderibrand 
et al (1983). The term "magnetic (mega) lineament" implies a raore-or- 
less linear alignment of anomalies and gradients, that extends for at 
least several hundred kilometers, and most likely represents a major 
structural and lithologic boundary or feature. In the U.S. midcontinent, 
the magnetic lineaments noted here would be Keweenawan-age (1100 million 
years) or older, and these zones of crustal weakness would help control 
the placement and extent of Keweenawan rifting and intrusion. 
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Figure 26. Region of main arm of CNARS, as outlined by 
MGA (Midcontinent Geophysical Anomaly) . 



uplift. This aulacogen is attributed to late Precambrian/ early Paleozoic 
rifting that formed the proto-Atlantic, along the "rifted continental mar- 
gin" noted on Figure 24. It is this margin where the Appalachian and Oua- 
chita mountains, formed in a later continental collision, are now located. 

iv) the Mississippi Embayment, with associated aulacogen features 
such as the Reelf oot rift, which also originated with the late Precambrian/ 
early Paleozoic rifting to its south. The rift was then reactivated in 
the Mesozoic, forming the Embayment over it as areal subsidence in pri- 
marily Cretaceous time. The latter development accompanied the (second) 
opening of the Atlantic. The Reelf oot rift and associated structures is 
presently active seismically, especially at the northern end. The rift 
has a central downdrop (graben) of 1.5 to 3 km. 

v) the enigmatic "South-central magnetic lineament" (see footnote 
on page 45), a belt of magnetic anomalies and trends/alignments, and grav- 
ity lows (about -30 milligals) , that extends from eastern Tennessee more 
or less linearly to the northwest. It begins at the "New York-Alabama 
lineament", that in turn extends SW-NE from Alabama up to central Pennsyl- 
vania/New York, and passes through the north end of the Reelf oot rift in 
the New Madrid region. It continues on to the dislocation of the MGA in 
southeast Nebraska (see Figure 24 again) (Arvids on et al, 1982; Zietz and 
Hatcher, 1983). Some suggest it may continue, although as a more obscure 
anomaly trend, to the northwest as far as Wyoming. There is some specu- 
lation that this apparent geophysical trend may be reflective of a trans- 
form fault, associated with late-Precambrian plate motions. At the north 
corner of the intersection of the South-central magnetic lineament and 
the New York-Alabama lineament, in central Kentucky/Tennessee, there is 
a region of large magnetic highs. The basement rocks are mafic volcanics 
and metamorphics . The former are similar to the Keweenawan volcanic rocks 
of Lake Superior, and may be of the same age (Hildenbrand et al, 1983). 

Aulacogens are typically directed into a continental craton from 
a rifted margin, and these failed (paleo)rift arms represent once-tensional 
flaws in the crust which can be later reactivated with sufficient tension 
or compression. The episodes of rifting, either arrested or successful, 
in the U.S. midcontinent (Keller et al, 1983), are listed on the next page. 


Tectonic episode 


Associated rift features 


Geologic time 

Keweenawan (1100 x 
10 6 yrs ago) 


precursor to opening of 
proto-Atlantic 


CNARS/MGA structure 


late-Precambrian/ opening of proto- 
early Paleozoic Atlantic 

(500-800 x 10 6 
years ago) 


Reelfoot rift; Southern 
Oklahoma aulacogen 


then closing of proto-Atlantic 


early Mesozoic opening of Atlantic 

(150-200 x 10 6 
years ago) 


Mississippi Erabaymentj^ 
Rio Grande aulacogen 


Rifts are often associated with the formation and development of 
sedimentary basins. For example, collapse (subsidence) of incipient or 
failed rift systems may be related to the development of the following 
basins in the U.S. midcontinent (see Figure 25), 


Tectonic development 

mid-Michigan (Keweenawan) intru- 
sion (eastern arm of CNARS) 

Midcontinent Geophysical Anomaly 
(MGA) intrusion (western main 
arm of CNARS) 

Southern Oklahoma aulacogen 

and possibly 

MGA/Nemaha uplif t/granitic plu- 
tonism in NE Kansas and NW 
Missouri 


Sedimentary basins 
Michigan Basin 

deep narrow basins flanking a 
central horst, in southern 
Minnesota and all of Iowa 

Anadarko and Ardmore Basins 
Forest City Basin 


Rifting, arrested rifting, reactivation of a paleorift, and deep struc- 
tural lineaments, can help control subsidence, uplift and erosion, and 
structural/ sedimentary history. This can have important implications 
for exploration for hydrocarbons. Rifting itself can also have a direct 
consequence for invasion by hydrothermal fluids, or magmatic intrusion, 
or other emplacement mechanisms for mineral deposits. It is for under- 
standing of these deep structures and major crustal features, that sat- 
ellite magnetic data have promise. 

The CNARS complex, or Midcontinent Rift System, is buried over most 
of its length. It was initially identified, and its extent mapped and 


Not to be confused with the Rio Grande Rift farther west in New Mexico. 
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modelled, using its characteristic gravity and magnetic anomalies ex- 
tending from Lake Superior down to Kansas. Since its gravity signature 
was the first known expression of the CNARS, the structure is sometimes 
referred to as the Midcontinent Gravity (or Geophysical) Anomaly, or MGA. 

The MGA arm of the CNARS is a failed paleorift zone of Keweenawan 
(late Precambrian, 1100 million years ago) age. It formed with attempted 
rifting of the North American craton, with an elongate seam of mafic in- 
trusives and layered extrusive basaltic rocks, up to 80 km wide. Asso- 
ciated basinal subsidence resulted (in the central portion, in southern 
Minnesota and Iowa) in deep flanking basinal troughs of downwarped sedi- 
ments. Reactivation there of the central basaltic seam has left it ele- 
vated as a fault-bounded horst with a vertical throw (uplift displacement) 
of up to 10,000 meters. This is in the major central portion of the MGA, 
in Iowa. 

The MGA's gravity expression of the CNARS is shown in Figure 27 , 
for the Bouguer gravity anomaly, and its aeromagnetic anomaly signature 
is shown in Figure 28 . The anomalous trend, outlining the buried struc- 
ture (central mafic core and flanking b as inn} , has a length of over 1000 
km and total width of up to about 150 km. 

The gravity anomaly over the most pronounced portion of the MGA 
has a central positive zone that exceeds +60 milligals in value, and 
flanking lows (over the adjacent sedimentary troughs) less than -100 milli- 
gals. The maximum difference, of over 160 milligals (on a profile in cen- 
tral Iowa less than 40 km long) is the largest such gravity gradient in 
North America. This indicates the very long, deep, and pronounced litho- 
logic inhomogeneity in the crust here. 

The aeromagnetic expression of the MGA is a belt of anomalies with 
a well-defined northeasterly trend, but having fairly modest amplitudes. 

Early studies of the MGA (e.g. Chase and Gilmer, 1973; Oc.ola and 
Meyer, 1973) focussed on the northern fragment of the rift system in nor- 
thern Minnesota and Wisconsin, where the basement rock outcrops. The 
central (southern Minnesota, and Iowa) and southwest portions could not 
be studied in detail because of burial under sediments, and a general lack 
of comprehensive geophysical information such as seismic surveys or gra- 
vity and magnetic mapping and modelling. Further, the structural develop- 
ment (especially as regards fault-bounding of the central basaltic seam, 


51 



Figure 27. Bouguer gravity anomaly (i.e. MGA) of CNARS paleorift in 
central midcontinent (map excerpt from Woollard and Joe- 
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and subsequent uplifting to form a great horst) is apparently different 
north of the Great Lakes Tectonic Zone (see Figure 25) than it is south 
of there, through central Minnesota and Iowa. 

There is now more understanding (e.g. Green, 1983; and here by 
Anderson and Black, 1981 and 1983) including detailed magnetic and gra- 
vity modelling to be described in a later Section. The complete CNARS 
extends in two long sinuous arms, for a total length of 2200 km from 
east-central Kansas northward to Lake Superior, then eastward through 
that basinal lake and down through the Michigan Basin into Ohio. The 
major western arm, or MGA, has a genesis and structure which most likely 
varies from north to south. In general, the northerly segment began its 
attempt to rift the continental craton about 1150-1200 million years ago, 
with fissure eruption of basaltic lavas onto a broad plain forming a trend 
of partially-overlapping flood plateaus. Extrusion was accompanied by 
progressive subsidence in the central region, perhaps because of regional 
tension or isostatic adjustment to the loading. The successive flows 
erupted over a period of 20-30 million years (Green, 1983) around 1100 
million years ago, producing plateau thicknesses of up to 7000 meters. 

With crustal subsidence, there was accumulation of thick sediments, which 
subsequently eroded with uplift of the central mafic core along steep boun 
ding faults. This later reactivation (completed by late-Precambrian time) 
was not everywhere present in this northern segment, and is attributed to 
regional crustal compression. While the total volume of flood basalts 
is estimated to be comparable (over 300,000 km ) to volumes in other 
great flood basalt provinces (Deccan/India, Parana/Brazil, Columbia River/ 
Washington) , there are major differences in character between this paleo- 
rift event and structure, and the ’’type" modern rift — the East African 
rift. This is for the MGA as examined in outcrop in the Lake Superior 
region. The differences include: 

CNARS/MGA failed paleorift East African rift 


plateau basalt rift plateau basalt rift 


plateau basalt rift 
but 


broad subsidence (downwarp) 
with the extrusion and rift 
attempt 

tholeiitic flood basalts, and 
basaltic andesites 
large positive gravity anomaly 


graben structure with 
rifting 

alkalic basalts, and central 
volcanoes developing 
fairly small gravity anomalies 
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At variance with models for the early stages of much rifting , 
which include crustal thinning by thermal and tensional influences, the 
central midcontinent here was apparently not thinned significantly. This 
• is because the crust is now anomalously thick, being 45-50 km thick under 

much of Lake Superior and at least as thick as the midcontinent average 
(about 40 km) toward the southwest. Interestingly enough, the region has 
apparently been in isostatic equilibrium since late Precambrian time, 
despite having a large gravity anomaly. 

The northern portion of the MGA has intrusives accompanying the 
basaltic lavas. These include diabase sills (e.g. Logan intrusives), 
and plutons and large cumulate bodies such as the Duluth Complex. This 
latter set of intrusives, on the Minnesota northwest edge of Lake Superior, 
is a crescent-shaped system over 200 km long, of layered gabbro. Its ba- 
sal zone has rich ore mineralization associated with it. 

As one proceeds southwards towards the central portion of the MGA, 
the origin and development of the rift system is associated more with 
block (graben) normal faulting during crustal extension and emplacement 
of the basalts as an infilling. Later, there was extensive reverse fault- 
ing on steep bounding faults, as the central horst developed (Craddock, 
1972; Anderson and Black, 1981 and 1983) . In southern Minnesota, the 
central uplifted block is called the St. Croix horst. The relative ver- 
tical displacement of the horst with respect to the deep flanking elastic- 
filled basins is up to 10,000 meters, and the sediments are thus up to 
at least this thick in places. A schematic of the hypothetical develop- 
ment of the MGA structure in southwest Iowa, based on geologic analogy 
and magnetic and gravity modelling (see a later section here) is shown 
in Figure 29 . This was formed largely 1100-1200 x 10^ yrs ago (Keweenawan) 
^ time) , but the central horst and its companion flanking basins continued 

having relative vertical movements into the Paleozoic era, based on map- 
_ ping of the post-Frecambrian sediments overlying. 

There are still some differences of opinion about whether the basal- 
tic and mafic core seam extends from the horst down to the lower crust as 
a coherent slab, or whether the horst basalt is underlain at intermediate 

■ f 

j depths by a crustal block of granite (as in Figure 29). Geophysical mod- 

elling of the MGA seems to support the latter interpretation. 
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Figure 29 . Schematic of development of MGA in southwest Iowa, 
1100-1200 million years ago (Anderson, 1981). Granitic rock has 
x T s, basalts have short dashes, sediments have dots. 



Farther to the southwest, into Kansas, the CNARS feature is pro- 
gressively deeper (to its top), and narrower in width, and more modest 
in structural extent. The sediment troughs are only about 3 km deep. 
Further, it is intrusive in nature and composition (being gabbroic, from 
drilling and geophysical modelling) rather than having much mafic (basal- 
tic) extrusives. As a consequence of these characteristics, the MGA gra- 
ity and magnetic signature is less pronounced there. To the east of the 
MGA ex-tension into Kansas is the Nemaha uplift, a subparallel but shal- 
lower and later feature. It is an uplift of Precambrian-age (1400 million 
years) granite, bounded on the east side by the Humboldt fault. This 
fault had its activation and displacement in Paleozoic time. 

The emplacement of the CNARS is a tectonic occurrence which could 
have consequences for an area which has considerable potential for econo- 
mic mineral resources. The geologic setting is favorable. Precambrian 
basement rocks representing abundant time (1100 to 2600 million years ago) 
underly a sedimentary section of younger rocks. There are major geo- 
logic structures — basins, uplifts, faults — which have had basement control 
and which have affected subsequent depositional environments. For example, 
there are large deposits of coal in Iowa, Illinois, and other midcontinent 
states, In the overlying sedimentary section. In the past year or two 
there has been initiated interest, and exploration efforts, directed at 
the oil and gas potential of the sediments along (on and adjacent to) the 
CNARS structure in Minnesota and Iowa, and in the Forest City Basin east 
of the MGA/Nemaha features. Analysis of the larger tectonic picture — 
thermal and structural origin, reactivation (neotectonism) , and possible 
effects on structures and stratigraphic conditions in the overlying sedi- 
mentary section — could help provide insight into the trapping environments 
for hydrocarbons. 

Further, where the basement rock is exposed at the ground surface 
or is shallow, as to the north (where, however, it is still often obscured 
by glacial deposits from the Pleistocene ice age) , there are rich and im- 
portant mineral deposits of iron, lead, zinc, and copper. Since the geo- 
logy of the basement rock to the south is analogous, in terms of structure 
and general evolution and the fabric of intrusive activity, the basement 
there has the potential to be the habitat of major orebodies, mineral pro- 
vinces, and similar mechanisms for economic mineralization. 




Much of the central midcontinent has a basement overlain by 
several hundred meters of sedimentary rock. It is geophysical methods 
of prospecting — including magnetic and gravity mapping — that will provide 
economical, efficient, and revealing means for reconnaissance interpreta- 
tion of the buried basement's features and bodies, composition Crock type) 
potential economic resources, structural trends and possible control of 
subsequent deposition and basin development, and of deeper crustal char- 
acter and conditions. A regional assessment (including long-wavelength 
satellite surveying) of deep structures, and their geodynamical signifi- 
cance, geologic evolution and development of an area, and consequences 
for emplacement of mineral resources, can form the basis for further ex- 
ploration to help fill the resource needs of the future. 

Although the buried basement rock of the central midcontinent is 
as yet largely unsampled by direct drilling, areal geophysical studies 
are setting the regional framework of geologic provinces and age relations 
This interpretation will be assisted by integration with a large-scale 
data set such as satellite magnetic anomaly maps. 

Figure 30 shows the geologic age provinces of the midcontinent, 
as determined by radiometric dating. 


Correlative Geophysical and Geologic information 


The regional (long-wavelength) satellite magnetic data can be com- 
bined with more detailed geologic and geophysical information to help 
interpret the geologic setting of a region and its deeper and broader fea 
tures. Such comparison can add new information on the structural and tec 
tonic framework of the area, and the relationship of geological provinces 
Figure 31 shows the thickness of the crust in the central midcon- 
tinent (Allenby and Schnetzler, 1983). The depth contours are in "kilo- 
meters below sea— level"; the total crustal thickness would be that value 
plus the elevation of the ground surface above sea-lev.'!.. In the central 
U.S., the latter is a fairly uniform couple of hundred meters. The crus- 
tal thickneje here is deduced from seismic refraction profile data. Such 
seismic coverage is quite sparse In the northern midcontinent, and where 
particularly doubtful the depth contours are dashed. 
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It can be observed (as noted by Allenby and Schnetzler, 1983) 

that, 

i) the average crustal thickness in the central midcontinent is 
about 40 km. 

ii) the crust is relatively thicker to the southwest, Oklahoma /north 
Texas (about 50 km)} and to the east, in Tennessee/Kentucky (over 
50 km) 

iii) the crust is relatively thinner to the south, in the entrant of 
the Mississippi Embayment (less than 30 km) 

iv) there is no obvious expression of anomalous crustal thickness for 
the MGA/CNARS tectonic feature that runs from Kansas up to Lake 
Superior; it has an apparent crustal thickness of 40-50 km, and 
lies along a trend that has general gradual thickening towards 
the west and to the north (above Lake Superior) . However, with 
the scarcity of seismic surveys in the region, there may yet be 
some more localized ("short-wavelength") variation in crustal 
thickness to be discovered. 

A long-wavelength gravity anomaly map of the region is shown in Fi- 
gure 32 (from Von Frese et al, 1982b) . This is from a gravity map of North 
America, with data average on a l°xl° grid, and filtered for wavelengths 
greater than 8°. We can note the gravity lows over Oklahoma/n. Texas, and 
Kentucky, corresponding in general to the areas of greater crustal thick- 
ness (and hence relatively more of the less-dense crustal material) in 
Figure 31. There is also a gravity low trending from Oklahoma to the north- 
east, toward Lake Michigan. There is a modest gravity high over the Missis- 
sippi Embayment, where the crust is relatively thinner (and where there is 
relatively more of the denser mantle rock) . 

An indication of lithospheric properties can be given by study of 
the propagation of seismic waves. Figure 33 shows the distribution of iso- 
seismal lines for equal observed (Mercalli scale) effects from two earth- 
quakes located in the central midcontinent (data from Docekal, 1970) . These 
earthquakes were of sufficient magnitude to be well reported, so that the 
pattern of isoseismal contours could be reconstructed. The earthquake 
events are described on the following page. 
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Figure 33 , 


Isoseismal contours for earthquakes in SE 
Nebraska (1935, Mercalli intensity VI) and 
NE Oklahoma (1952, Mercalli intensity VII) 
(data from Docekal, 1970) 




1 

■ The earthquakes are, 

a) 1935; epicenter near Table Rock/Tecumseh, southeast Nebraska; 
maximum Mercalli intensity I = VI. The solid dots are the re- 
cording (observing) sites, and the solid contours are isoseis- 
mal lines of equal Mercalli intensity. 

b) 1952; epicenter near El Reno, northeast Oklahoma; maximum Mer- 
calli intensity 1^ = VII (and Richter magnitude about 5.5); focal 
depth about 4-12 km; was felt in seven states in the vicinity. 

The crosses are the observing sites, and the dashed contours are 
the isoseismals. 

It can be noted that there is strong azimuthal dependence of the energy pro 
pagation. The isoseismal patterns are elongated in the northeast direction 
with similar effects being felt twice as far to the northeast as in the 
orthogonal (NW, or SE) direction, for example, for the 1952 event there 
were Mercalli intensity XV effects along the axis of the MGA up to 800 km 
away (to the northeast, in Iowa) , but extending only 325 km away in either 
the northwest or southeast directions from the epicenter. The preferred 
direction of energy propagation (with less attenuation) is parallel to the 
major NE-SW structures, the MGA and Nemaha features, and to the general 
trend of Moho (crustal base) "topography" (see Figure 31) . Seismic wave 
travel depends on rock type (lithology), physical properties, and regional 
structure. The tendency to travel with reduced attenuation in the north- 
easterly direction suggests an influence of 

i) reduced absorption of seismic energy, because of denser and 
more "rigid" mafic rock in the crust and perhaps upper mantle. 
This could be associated with a geologic province trending to 
the northeast, or the MGA structure with its major intrusive seam 
of denser (basaltic) rock of great depth extent, 
or ii) channelling of energy to the northeast, because of intracrustal 
structure and the Moho-boundary trend. 

For prospective analysis and interpretation of any component of 
remanent magnetization in the crust, contributing to satellite magnetic 
anomaly fields, we would benefit from paleogeographic reconstruction. That 
is, rearrangement of crustal blocks (geologic provinces, continental frag- 
ments) into their configuration at the time the remanence was acquired by 
the rocks. This will typically involve both rotation of landmasses, and 
translation to the appropriate geographic latitude. This can be done by 
paleomagnetism, using the magnetic inclination and declination (azimuth) 




of the remanent magnetization in rocks. This magnetization was acquired 
in the direction of the earth r s field at the time the igneous rock (e.g. 
basalt, andesite, granite, gabbro) was cooled, or the sedimentary rock 
(e.g. sandstone) was deposited. The inclination gives ancient magnetic 
latitude, and the declination gives the amount of rotation of the landmass. 
The paleomagnetic information can be buttressed by other paleolatitude and 
paleoclimatic indicators . 

The MGA complex of mafic intrusions (e.g. gabbro) and layered mafic 
(basaltic) volcanic rocks, is the largest upper- to mid-crustal geologic 
body which could be expected to have major remanent magnetization in the 
U.S. midcontinent . This is because the central seam of igneous rocks 

i) is very large: length of over 1000 km, from Kansas to lake Sup- 
erior; width up to 80 km; large depth extent, from the surface 
(in Minnesota) or a top at a few hundred meters (in Iowa and 
Kansas) down to as deep as at least 7 km (from modelling) and 
some authors would suggest much more than that 
ii) has the potential to be much more anomalously magnetic than the 
surrounding basement rock, being mafic (ferromagnesian, including 
higher in magnetite) compared to granite of the surrounding ter- 
rain. 

For interpretation of the potential magnetic signature (see later 
Section) , Figure 34 shows the paleogeographic location of the North Amer- 
ican (Precambrian) craton in Keweenawan time, 1100 million years ago. This 
is based (in part) on the paleomagnetic data from the Duluth Complex north- 
west of Lake Superior. It has a magnetic declination of remanence of 290°, 
and inclination of 42°. Thus at the time, the central midcontinent was 
oriented about 70° east of north (with respect to present geographic orien- 
tation) and at 15-25° north latitude. 


j 

i 

f 

i 


65 


DULUTH GABBRO IIOOma 


Declination 290° Inclination 42 4 


from Jahren , 64^CracJdccK'72 




msm 




■■Hi 





axinrnlha] thud-area projection 


ANDERSON 1981 


Figure 34 . Paleogeographic reconstruction of North American 
craton in late Precambrian (Keweenawan) time, 1100 million years 
ago, during formation of CNARS structure. From paleomagnetic 
data. 
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Inversion of magnetic anomalies to crustal sources 


In order to interpret crustal structure and composition, it is 
helpful to be able to invert the magnetic anomaly map to the causative 
magnetic sources. There are some differences to be considered in compar- 
ing the inversion process for satellite versus the conventional aeromag- 
netic data. 


Satellite anomaly map 

- data acquired at altitudes of se- 
veral hundred km 

- need to transform data to a common 
altitude (datum) level 

- interpretable objective: broad- 
scale (long-wavelength) sources and 
features, typically deeper 

- to invert for large areas: best to 
use spherical-earth modelling, and 
use an Earth's field that varies 
in inclination and magnitude with 
latitude 


Aeromagnetic map 

- data acquired at heights of 
a few hundred to a few thou- 
sand meters 


- objective: localized (short 
wavelength) sources and fea- 
tures, typically shallower 

- to invert for local region: can 
assume flat earth, and an Earth's 
field constant in inclination 
and magnitude over the area 


As a first approximation for modelling of crustal sources by in- 
version, it can be assumed that the magnetization is all induced. That is, 
in the direction of the local Earth's field. One can then compute, by some 
statistical best-fitting procedure, the "equivalent sources" in the erupt 
that will give the observed magnetic anomaly map. The sources could be a 
regular grid of dipoles, each oriented in the local direction of the geo- 
magnetic field. The model can be constrained by either 

i) assuming a constant thickness of magnetized layer (e.g. the crust), 
and calculating the different dipole moments (strengths) required. The ef- 

_ JL 

fective average rock magnetization, in emu/cni , can be calculated by di- 
viding each dipole moment by its grid element volume. A typical presumed 
thickness of the "magnetic crust" could be 30-40 km, which could be con- 
sidered as the effective depth from the surface to the Moho (crust/mantle 
boundary) or to the Curie-temperature isotherm. At the latter, magnetiza- 
tion is lost due to thermal effects. 


3 3 

1 emu/cm (cgs units) = 10 amp/m (SI units) 
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or 


ii) assuming a constant magnetization intensity for the crust, so 
that different calculated dipole strengths relate to different thicknesses 
of magnetized layer 
or 

iii) some combination of the above. 

The "crustal magnetization" as deduced, would be the net bulk mag- 
netization. This could be the (vector) sum of the following, 

i) induced magnetization, because rock of magnetic susceptibility 
"k" is sitting in the Earth's magnetic field H^, i.e. 

J i =kE E 

ii) viscous remanent magnetization (VRM) , acquired slowly in the 
direction of the present Earth's field, because of thermal activation. 

This is particularly the case at temperatures above, say, 300°G, in the 
mid- to lower-crust; i.e. 

°VRM 

iii) paleomagnetic remanence, such as a thermo-remanent magnetization 
(TRM) or chemical remanent magnetization (CRM) in basement rock bodies. 

This would be in the direction of the Earth's field at the time the rock 
as formed, and where the rock mass was then located geographically. Be- 
cause of plate tectonic rotations and lateral translations of continental 
blocks, the paleoremanence direction could be in a direction very different 
from that of the present Earth's inducing field locally. Call this rema- 
nence, 

J 

rem. 

Efforts to invert satellite magnetic anomaly data to crustal sources 
have suggested that, for the long-wavelength fields being modelled, paleo- 
magnetic remanence can be neglected. That is, one can reproduce the observed 
satellite fields satisfactorily well by assuming only "induced" magnetiza- 
tion (i.e. crustal dipoles in the direction of the present local Earth's 
field) (Galliher and Mayhew, 1982) . 

The net magnetization is thus 

J — J. + 
net x VRM 

= kH E + W tH E 
- k (1 + 

- k H_ 

app 
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where k *=» true average magnetic susceptibility of the crustal 
rock being modelled 

k = e ^^ ect ^ ve > or apparent, susceptibility, combining 
the effects of induced magnetization and viscous re- 
manence 


+ q vem- > 

^VRM ” Kaenigsberger ratio for the VRM, 


x.e. 


n 0 J VEM 

m vem 


where 




= magnitude of local Earth's field. 


The viscous remanence contribution could, for lower-crustal rocks 

having resided at elevated temperatures for a long period of time in the 

A 

present Earth's field , exceed the induced magnetization in magnitude. That 
is, could be greater than unity in magnitude. This makes it difficult 

to make definitive petrologic (rock type, composition) identifications of 
major crustal source bodies based on interpreted "effective" susceptibili- 
ties. It is the "true" susceptibility which would be more diagnostic, but 
this is included within the modelled "apparent" susceptibility. One is 
thus limited to differentiating rock types by general gross category of 
i) mafic character (relative proportion of ferromagnesians — iron- 
rich minerals, including magnetite, the most magnetic of the rock- 
forming minerals) , 

petrologic category typical crystalline rocks 


more mafic, 
and typically 
higher suscep- 
tibility 

v 


acidic, or sialic 
intermediate 
basic, or simatic 


granite, rhyolite 
diorite 
gabbro, basalt 


in combination with 

ii) degree of metamorphism, for continental Precambrian cratons and 
lower-crustal rock provinces. Intermediate-pressure granulite- 
grade metamorphic rocks, especially if mafic, typically have the 
strongest-intensity magnetization . 


The last major polarity reversal of the geomagnetic field was about 
700,000 years ago. 
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In the inverse magnetic modelling, it should also be recalled 
that it is the "anomalous” magnetization that is deduced. This is based 
on the contrast in properties between the source body and the surrounding 
rock. A uniformly- thick layer of constant susceptibility gives no anoma- 
lous field. In the foregoing, therefore, it is the lateral susceptibility 

contrast, ^k<or £»k , which is calculated. 

app 

After the anomaly map has been inverted to give an equivalent- 
source array of dipoles, one can then reduce the data to— the— magnetic— pole 
by reorienting the "crustal" dipoles vertically (i.e. to simulate radial 
polarization) and computing the new magnetic anomaly field. 

Experience has indicated (e.g. Mayhew, 1982; Mayhew and Galliher, 

1982) that the POGO satellite data can be reliably inverted to a grid of 
dipoles down to about 300 km apart, whereas the higher resolution Magsat 
satellite data can be inverted to a set of dipoles 220-240 (i.e. 2°x2° grid) 
km apart. That is, the inversion technique has an associated resolution 
of about half the data altitude. 

A refinement of the "constant- thickness, varying magnetization layer" 
or the "constant-magnetization, varying thickness layer" inverse modelling, 
is constraining the crustal magnetization to originate from some plausibly- 
defined source layer. Based on considerations of susceptibility magnitude, 
magnetic mineralogy. Curie- temperature depth, and magnetic modelling (see 
later Section) , this is best taken as the mafic lower-crustal region or layer 
defined by higher seismic velocities. This layer is generally found, and 
conventionally taken as, lying between the Conrad and Moho discontinuities. 

A map of the crustal thickness of the continental U.S., based on 
seismic refraction data, is given in Figure 35 (Allenby and Schnetzler, 1983) 
An excerpt of this, for the central U.S., was given here in Figure 31. The 
map shows contours of depth- to-Moho below sea level, in kilometers. The 
seismic coverage is sparse in the midcontinent region. The Moho is deeper 
(i.e. the crust has a root, more than 50 km below sea level) under the main 
belt of the Rocky Mountain (continuing down to Oklahoma/n. Texas), the Sierra 
Nevada range to the west, the Appalachian mountains, and in central Ken- 
tucky/Tennessee. The Moho is shallower (less than 30 km deep) under the 
Mississippi Embayment/aulacogen, the Rio Grande rift in New Mexico, and 
the intermontane region in the northwest between the Cascade and Rocky 
mountains . 






The thickness of the lower-crustal layer, again as identified 
from seismic refraction data, is typically about 18 km for the continental 
U.S. It is about 20 km in the midcontinent (see Figure 36) . The basal 
layer is anomalously thick (over 30 km) under Kentucky /Tennessee, and 
under western Lake Superior. Combining the thickness of the magnetic 
lower-crustal layer, with a map of dipole moment strength (see Figure 37 , 
from equivalent-source inversion of satellite anomaly data) , one can es- 
timate the effective magnetization intensity of the lower-crustal rocks 
(Schnetzler and Allenby, 1983) . 

The variation in intensity of magnetization, as deduced, is shown 
in Figure 38 . This is a map of magnetization calculated from the modelled 
dipole moment strength divided by the grid volume (surface grid area times 
the seismologically-determined basal "magnetic layer" thickness) . The ave- 
rage net magnetization of the lower crust for the continental U.S. is deter- 
mined to be 3.5±1.1 amp/m (i.e. 3.5+1. 1 x 10”^ emu/cm ) by Schnetzler and 
Allenby (1983). 

Although the seismic data coverage is not sufficiently complete in 
the midcontinent region, there are indications for the lower crust there 

■ \ 

of, 

-- more highly magnetized rock under central Kentucky /Tennessee , where 
both the total crust (Figure 35) and lower-crust (Figure 36) are also 
thicker, and where there is a large magnetic high on the satellite 
anomaly map (Figure 21) . 

— average magnetized rock under Oklahoma /north Texas, where the total 
crust is thicker in general (but not anomalously-so locally) , and 
where there is a moderately large magnetic anomaly high 

— more highly magnetized rock under central Missouri, where the crust 
is of average or less-than-average thickness, and where there is 
a satellite anomaly high that is part of a trend high to the north- 
east 

— less magnetized rock under the southern Mississippi Embayment, 
where the crust is thinner and the satellite anomaly map shows a 
large magnetic low 

— progressively less-magnetized rock as one proceeds from the central 
midcontinent (Nebraska/Iowa) towards the north (the Dakotas, and 
Minnesota) . This is accompanied by a relatively thicker total 
crust, and a satellite magnetic low. 
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Figure 36. Thickness of the lower-crustal (mafic) layer 
for continental U.S., in km. (Allenby and 
Schnetzler, 1983) 
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Figure 37 . Equivalent-source dipole moment strength, for 
continental U.S. Derived from inversion of POGO satellite mag 
netic anomaly field. Dipoles are on a grid 150 x 150 km; con- 
tour values in 10^ amp-m 2 (from Allenby and Schnetzler, 1983, 
as modified from Mayhew, 1982) . 
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Based on the assumptions, which seem reasonable in light of cur- 
rent studies, that 

i) the satellite (POGO, and Magsat) magnetics present anomalies 
having wavelengths greater than 250 km 
ii) the crustal magnetization source is primarily the mafic lower 
crust 

and iii) the magnetization is all "induced" (or at least is in the direc- 
tion of the local Earth's field), 

Schnetzler and Allenby suggest that satellite magnetic anomalies owe their 
presence, on a continental average, more or less equally to the effects 
of 

— differences in lower-crustal layer thickness 

and 

— lateral variations in rock magnetization intensity. 
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Resolution of satellite magnetics 


With the characteristic magnetic-field wavelengths ’’seen" by 
Magsat at satellite elevations of 350-550 km, it will be of interest to 
estimate the resolution of the resulting satellite anomaly maps. The 
best data set would include magnetic "quiet-time* 1 data, in sufficient 
quantity to give reliable and consistent results for a given data point, 
and from the lowest altitude possible. Figure 3 shows how the resolution 
(i.e. ability to distinguish between different sources) increases for a 
theoretical model field as one decreases in altitude. 

For the Magsat anomaly-field data set, as pre-processed by NASA 
and distributed to investigators. Sailer at al (1982) estimate that one 
can resolve crustal-source anomalies down to 250 km. This was based on 
spectral analysis of magnetic signal and noise for the Indian Ocean region, 
but the analysis was of general utility. By equivalent-source inversion 
of anomaly data, Mayhew and Galliher (1982) suggest that a regular grid 
of dipoles spaced down to 220-240 km apart can be used to represent the 
observed satellite anomaly field. 

The resolution of individual actual crustal sources will depend 
on factors such as 

i) relative magnetization of the source with respect to neighboring 
anomalous sources, and of the surrounding crustal rock 

ii) size and shape of the source, including the length-to-width ratio 
of the body or feature 

iii) orientation of an elongated source with respect to the azimuth of 
the satellite track (here, roughly north-south) 

iv) orientation of an elongated source with respect to the direction 
of the local earth’s field, in both 

— inclination: anomaly shapes and relative magnitudes will vary 

as a function of magnetic latitude, unless the data are "reduced- 
to-the-pole" 

— declination: for example, east-west borders or edges of fea- 

tures give H and H„, anomalies predominantly; north-south borders 
or edges give expression as R anomalies* 


This will be the outer layer of the Earth that has as its lower bound the 
Curie temperature (T c ) isotherm of magnetic minerals. Below this levels 
the temperature is too high for minerals to retain magnetisation, and they 
have negligible susceptibility and no magnetic remanence. This depth will 
depend on the local geothermal gradient. A higher gradient, as beneath 
regions of volcanic and geothermal activity such as hot spots (mantle 
plumes) or active rifts, will result in a thinner crustal magnetic layer. 
The depth will also depend on the Curie temperature of the dominant mag- 
netic mineral species. This is titanomagnetite, which is almost ubiqui- 
tous in crystalline crustal rocks and has a relatively strong magnetiza- 
tion. Pure magnetite, Fe^O^, has a Curie temperature of 575°C., and as 
the amount of impurity titanium increases this temperature is lowered. 
From petrologic considerations and analysis of typical (lower) crustal 
rocks, a representative Curie temperature of 555-575°C. is often ac- 
cepted (e.g. Schlinger et al, 1983) . More will be said on the variation 
of Curie temperature with magnetic mineralogy, in a later Section. 

Theoretical work, considering typical geothermal gradients and 
the magnetic mineral series titanomagnetite, 

ulvospinel - magnetite 
xFe^TiO^ . ( 1-x) Fe 3 0^ 

suggests a depth to the T -isotherm of the following (Carmichael, 1977): 

c 


for mineral titanomagnetite: 


T,C°C) 


T c depth (km) 


Fe«0 , (i.e. x = 0) 
J * x = 0.1 
x = 0.3 


m 





Depths for the effective "bottom" of magnetic sources, calculated 
and interpreted from magnetic .modelling of localized anomalies on con- 
tinents, are summarized by Green (1976) as: 

magnetic "bottom" (km) 


5 studies. North America 
3 studies, USSR 
1 study, Britain 


20 - 48 
15 - 35 
15 - 26 


Depending on the geothermal gradient, the could be above the 
Moho boundary between crust and mantle, or even above the "Conrad dis- 
continuity" locally. The latter is a seismic/petrologic boundary or 
gradational zone in the mid- to lower-crust, but not everywhere iden- 
tifiable. It divides the crust into an upper ("sialic" ) layer and a 
more mafic lower ("simatic" ) layer. The Curie-temperature could be 
shallower than this, for example, in the Yellowstone (NW 'Wyoming) re- 
gion, or the Basin and Range in southern Nevada. By contrast, a typical 
of 560°C. would be below the Moho in the"eastern U . S . "generally (geo- 
logically, east of the Rocky Mountains) . 

It has been suggested that the Moho physicochemical boundary 
is also a lower bound of the magnetic layer (s.g. Wasilewski et al, 1979). 
That is, where the Curie-temperature isotherm depth is below the Moho, 
as for some regions of oceanic crust and thin continental crust, the mag- 
netic "bottom" is up at the shallower Moho. The evidence for this de- 
rives from analysis of rocks formerly resident in the lower crust and 
upper mantle. The inference is that the upper mantle is characterized 
by non-magnetic, complex-composition spinels of "normal" structure, rather 

A 

than by the ferrimagnetic "inverse" spinel titanomagnetite of the crust. 

ii) the relative partitioning of the crust into a sialic upper- 
crustal portion , or "layer", and a more mafic lower-crustal layer . As 
mentioned, these layers or zones can be differentiated by their seismic 
expression (e.g. velocity, from seismic refraction), or petrologic comp- 
osition if the samples are available. The latter can occur if the lower 
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characterized by rocks of sialic (Si + Al) composition (lower density; 
less mafic; less magnetic) such as granite, as contrasted to simatic 
(Si + Mg) composition (higher density; more mafic — i.e. more magnesium 
and iron; more magnetic) such as gabbro or basalt. 

A contrary view, suggesting that non-oxide phases may contribute to mag- 
netization of the lower crust and upper mantle, is expressed by Haggerty 
(1979) . y9 


crust is lifted, thrust, or extruded upward by major tectonic activity. 

It is being recognized that, for modelling and interpretation 
of long-wavelength (e.g. satellite) magnetic anomalies, the distinction 
between these two crust zones— and in particular, the contribution of 
the mafic lower-crustal layer — may be very important and often dominat- 
ing. The large scale magnetic anomalies may be more dependent on this 
lower crustal layer than on the thickness of the crust as a whole. In 
one of the earlier such studies. Hall (1974) analysed intermediate-wave- 
length aeromagnetic anomalies over south-central Canada, to interpret 
the effect of major crustal structure. He excluded wavelengths less than 
60 km, attributable to localized structures and bodies in the upper crust 
and greater than 4000 km, presumably arising from geomagnetic (core) 
fields. The intermediate wavelengths correlated with the thickness of 
the lower crustal layer (having a magnetization of 5.3 amp/m) , rather 
than to lateral (compositional) inhomogeneity in the crust or to the 
crustal thickness as a whole. 

Typical properties of the crustal layers are given below (and 
see other data assemblies in Wasilewski and Mayhew, 198?.; and Schnetzler 


approx. Magnetic Properties of major crustal zones 

Magnetization J is "net” (total) magnetization 
1 emu/cm^ (cgs) = 10^ amp/m (SI) 

Upper-crust Lower-crust 

J (amp/m) J (amp/m) 


continental U.S. (Schnetzler and 
Allenby, 1983) 

continental crust, in general 
(various authors: Hinze, Hall, 
Shuey et al, Elming & Torne, 

Hahn) 

(Hinze) 

(Schnetzler) 

(Krutikhovskaya & Pashkevich,1977) 

localized intrusive (mafic) bodies 
(various authors: Mayhew et al, 
Hinze, Coles, Riddihough) 

oceanic crust (Hinze) 


3.5 ± 1.1 


2-5 


0.2 

0.1 - 0.5 

0.1 - 0.6 typically 4, 

locally up to 6 

3-5 


5 
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and Allenby, 1983 for a map of the continental U.S. shown here as Fi- 
gure 38) . The values in the foregoing table are derived and interpreted 
from magnetic modelling} e.g. by inversion of long-wavelength aeromag- 
netic or satellite anomalies, and measured in lab study of rock proper- 
ties. 

In general, then, it is interpreted that the lower-crustal layer 
has an effective magnetization, and thus an apparent magnetic suscepti- 
bility, an order of magnitude larger than for the upper crust. The 
mafic (ultrabasic) lower crust is more magnetic, and is also more mag- 
netizable in the Earth's field because it probably has relatively lower 
magnetic coercivity. The latter would be due to coarser grain sizes, and 
the elevated temperatures. 

iii) contribution of different components of magnetization . These 
will include (and see p. 68-69 here), 

- induced magnetization, due to rock with magnetic susceptibility 
being in the ambient Earth’s magnetic field 

- viscous remanence, acquired with time and particularly at elevated 
temperature. This would be in the direction of the present Earth's 
field, and could be expected to be most prevalent for the lower 
crust (deeply buried, higher temperature, more magnetizable mafic 
rock). The higher temperatures, say 300-400°C, would result in 
thermally-activated acquisition of the viscous remanence. The 
elevated temperature could also result in enhancement of magnetic 
susceptibility (the "Hapkinson effect") near the Curie temperature, 
although the relative importance of this in the lower crust is not 
yet known. 

- other remanence, e.g. thermal (acquired at the time the crystalline 
rock cooled below its Curie temperature) or chemical (acquired with 
later physicochemical changes, as during metamorphism) . This could 
be in any direction, reflecting the ancient direction of the paleo- 
field. This would depend on the orientation (azimuth) and paleo- 
latitude of the local crustal segment at the time the remanence 
was acquired. It appears that this remanence is not significant 

in the lower crust, presumably not being stable against the higher 
temperatures and for the long residence time there. 
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iv) relative metamorphic grade of the crustal (basement) rocks. 
Metamorphism is the change of rock type, and their constituent minerals, 
under the influence of higher temperature, pressure, fluid solutions, 
and long time. This can occur for localized bodies in the crust, or re- 
gionally over large areas as for lower-crustal rocks which reside in 
an elevated- temperature, higher-pressure environment. The physical and 
magnetic properties of the rocks change along with the physicochemical 
metamorphic changes to the minerals. Long-wavelength magnetic anomalies 
can thus depend in part on the metamorphic character (grade) of 

- geologic terraces of polymetamorphic character, due to experiencing 
different geologic histories and changes, as in the large ancient 
cratons (continental shields) . It is possible that over ancient 
cratons, such as the Precambrian (Canadian) Shield in north-central 
North America, the magnetic anomalies associated with the upper 
crust may be determined more by metamorphic grade and structure than 
by the inherent rock type (lithology) . This would apply for geolo- 
gic provinces of some lateral and depth extent, and which contri- 
bute to the "long-wavelength” anomalies as mapped aeromagnetically 
(or correspondingly, the "shorter-wavelength" ones of those napped 
by satellite) • 

- mid- to lower-crustal rocks, as in the basal more mafic layer. 

Since analysis and interpretation of long-wavelength satellite 
anomalies generally suggests that much of the causative magnetization 
originates in the lower crust, the degree of regional metamorphism there 
can be important. Schlinger et al (1983) suggest that different meta- 
morphic grades are characterized by: 

Metamorphic grade Typical rock type Magnetization 

low-pressure greenschist less magnetic, and less stable 

magnetically 

intermediate-pressure granulites most magnetic (especially if 

mafic) , most stable magne- 
tically 

high-pressure eclogites less magnetic, and less stable 

magnetically 
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One of the reasons that granulite-grade metamorphism in the 
lower crust results in higher magnetization may be that there is exso- 
lution of the titanomagnetite/ilmenite to produce fairly pure magnetite 
as one component. 

v) lateral inhomogeneity and differences in crustal character. 
Localized intrusions of rock having anomalous magnetization (for example, 
a mafic pluton or dike swarm) produce magnetic anomalies. Larger blocks 
of crust also produce anomalies, and these are more representative of 
the wavelength distribution "seen” at satellite elevations. Such blocks 
could be brought in conjunction by the ongoing process of continental 
rifting, translation and rotation of continental (plate) fragments , 
and then collision and suturing together, that is associated with plate 
tectonics. The lateral conditions, geologically and thus magnetically, 
that result in anomalies include variations in: 

— crustal "thickness" (of the total crust, or crustal zone down 
to the T c isotherm, or lower-crustal mafic layer) 

— geologic processes, with rocks of different petrology (composi- 
tion) , age, structural texture (e.g. intrusions, faults), mag- 
netic properties (e.g. susceptibility), or depth to T c -isotherra 
(because of the crustal blocks’ different ages and thermal his- 
tories) . 

In general, one expects regional-scale magnetic "highs" (allow- 
ing for variations of anomaly appearance with magnetic latitude) to be 
associated with one or more of: 

— thicker crust 

— thicker lower-crustal layer (since is more mafic, and resides 
at moderately elevated temperatures) 

— deeper Curie-temperature isotherm (and lower geothermal gradient) 

— higher magnetite content of rocks 

— intermediate-grade (granulite) metamorphism 

— borders or boundaries between crustal terranes of different 
thickness, rock type (composition, age), structural style. 
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The characteristics of a magnetic anomaly map, that can be used 
for magnetic interpretation of crustal character, are: 

-- magnetic highs and lows (or positives and negatives, if the 
zero-level is arbitrarily set at a median level) 

— magnetic gradients; these can indicate geologic and structural 
trends, and boundaries. Gradients can be accentuated by second- 
derivative analysis. 

— areas of different magnetic character and pattern. Different 
textures and "morphologies" of magnetic pattern could indicate 
different tectonic blocks (now sutured together) , different 
lithologies, or different regional metamorphic grades. 

Satellite magnetic anomaly maps provide broad-scale (continental- 
wide) coverage with observed wavelengths characteristic of geologic pro- 
vinces and features of crustal dimensions. That is, one can hope to 
"see" anomalous crustal sources of regional lateral extent, and of great 
depth extent. The latter would extend down to the Curie-temperature iso- 
therm, or Moho boundary, whichever is shallower. These crustal sources 
would include continental fragments (or blocks, or microplates). Thus if 
the good-resolution Magsat data can resolve and distinguish these blocks, 
and their extensional or collisional boundaries, then satellite surveying 
can help interpret the patterns of plate tectonics. This would apply for 
both: 

— modem plate tectonics. Magnetic anomaly mapping might help 
interpret the thermal regime (geothermal gradient, variation in depth to 
T c -isotherm) and any associated crustal or lithospheric thinning, at 
active rifts at extensional plate boundaries, or at localized hot spots. 
Sample study areas would include the midocean spreading rifts, and Ice- 
land, Hawaii, Yellowstone, and the Red Sea. 

and 

— ancient plate tectonics. This would involve old paleorifts 
where continental rifting failed (i,e. aulacogens) , or collisional zones 
where continental blocks have been welded together. 

The expected relationship between magnetic anomalies and the 
depth to T c ~isotherm in volcanic/geothermal/rif ting areas is that the 
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elevated Curie temperature isotherm will give a thinner "magnetic crust" 
and thus a more subdued magnetic anomaly. That is, a magnetic "low" 
compared to the surrounding crust of normal, or average, thermal charac- 
ter. This relation may not apply, however, for: 

— very old thermal tectonism. After the original thermal pertur- 
bation and associated effects, the crust would have equilibrated 
thermally and the T c -isotherm returned to its normal lower level. 
There would thus now be no anomalous thermal signature to con- 
tribute to a magnetic anomaly. Any anomaly would be due to the 
residual differences in crustal composition, or magnetic prop- 
erties, or thickness of crustal layers. 


or 

— very young thermal tectonism. With recent volcanic or plutonic 
activity, the intrusion upward of molten or semi-molten rock 
and geothermal fluids would be accompanied by an elevated T^- 
isotherm. This would be a local phenomenon, however, since 
rock is a poor thermal conductor. It would take some time for 
the T c isotherm to become established and equilibrated shallow 
over a broad-enough region and result in changes in the long- 
wavelength anomalies seen by satellite. 


(b) Related interpretations 

Now let us consider general interpretations that can be made re- 
lating satellite magnetic anomalies to causative crustal sources. A 
number of major crustal geologic features are expressed on NASA’s global 
Magsat anomaly map, as interpreted by Frey (1982a and b) . The spatial 
relationship of anomalies to crustal sources is even better defined if 
one uses a reduced- to- the-pole version of the magnetic map, as done for 
POGO data by Langel and Frey (1983). In general, with respect to the 

A 

observed magnetic anomalies.- , 

i) ancient (i.e. Precambrian-age) continental shields are generally 
magnetic highs (positive). Some are not, however (e.g. southern Superior 
province in North America, and southern India) , and there can be varia- 
tions in the magnetic texture and signature due to differences in regional 
metamorphic grade of the crustal rocks. In general, a greater degree of 
metamorphism produces a higher (more positive) anomaly expression. 

* The largest magnetic anomalies measured by Magsat, after "normalizing" 
for local geomagnetic field intensity (i.e. induced-magnetization ano- 
malies are greater at higher latitudes, because the earth’s field in- 
tensity is stronger there) are, in order, Bangui (west Africa), Kursk 
(USSR), Kentucky /Tennessee, and Tibet. 
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ii) ancient structures seem to have anomalies arising mainly from 
intracrustal lithologic variation. An example of this might be the large 
Bangui anomaly in the Central African Republic. Here (Regan and Marsh, 
1982; Hastings, 1982) , magnetic modelling suggests a central Precambrian 
shield (with depth of about 40 km) flanked by younger and less-magnetic 
sedimentary basins. To the north is an arm of the Chad basin, with depth 
about 13 km, and to the south is the Congo basin, with depth about 10 km. 
A possible alternative model is a paleorift. However, modelling of the 
observed anomaly apparently requires effects of more than a single body 
having a near-horizontal magnetization (with such inclination, because 
the feature is near the magnetic equator) . 
as compared with 

younger structures and tectonic features (typically, of Ceno- 
zoic age — i.e. younger than 70 million years), which seem to have ano- 
malies resulting from variations in Curie-temperature-isotherm depths. 

iii) rifts should have a magnetic expression, depending on the age 
and maturity of the tectonic development. Although the process could 
be variable, one could perhaps make a tectonic distinction between 

— "aulacogen-style" rifting, to yield failed paleorifts. This would 
have doming and arching of the lithospheric plate, due to the 
thermal upwelling in the upper mantle, but leaving the lithosphere 
largely unfractured. 

and 

— true rifting, with extension. This would have thinning and frac- 
turing of the lithosphere, leading to lateral movement apart of 
the plates (i.e. continental drift) . 

For an early stage of rifting of continental crust accompanied by crus- 
tal thinning and elevation of the T^-isotherm, both effects would pro- 
duce a magnetic anomaly low. This characteristic should endure with time 
except that the thermal regime would equilibrate back to normal and the 
T^-isotherm contribution would diminish. There are (on reduced- to-the- 
pole anomaly maps) magnetic lows over the Rio Grande rift (New Mexico; 
an active "rift" feature but with little associated volcanism) , East 
African rift, and over such (failed) paleorifts as the Mississippi Em- 
bayment/aulacogen, the Benue trough (an aulacogen rift, north of the Ban- 
gui anomaly, presumably from a triple junction associated with the 
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opening of the present Atlantic) , Amazon River valley, and several in 
the Siberian continental platform (Frey, 1982a and b) . With tectonic 
evolution of a rift, there may be injection of mafic rock in the axial 
crustal zone, and emplacement of a rift "pillow" below the rift f s cen- 
tral thinned crust. The latter rock, a mafic differentiate of the up- 
per mantle, would be ultrabasic and highly magnetic (high susceptibi- 
lity) . Such material, because it offsets the crustal thinning and con- 
tributes more-magnetic rock, results in a thicker "magnetic lithosphere" 
locally. This would yield a magnetic anomaly high, or at least partially 
negate the pre-existing magnetic low. There can thus be a combination 
of effects. For example, the magnetic low over the Mississippi embay- 
ment has been attributed to crustal thinning, elevated isotherms, and 
if it does have an axial lower-crustal or sub-crustal rift pillow (the 
crust locally is 35-40 km thick) then this is anomalously less magnetic. 
This may be (Thomas, 1983) because of the presence of higher-titanium 
titanomagnetite in the lower crustal rock, which may be gabbroic in 
composition. Titanium reduces both the magnetization and Curie-temp- 
erature of magnetite; this also raises the T^-i so therm depth without 
an increase in geothermal gradient. 

A related aspect of rifting would be the thermal (and possibly 
physical) thinning of the lithosphere as a tectonic plate moves over 
a relatively fixed thermal disturbance, such as a mantle plume (hot spot) . 
An example would be the Hawaiian islands, where the Pacific oceanic 
plate moves west-northwest over a plume. Depending on the plume’s 
areal extent, temperature, and magmatic activity, as well as the ini- 
tial thickness of the lithosphere, the rate of plate movement, and the 
thermal conductivity of the lithospheric rock, there would be a lag 
in the recovery of the thermal thinning back to normal. For a fast- 
moving plate, there could be a small magnetic anomaly low displaced 
a little "downstream" from the source plume. 

iv) most submarine platforms and rises are magnetic highs, as 
contrasted to many oceanic abyssal plains and basins which are lows. 

v) several subducting lithospheric slabs are magnetic highs 
(Langel and Frey, 1983). This would be due to the cooler slab penetra- 
ting down into the upper mantle, carrying the isotherms (including the 
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T c ~isotherm) with it . The magnetic anomaly would be expected to be 
greater for more rapid subduction , and would reduce with time as the 
subducted lithosphere becomes thermally equilibrated to the higher temp- 
eratures. 

vi) continental crust has larger-amplitude anomalies than oceanic 

crust. This is because the latter is thinner* and has a shallow T -iso- 

c 

therm* even though oceanic crustal rock is more magnetic (per unit volume) 
than typical continental crust. 

(c) Interpretation of study area 

The processed Magsat magnetic anomaly map (scalar H^, redueed- 
to-the-pole) for the central II. S. midcontinent was shown in Figure 21, 
and is repeated here as Figure 39 . Same of the related geologic and geo- 
physical features are given in Figure 24 (major tectonic structures) 
which is included here as a transparent overlay on Figure 39, and also 
Figure 30 (basement age provinces) , Figure 31 (crustal thickness) , and 
Figure 36 (thickness of lower-crustal layer) . 

Several major anomalous features are observed on the anomaly 
map of Figure 39. They are: 

i) a bullseye magnetic high over east-central Kentucky /Tennessee 
(centered at 37%°N. Lat., 84^°W.long.) . This is one of the most promi- 
nent anomalies on the Magsat anomaly map of the world, and was identified 
also on the earlier P0G0 satellite data. 

On the reduced— to-the— pole map here, the anomaly has a central 
maximum of about 15 nT with respect to the general regional field value. 

It is generally attributed to a mafic basement rock complex in the upper 
crust, and a lower crust intrusion which is also anomalously ultrabasic. 
This occurs at the intersection of the southeast-trending "south-central 
magnetic lineament" and the southwest-northeast-trending "New York-Alabama 
lineament" (see Fig. 39). It may also be along the southward extension 
of the Grenville front, which comes down from Lake Huron and separates 
older rock on the west from younger rock (approximately Keweenawan age, 
about 1100 million years old) on the east. 

* 

Subduction rates are typically several cm per year. By modelling of the 
thermal regime and rock properties, it is deduced that for very small 
subduction rates (well under 1 cm/yr) , temperature isotherms are "dragged 
down" in the lithosphere slab faster than heat can be conducted upward 
to maintain lateral equilibrium. This is for the case of heat trans- 
port by conduction alone. 
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The basement rocks here are mafic volcanics and metamorphics . 

The former are similar to the Keweenawan (late-Precambrian) volcanic 
rocks of the western Lake Superior region, and may also be of the same 
age (Hildenbrand et al, 1983) . The high-resolution aeromagnetic map 
for the region is a jumble of high-amplitude anomalies, in contrast to 
the satellite expression (in effect a high-pass wavelength-filtered ver- 
sion) of a large, single high. 

Geophysical modelling of the anomaly (Mayhew, Thomas, and Wasi- 
lewski, 1980; J. Phillips/U.S .G.S .) suggests a feature such as portrayed 

in Figure 40 . A mafic intrusion such as gabbro in the crust would give 

3 3 

the increased density (3.0 gm/cm , versus crustal value of 2.8 gm/cm ) 

and magnetization (anomalously greater by 4 amp/m, or 0.004 emu/cra ) to 
yield the positive gravity and aeromagnetic anomalies. The intrusive 
body is interpreted as extending from 5 to 25 km down (with the Moho at 
about 30 km) , and has a width extent of about 40 km. This modelling is 
thus based on an effectively thicker localized lower crustal "layer", 
for an otherwise constant-thickness (i.e. as gauged by depth-to-Moho) 
crust. This is consistent with seismically-determined crustal thick- 
ness in the region (see Fig. 31, Allenby and Schnetzler, 1983) of 35 km 
in east-central Kentucky/Tennessee, although the crust thickens to the 
southwest and to the east to 50 km. There is a thickening of the lower- 
crustal layer from a regional value of about 20 km to a thicker 30 km, 
in an area extending from central Kentucky through central Tennessee to 
northern Alabama (see Fig. 36) . There is in fact a corresponding elonga- 
tion of the magnetic anomaly high in this direction (Figure 39, again). 

In the Interpreted/calculated map of lower-crustal magnetization (Fig. 
38), an anomalous value of 4.5 amp/m is assigned for the localized re- 
gion of east-central Kentucky/Tennessee. This is close to the value of 
4.0 amp/m deduced for the simple modelling of the aeromagnetic data. 

ii) a magnetic high over NW Texas /Oklahoma (centered at 36°N. 
lat., 101%°W. long.), then extending in an arcuate moderate high up to 
the northeast through Missouri to Lake Michigan. While this anomaly 
could also be noted on P0G0 data, its continuation up to the northeast 
is better evidenced in this data here. 
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THROUGH MODEL 


Figure 40 . Two-dimensional magnetic and gravity modelling of 
. Kentucky /Tennessee magnetic high. Profiles are west-east; 
solid lines are measured anomalies, dashed lines are modelled. 
Densities in gm/cnr', magnetization contrast is 4 amp/m (0.004 
emu/ cm^) a (Mayhew, Thomas, and Wasilewski, 1980) 



The magnitude of the anomaly over NW Texas/Oklahoma is about 
9 nT. The origin of this anomaly is not well known. There is a gravity 
anomaly low centered over the area (Fig. 32), implying a thicker crust, 
and the crustal thickness map (Fig. 35) shows anomalously thick crust 
(50 km) in the larger region. There is insufficient seismic coverage 
to define the local crustal thickness there in detail. If one attri- 
butes the localized magnetic high to crustal thickening (i.e. a thicker 
magnetizable crust) , there is less need to invoke anomalously more-mag- 
netic crustal rock. In fact the lower-crustal magnetization map (Fig. 

38) interprets little if any such anomalously magnetic rock in the area. 

The source of the anomaly high trend towards the northeast is 
also uncertain. There are a couple of possibilities: 

1)- the origin is a zone of greater crustal thickness, perhaps deve- 
loped in late Precambrian time (Black, 1981) either during the Mazatzal 
orogeny, or where the subsequent and subparallel Granite-rhyolite terrane 
of the upper crust now is emplaced. The Mazatzal was a tectonic/deforma- 
tion episode, possibly at a convergent plate boundary, of 1600-1700 mil- 
lion years ago, and the Granite-rhyolite terrane has an age of 1400-1500 
million years. These geologic provinces of the Precambrian basement are 
shown in Figure 30, as developed from radiometric age-dating and basement 
petrology (rock types) . The Mazatzal belt ("M") runs from NW Texas/Okla- 
homa up through Missouri to Lake Michigan, and the Granite-rhyolite ter- 
rane ("G-") sweeps up on the southern and eastern margin. Additional and 
more recent borehole and age-date information suggests that the northern 

boundary of this latter terrane should be placed farther north into NE 

* 

Missouri, SE Iowa, and northern Illinois. Thus either geologic trend 
(Mazatzal, or Granite-rhyolite), or both, could contribute to the elon- 
gate magnetic anomaly high. 

The gravity anomaly map (Fig. 32) does show a low along the trend 
of the magnetic high and the Mazatzal belt, implying a thicker crust. 
However, the crustal thickness map of Fig. 35 shows no anomalous thick- 
ening in a fairway up to the northeast, but rather a general thickening 
in the northern midcontinent from east to west. Unfortunately, the seis- 
mic coverage is incomplete. On the other hand, there is also the fairly 
remarkable tendency for seismic energy in the crust to be propagated pre- 
ferentially in the northeast direction from Oklahoma, as seen in Fig. 33 
for earthquake isoseismal con tours . 

* although there is a possibility that the rock bodies sampled, and whose 
younger ages would extend the Granite-rhyolite terrane farther north, 
represent later events or intrusions in the older Mazatzal terrane. 
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2) the origin of the anomaly high trend toward the northeast is 

A 

due to a zone in the lower crust which is more magnetic. For example, 
in central Missouri the crust is of average thickness (about 40 km) , but 
the interpreted lower-crustal magnetization is higher than average (5 to 
6 amp/m; see Fig. 38). 

iii) a magnetic low over the Mississippi Embayment/aulacogen , 
over the Gulf Coast and extending from about 91—94° W. longitude. The 
anomaly has a minimum (low) of about —15 nT over the Gulf of Mexico coast 
of Louisiana, and extends north-northeast up the Mississippi River (and 
the Reelfoot rift) to the New Madrid seismically-active zone in southeast 
Missouri. 

Gravity and magnetic mod filing of the region, in addition to 
what is known of the sub sit? .ace geologic structure from seismic, drilling, 
and other information, suggests the probable origin of the regional mag- 
netic low (Von Frese et al, 1980 and 1982b; Keller et al, 1983) . The 
Embayment is a broad reentrant of Mesozoic and Cenozoic— age sedimentary 
rocks extending into the Pracambrian craton. It originated with the late 
Precambrian/early Paleozoic rifting of the proto-Atlantic at the "rifted 
continental margin" (see Figure 39) to the south. Associated with this 
were aulacogen features such as the Reelfoot rift. With the (second) 
opening of the Atlantic in Mesozoic time, the Embayment formed with areal 
subsidence. The central rift and associated structure, particularly at 
the northern end, are presently still active seismically . 

There are major and broad anomalies — a gravity high, and a (FOGO, 
and Magsat) magnetic low — over the Embayment. This negative correlation 
between long-wavelength magnetic and gravity anomalies for large crustal 
sources, and the positive correlation generally observed between sat- 
ellite magnetics and continental crustal thickness, suggests that the 
anomalies here are due to crustal thinning. This would presumably have 
resulted from the paleorift development. The thinning would give the 
gravity high, due to upwarped denser mantle, and the magnetic low due to 
a thinner magnetizable crust. 

The crustal thickness map (Fig. 35) shows that there is some 
localized thinning of the crust, to 25 km, compared to the general 

a 

An ultramafic and more magnetic lower crust here might have originated, 
or been associated, with a continental collision/ subduction which helped 
form the miacontinent craton in Precambrian time. 
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vicinity which has a thickness of about 40 km. However, the thinning 
does not appear to reflect the elongated nature and length-extent of the 
observed magnetic low. Another contributing factor could be a shallower 
T c ~isotherm depth (to produce, for example, the localized magnetic ano- 
maly low over Yellowstone) , but there is no substantial evidence for a 
higher regional geothermal gradient in the Mississippi embayment. 

Another possibility for the magnetic low, and suggested by the 
geophysical anomaly modelling, is that the local lower-crustal layer 
has the usual higher density (to contribute to the gravity high) but an 
anomalously low magnetization. To explain this, one could invoke miner- 
alogic/petrologic conditions, such as more titanium impurity in the mag- 
netite, to reduce the magnetization. The interpreted lower-crustal rock 
magnetization (Fig. 38) shows a value under Louisiana of about 2-3 amp/m, 
which is below the regional (and continent's average) value of 3.5 amp/m. 

iv) a magnetic low over South Dakota , with magnitude of about 
-6 nT. This could be due to some combination of: 

— elevated Curie- temperature Isotherm; however, regional geothermal 
gradients are not well known in the area. Preliminary recent heat 
flow data suggest that northern Nebraska (and perhaps southern S. 
Dakota) may have higher heat flow (at least in the upper crust) 
than the U.S. midcontinent in general (Goswold, 1980). 

— crustal thinning; however, the crust is about 45-50 km thick under 
South Dakota itself, which is at or a little thicker than the ave- 
rage for the northern midcontinent. On the other hand, the crust 
is even thicker, over 50 km thick, to the north and south of S. 
Dakota where the magnetic low is centered. 

— anomalous magnetization deficiency in the lower-crustal layer. 
There is a hint of this in Fig. 38, suggesting a local magnetiza- 
tion of 2. 4-2. 8 amp/m, which would indeed be below the U.S. conti- 
nental average of 3.5 amp/m. Such a magnetization deficiency 
could be due to less-mafic rock having less magnetite (e.g. being 
diorite, rather than gabbro) , or rock with a lower degree of re- 
gional metamorphism than granulite-grade, or magnetite with more 
titanium as an impurity (and thus lower magnetization and lowered 
Curie temperature) . 

* 

and the lower crust here may also have an anomalously higher density 
for this basal layer, as a relic of the aulacogen rifting 
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v) a magnetic low over west-central Georgia (south-southeast 
of the Kentucky /Tennessee anomaly) , with magnitude about -8 nT . This 
is probably attributable to a localized region of anomalously less-mag- 
netic lower crust, as interpreted in Figure 38. The value assigned 
from analysis is 1. 5-2.0 amp/m, again well below the continental ave- 
rage of 3.5 amp/m. 

vi) a general tendency for the magnetic pattern to follow and 
reflect the southern "rifted continental margin ” of the North American 
craton (see Fig. 39) in this area. This margin was formed in late Pre- 
cambrian/early Paleozoic time, with the opening of the proto-Atlanti - , 
and had associated with that rifting event a number of failed rift arms 
(aulacogens) directed into the craton. These include the Reelf oot rift, 
and the southern Oklahoma aulacogen (Keller at al, 1983). The margin 
was then the site for the Appalachian and Ouachita mountains, formed in 
the later collision closing the proto-Atlantic. 

Referring to Fig. 39, it is seen that the ancient continental 
margin approaches from the east, wraps concave-northward around the con- 
tour lines of the large Kentucky/Tennessee magnetic anomaly high (and 
also the south end of the Appalachians) , then wraps concave-southward 
around the contour lines of the magnetic low associated with the Missis- 
sippi embayment, and then wraps concave-northward around the south side 
of the large NW Texas /Oklahoma anomaly high. 

vii) a lack of apparent magnetic expression for the Central North 
American Rift system/Midcontinent Geophysical Anomaly (CNARS/MGA) . This 


long, deep, and major paleorift — an immense seam of mafic rock in the 
upper crust — would be expected to have a magnetic anomaly, even if due 
only to induced magnetization. That it does not — there is no satellite 
anomaly, and even the aeromagnetic signature is modest in amplitude — may 
be due to a combination of factors: 

— there is a combination of normal and reversed remanence, the latter 
being of opposite polarity and acquired during a time period when the 
Earth’s geomagnetic field was reversed. The adding of opposing magnetic 
effects serves to reduce the net total magnetization of the anomalous ma- 
fic (largely basaltic) seam, in terms of the external anomalous field pro- 
duced. Geophysical modelling, to be presented in the next Section, sup- 
ports this interpretation. 
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— the CNARS main arm, from Kansas over 1000 km up to western 

Lake Superior, is oriented along the north- south azimuth of the satellite 
orbit track. It is thus less likely to show an anomaly than an east- 
west trending structure or boundary. This is because a narrow, along- 
track anomaly may be degraded in the data processing, by corrections 
for geomagnetic field removal, for base level adjustment due to external- 
field transient effects Ci»e» taking out long-period transient varia- 
tions, which can "look" like crustal anomalies on a track profile) , or 
for adjustment to observed fields made for altitude differences of the 
orbit along a track or between different passes. 

— the feature is aligned with its long axis roughly along the azi- 
muth of the present Earth's field. Thus the anomaly is less well defined 
(e.g. in terms of distribution of "free poles" on the causative source 
body) than if it were oriented east-west. 
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Geophysical modelling of CNARS /Midcontinent Geophysical Anomaly 
— magnetics and gravity 

As described previously ("Geologic Setting", page 50-57), the 
major arm of the Central North American (paleo)Rift system has its 
greatest tectonic development, and largest gravity anomaly expression, 
in the central portion. This is in the state of Iowa, between the 
MGA extremities in Kansas to the southwest and northern Minnesota /Lake 
Superior to the north-northeast (see Fig. 24-26, and Fig. 39 and overlay). 
The basement rocks for this late-Precambrian aulacogen are covered by 
younger sediments, and have been penetrated by drilling only infrequently 
and not in a geologically-sy sterna tic way. The best evidence for the 
nature, extent, and properties of this huge crustal structure — a seam 
of intrusive and extrusive mafic rock in an otherwise granitoid upper 
crust— remains the interpretation of gravity and magnetic anomalies. 

Beep probing by reflection and refraction seismology is rare in the 
northern midcontinent, although detailed reflection work 

— has been done in the past couple of years by the COCOHP univer- 
sity consortium, on profiles in northern Minnesota (because of 
interest in the tectonic westward extension of the "Great Lakes 
Tectonic Zone") and northeast Kansas (for interest in the south- 
ern extension of the CNARS and the subparallel Namaha uplift) 

— is being done commercially at present (for the first time) in 
Iowa, as "frontier" reconnaissance exploration for oil and gas 
potential associated with the MGA structure and related basin 
development . 

The MGA' s gravity anomaly map is shown in Fig. 27, the lowa-region 
aeromagnetic (total-field) anomaly map in Fig. 28, and a schematic of 
its hypothesized development in the area in Fig. 29. The anomalous ex- 
pression for the MGA arm, a deep-seated failed paleorift, is over 1000 
km long and up to 150 km wide. The maximum gravity difference between 
the central uplifted basalt-dominated horst, and the flanking deep sedi- 
mentary basins, exceeds 160 milligals. The aeromagnetic expression for 
this major crustal inhomogeneity with highly anomalous magnetic prop- 
erties compared to the adjacent "country rock", is a belt of elongate 
and linear anomalies with a well-defined trend but very modest ampli- 
tudes. Anomaly variations over the MGA, in the central portion in Iowa, 


are generally less than several hundred nanoTeslas, whereas anomalies 
due to localized basement intrusions nearby can exceed a couple of 
thousand nanoTeslas. 

Researchers here and at the Iowa Geological Survey 'have been 
doing gravity and magnetic modelling of the MGA feature, in order to 
understand the dimensions of the source structure, the composition and 
physical rock properties of it and its crustal environment, and its 
evolution as a paleorift (Anderson and Black, 1981) . This work will 
also help interpret a magnetic expression seen (or not seen) at Magsat 
satellite altitude. A model for the feature's crustal evolution may 
also provide clues to the potential emplacement and distribution of 
mineral deposits or other resources. 

The above authors modelled for over a dozen profiles, across 
the MGA feature in southwest Iowa (see Figure 41) . Typical profiles 
are shown in Figures 42, 43, and 44 . In this region, the maximum gra- 
vity anomaly difference is about 140 milligals, and the variation of 
aeromagnetic anomaly (scalar) field is about 600 nT. For a given pro- 
file, first a rough density model was created to approximate the ob- 
served gravity anomaly. This was then modified to give a calculated 
magnetic anomaly matching the observed magnetic anomaly. The model was 
then adjusted to have its gravity anomaly more closely fit the observed 
anomaly, and finally refined to match the magnetic anomaly even more 
closely. Using plausible lithologies and rock properties to fit the ob- 
served gravity and aeromagnetic fields, they suggest a central uplifted 
block of mafics (primarily basalt), about 45-60 km wide, extending from 
a top (i.e. of the crystalline basement) at a few hundred meters, down 
to a bottom that in places exceeds a depth of 7,000 meters. This horst 
is floored by normal-crustal granitic rock, and bounded laterally by 
high-angle faults. It is flanked by sedimentary troughs of Precambrian 
clastic rocks, which can have depths as great as 10,000 m. The central 
basaltic seam should have physical property contrasts to the adjacent 
sediments and granite. Typical ones used for modelling are listed on 
page 103. 
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ER GRAVITY ANOMALY MAP 
OF SOUTHWEST IOWA 


Location of profiles for geophysical modelling 
of CNARS/MGA paleorift structure, in southwest 
Iowa. On Bouguer gravity anomaly map; C.I. = 

5 milligals. 
























_ £ Bouguer Gravity Anomaly (mgols) 

Depth (Km.) jje » (minus - 45 mgal regional) 


original page iS 

OF POOR QUALITY 


MGA Profile E-E' 


\\ 

^ model 

gravity \ 


observed 
/ gravity 


model / \\ 

magnetics y ^ 

^-^•‘ : ^--<l^)bserved 

magnetics 


\r\ 




:o: 


Total Field Magnetic Anomaly 





Parameters for Modelling of CNARS/MGA In Iowa 


Structure component 


Central horst 


Flanking 

Sediments 


Rock type: 


Property: 

density (gm/ cur) 

susceptibility contrast 

(emu/ cm^) 

remanence gontrast 
(emu/ cm; 


basalt-dominated elastics 


2.99 
4-4.5 x 10 


2.47 


Surrounding 

Basement 


granite 

2.67 


-4 


1-1.2 x 10 


The magnetic susceptibility alone of typical basaltic rock should 
result in a large magnetic anomaly for a body of this size. Since it 
does not, the implication is that either, 

— the mafic block is anomalously non-magnetic . This is an unlikely 
prospect for basalts here, especially since where analogous rocks 
are accessible or exposed at the surface (to the north) , they 
have typical mafic magnetic properties 

or * 

— the rock "pile 11 has both normal and reversed magnetic remanence, 

which combine to give a net reduced magnetization (induced, plus 
the normal and reversed remanence) for the crust as a whole. 

The latter possibility is very plausible. This is because there 
is a major geomagnetic field transition from (older) reversed to (younger) 
normal in Keweenawan time and dated at about 1100 million years ago, as 
measured for rocks of the same age and geologic/tectonic origin exposed 
around western Lake Superior. Since the central block of mafic rocks 
is largely layered, with the oldest rock layers on the bottom, the pre- 
sumption that it is the bottom portion of the block that has reversed- 
polarity remanence is a reasonable one. The emplacement of the MGA in- 
trusives and extrusive lavas thus apparently "caught" a magnetic polarity 

& 

"normal" polarity remanence does not imply greater normalcy, but 
rather that it represents and reflects the configuration of the geo- 
magnetic field as as present, i.e. with the "north" magnetic pole 
near the north geographic pole. 
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transition of the Earth’s field. 


The direction of remanence in the rock is introduced as appro- 
priate for Keweenawan time. This is a declination (azimuth of the paleo- 

- -JL 

field) of 290 (i.e. west-northwest presently) , and paleoinclination of 

40°. This information is obtained from lab study of CNARS rocks where 
exposed farther north, and the paleomagnetism thus says that in Keweena— 
wan time, the U.S. northern midcontinent was located geographically closer 
to the equator, and was rotated (70°) clockwise, compared to now. 

The rock fractions with normal and reversed remanence are as- 
signed the identical value for susceptibility, since there is no reason 
to think that petrologic composition varies with changes in the terres- 
trial magnetic field. The relative contributions of the normally and 
reversely magnetized zones can be adjusted to give the net observed mag- 
netic anomaly. This is done by changing the relative volumes of them 
in the central mafic horst. A further refinement to modelling that could 
be considered is that the remanent magnetization acquired when the field 
was reversed may be greater (of stronger intensity) that the remanence 
acquired when the field was of normal polarity. This is because paleo- 
intensity determinations for Keweenawan— age rocks from the western Lake 
Superior region (Pesonen and Halls, 1982) indicate that the Earth’s 
field (there, at that time) was about 40% stronger when it was of reversed 
polarity, than when it was "normal”. Because of the stronger remanence 
that would result, the relative volume fraction required of reversely- 
magnetized rocks would be reduced somewhat. The authors further note 
that, based on experimental study, in late Precambrian (Keweenawan) time 
the earth’s field in general was about 40% stronger than it is at the 
present time. 

A schematic diagram for the development of the MGA in the re- 
gion of southwest Iowa was shown in Figure 29. The geologic setting 
has also already been described. It is important to understand the evo- 
lution of continental rifts, and ancient aulacogens, so that we can 
better understand tectonic regimes on continental cratons, plate tectonic 
processes in the past and any subsequent resurgent activation (such as 
seismic activity, and uplift or subsidence which has consequences for 
basin development) , and relation to emplacement of mineral deposits or 

* i.e. N.70°W. 
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other resources. 

The detailed modelling by Anderson and Black here clarifies 
several matters: 

i) while the CHARS /Midcontinent Geophysical Anomaly in the nor- 
thern Minnesota/Lake Superior region may have a central seam of fissure- 
fed successive flood basalts deposited along a subsiding, trend (Green, 
1983) , the structure farther south (as in Iowa) was characterized during 
formation by major graben development with steep bounding faults. The 
axial rift valley progressively filled with successive volcanics and in- 
trusions, perhaps brought up along the fault zones bounding the central 
graben block. 

ii) the central block, now uplifted with reactivation to form a 
horst, is floored (at a depth of 3-7 km, in the model profiles of Figs. 
42-44) by crustal granite. This is the crustal block that subsided with 
the rift-initiated graben development. This is in contrast to speculation 
that the central mafic seam extends all the way down to the middle or 
lower crust. If the latter were the case, the mass anomaly (due to the 
density contrast of the mafics /basalt compared to the upper-crustal gran- 
itoid rock) would give a gravity anomaly over the MGA much greater even 
that it is now. 

iii) retaining reasonable petrologies (analogous to CNARS rocks where 
they are exposed farther north) and corresponding density and magnetic 
properties, one can explain the observed gravity anomaly by the structure 
proposed, and the observed modest magnetic anomaly by partitioning the 
central mafic block into zones of opposite remanence. This involves a 
lower, older zone of reversed remanent magnetization, and an upper 
younger zone of normal remanent magnetization. In support of this, there 
is a known polarity transition from reversed to normal paleofield (as 
reflected in remanent magnetizations) in rocks of Keweenawan age, 1100 
million years old, in the region. 


Resources potential 


The U.S. midcontinent has some important economic resources, 
from the major oil and gas provinces of the Texas/Louisiana/Mississippi 
Gulf coast (both onshore and offshore) in the south, and substantial 
reserves in sedimentary basins farther north in Oklahoma and Kansas, to 
the mineral deposits of copper, zinc, lead, and iron ranging northwards 
from Missouri through Wisconsin to northern Minnesota and upper Michi- 
gan. The region has considerable potential for discovery of new re- 
sources, particularly of mineral deposits associated with the buried 
basement rock. This is because much of the northern Midcontinent is 
covered by a veneer of glacial deposits from the most recent (Pleisto- 
cene) period of glaciation and, except for the craton exposed in northern 
Minnesota, Wisconsin, northern Michigan and farther north, the region's 
basement is overlain by a sequence of sedimentary rocks deposited on an- 
cient marine platforms and in subsiding basins. This "overburden" ham- 
pers exploration and development of mineral deposits in, or immediately 
overlying, the crystalline basement rock. 

The region, particularly the northern half where the basement is 
not too deeply buried, holds much promise for discovery of new resources 
for the future. This is because the basement rock is often analogous in 
age, geologic setting, tectonic evolution, and geophysical "texture" 
(aeromagnetic and gravity anomaly patterns of intrusives and structures) , 
to terrane which is exposed or accessible and which is economically min- 
eralized. For example, the Precambrian basement rock of Iowa, generally 
buried beneath a few hundred meters of sedimentary rock and which has 
seldom been drilled to any appreciable penetration, should be similar 
in many ways to the basement rocks of northern Minnesota and Wisconsin, 
or southern Missouri, which have exploitable mineral deposits. 

Major factors hampering exploration programs have been, 

— the economics of deep mining of ores from buried basement rock. 
Hence the characterization of many of the potential natural resources 
here as being "resources for the future". 


The basement rock is the crystalline complex of igneous and metamorphic 
rocks that underlies sediments and sedimentary rock, and extends down 
to the crustal base (Moho) . It is typically, and in the U.S. midcon- 
tinent, of Precambrian age, and represents part of an ancient continen- 
tal craton (shield) . 
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— the past absence of sufficient encouragement from remote methods — 
gravity, magnetics, seismic, satellite remote sensing — to pursue pros- 
pects in otherwise unprobed terrain. In the past decade there has been 
much new data gathered and areally compiled, by state geological sur- 
veys, federal agencies, and others, so that a regional picture of sub- 
surface geology can be developed in considerable detail. It is the 
regional geophysical surveying that allows "mapping", or at least in- 
terpretation, of the buried geology. Examples of this would be the U.S. 
gravity and magnetic anomaly maps (U.S. Geological Survey/Soc. of Explor. 
Geophysicists, 1982), and the magnetic anomaly map for the central U.S. 
(Hildenbrand et al, 1983). 

Interpretation of the long-wavelength Magsat anomaly map, in 
conjunction with the other more detailed geophysical and geological data 
sets, may aid in encouraging future exploration for resources. This might 
be particularly the case in reconnaissance exploration in frontier re- 
gions, that is where there has been little study to date or at depths 
that had not previously been considered. It is for such work that the 
satellite magnetic data could contribute in establishing the crustal 
framework for exploration strategy: the gross crustal structure and 

tectonic evolution, boundaries between terrenes of different geologic 
character, possibly the thermal regime in the crust, and anomalous zones 
in the lower crust which may have consequences for mineralization in the 
upper crust or overlying sediments. 

Two crustal tectonic features, and which can have major geophy- 
sical expression, are of special importance: 

— basins, formed by subsidence of crust locally and accompanying de- 
position of sedimentary rock. These can provide the source rocks, reser- 
voir rocks, and traps for the accumulation of hydrocarbons (oil and gas) . 

— major boundaries between crustal blocks. These could represent 
either compressional suture zones, where continental fragments (plates, 
or microplates) have collided and been welded together, or transverse/ 
shear faults where crustal blocks have moved sideways against one another, 
or tensional (failed) rift zones where continental blocks have tried to 
separate. These boundaries are zones of major crustal weakness, can be 
identified as deep structural lineaments, and can be the locus of much 



geologic "action" . The latter can involve enrichment and emplacement 
processes that can in turn lead to economic mineralization. This geo- 
logic activity would include deformation (folding, faulting, mountain- 
building), intrusion and plutons at depth with upward injection of dikes 
and sills, perhaps extrusion of magma at the Earth’s surface, circula- 
tion and invasion of hydrothermal fluids and minerals from depth up 
through the network of fractures, and reactivation later along the pre- 
existing zone of crustal weakness. 

For example, the paleorifting that resulted in the development 
of the Mississippi Embayment might be causally related, thermally and 
structurally (via uplift and subsidence) to the lead-zinc deposits in 
the "Mississippi Valley" region (Arkansas, Missouri, up to Illinois). 

Also, there are many faults in the Duluth complex, a crescent-shaped 
intrusive system 200 km long on the northwest shore of Lake Superior, 
and which developed as part of the CNARS/MGA paleorifting. These may 
have helped control the later intrusion by hydrothermal fluids to de- 
posit the copper/nickel sulphide orebodies at the base of the complex. 

Continental rifts — with their tendency for tensional and ther- 
mal thinning of the crust, and doming or graben development — can be 
associated with the development of some sedimentary basins. This is 
because early rifting, failed rifting, and later reactivation can all 
affect the formation of basins by controlling subsidence, uplift, and 
the depositional/sedimentary/structural character of the sediment package 
being deposited. For example, regional subsidence or uplift can affect 
deposition versus erosion, affect cycles of marine transgression and 
regression, and introduce facies changes. Such changes in depositional 
environment are very important for developing hydrocarbon accumulations. 
Faults can also propagate up through younger sediments, and affect basinal 
structures, from reactivated faults in the older underlying Precambrian 
basement . 

Thus for many aspects of upper-crustal geology, 

"... the Precambrian basement calls the shots . " 

Anonymous 

4th Internat. Conf . on Tectonics 

Aug. 1981 


It is in the interpretation of major structure, composition, 
and tectonic development of this crustal basement rock, that satellite 
magnetic anomaly data can hope to provide an assist. 

Let us consider now in more detail the resource prospects for 
the study area — the central Midcontinent. 

(a) Oil and gas 

To date, the potential for hydrocarbon resources in the northern 
midcontinent has not been considered to be high. Significant, but still 
fairly modest, exceptions are the oil and gas found in localized basins 
in Kansas and around the periphery of the Michigan Basin in lower Michi- 
gan (see Figure 25) . There have been, for example, no oil finds or pro- 
ducing wells in Iowa. 

Although it is not widely known, due to the proprietary nature 
of oil exploration, there has arisen in the past year or two a serious 
interest in the oil and gas potential of the Iowa region of the midcon- 
tinent. This activity comes to the surprise of some, and to the disbe- 
lief of others. The interest is in the form of geophysical (seismic) re- 
flection work, leasing of land, and plans for test drilling in the near 
future. While we are not involved in any of this exploration in what 
could be considered, geologically, to be a "frontier' 1 area as regards 
prospects, the state Geological Survey does monitor the activity. Such 
exploration probably has its impetus from a combination of: 

— thoughts, or hopes, that improved economics of oil and gas produc- 
tion will in the future allow exploitation of what might now be 
marginal reserves 

— the desire for domestic and more secure sources of oil and gas 

— interest in the geologic nature and development of the CNARS /Mid- 
continent Geophysical Anomaly and its associated structures. The 
latter would be particularly the sedimentary basin development 
along the flanks of the central horst, and on top of the horst 
block, and in the vicinity (e.g. the Forest City Basin, in south- 
west Iowa and northwest Missouri) . The location of the MGA in 
the Iowa region is shown in Figure 26, and geophysically inter- 
preted cross-sections showing the CNAB.S structure (horst and ba- 
sins) are given in Figures 42-44. Note that the latter models 
have a vertical exaggeration of 9X. 


The basins of greatest interest, but as of yet unknown prospects, 
are thus, 

i) the Forest City Basin (Figure 26) . This has over 1500 meters of 
Paleozoic sedimentary rocks. It is not clear if these sediments have been 
buried deeply enough to mature organic material to oil, or even if there 
was sufficient organic source material, or if the rocks (and oil) might be 
too old to have retained hydrocarbons to the present. If the organic mat- 
erial was present at depths not generally great enough to encounter temp- 
eratures suitable for maturing it to oil, a thermal pulse could have helped 
Such a change in the thermal regime, perhaps associated with major deep in- 
trusion or a paleorifting attempt, would raise the geothermal gradient and 
elevate isotherms for a period of time. Apart from the thermal rift tec- 
tonism of late-Precarabrian time here, there is some speculation that there 
may have been Cretaceous-age intrusion of some older basement rock upward, 
and injection of kimberlites, at least into eastern Kansas. Such a thermal 
event, with associated hydrothermal fluids, might have provided sufficient 
localized heating to mature hydrocarbons. 

ii) basins, structural or erosional in origin, on the top of the ma- 
fic central horst of the MGA. These would be relatively more accessible 
to drilling. Such basins are shown schematically in the hypothetical "mo- 
del" of Figure 29, and the gravity- and magnetics-calculated cross-sections 
of Figures 42-44. The latter suggest that basins (on these three profiles) 
could have widths up to 15 km, and depth extents up to 1500 meters or so. 
These basins contain Precambrian clastic sediments. The rationale for 
the potential for "Precambrian" oil or gas would be that there is hydro- 
carbon-rich oil shale in Precambrian sediments in northern Michigan (e.g. 
in the White Pine copper mine, on the Keweenaw peninsula) , and there may 
be the prospect for gas in these old rocks. The source for hydrocarbons 
along the central portion of the MGA might have been from algae and bac- 
teria, in sediments in ponded depressions along the subsided rift zone 
in Precambrian time. A modern analogy might be the algae similarly located 
in lakes and ponds along the East Africa rift system today. Other fac- 
tors to consider would be the thermal regime of the CNAK.S paleorift in 
late-Precambrian time, and subsequent geothermal history, to assist mat- 
uration of the shalloer organic material to oil, and also whether appro- 
priate trapping structures exist to retain oil or gas. 

iii) basins, flanking the central horst. As modelled (as in Figs. 
42-44) and by analogy to geologic sections exposed to the north in 
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Minnesota, these long downwarped troughs have widths up to 50 km, and 
depths from 5 to as much as 10 km. The sedimentary rocks are Precam- 
brian elastics. These basins have been virtually unsurveyed seismically, 
and are untested by deep drilling. Speculations on their hydrocarbon 
content must again consider whether they are too old to have retained 
oil and gas, and whether their greater depth (especially for depths 
over 5000 m.) would result in the higher temperatures having ’’cooked" 
oil to gas completely. 

(b) Mineral deposits and provinces 

Figure 45 shows the location of known copper, zinc, and lead 
deposits and provinces in the central midcontinent (from Tooker et al, 
1980) . "Provinces'* are areas where there are groups of economic depo- 
sits, and sub economic occurrences, and are thus favorable search areas. 

If mineral provinces can be causally related to major structural or tec- 
tonic features, and the latter can be identified using remote geophy- 
sical methods (including satellite magnetic mapping as an aid), then 
the initial exploration process can be improved. For example, from 
observation of the distribution of mineral occurrences in the continental 
U.S . , Tooker et al (1980) suggest that, 

— "copper provinces seem to be closely identified with characteristic 
geologic environments and host rocks in the two main types of 
crustal plates on the (U.S.) map — accreted oceanic and island arc, 
and old continent." 

— 2; i- n c provinces seem to be identified with old continental crustal 
rocks, especially uplifted regions where Precambrian rocks are ex- 
posed. The deposits are often localized in the sedimentary or 
other younger rocks covering or adjacent to the ancient crystalline 

basement rocks. "Stratabound zinc provinces, mostly of the Missis- 
* 

sippi Valley type , are spaced at Intervals along two broad north- 
east-trending zones that seem to parallel axes of uplifted and 
locally exposed basement rocks . Tf The most prominent zone lies in 
a band from Oklahoma/Arkansas to northern Illinois/lndiana; it 


i.e. in carbonate sedimentary rocks along the basin platform margins 
adjoining areas of uplift 
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lies roughly between provinces Z3 and Z9 on Figure 45, and is closely 
associated with faulted uplifted basement rocks. 

— lead provinces are found mostly in old continental crust (rather 
than accreted oceanic and island-arc terranes) , typically in Phanerozoic 
sedimentary rocks over Precambrian basement rocks and especially if the 
latter have been uplifted. 

The mineral resource localities shown in Figure 45 are listed 

below: 


Description of mineral deposits and provinces 
in the TJ. S. Midcontinent 


Number Location 


Type of deposit 


Cl 


C2 


C3 


C4 


i 

j 


C5 


Copper 

northern Michigan stratabound in (volcanogenic or metasedi- 

(in late-Precambrian mentary) sedimentary rock; disseminated 

terrane) deposits in conglomerate & basaltic rocks 

- has one active copper mine, plus other mines presently 
inactive, plus occurrences of unproven value 

Minnesota, Wisconsin, stratiform, in volcanogenic sediments, and 

n. Michigan maybe greenstones; massive sulphide 

(in mostly Precambrian deposits 

terrane) 

- southeast lobe has occurrences of unproven value (there 
are disseminations of copper and nickel in mafic rocks 
at the base of the Duluth Complex) ; region may have 
"large but unevaluated" resources of copper (Tooker et 
al, 1980) 

s. Wisconsin, n. Ill- hydrothermal, in sedimentary rock; 
inois ("Upper Missis- vein or replacement deposits 
sippi Valley lead 
district") 

- occurrences of unproven value 

southeast Missouri hydrothermal, in sedimentary rock; includes 

("Lead district") deposits in early-Paleozoic carbonate reef 

rocks which are peripheral to exposed late- 
Precambrian basement rock 

- has one active deposit (copper is a byproduct) plus 
other occurrences of unproven value 

Kansas, Oklahoma, Texas stratabound in sedimentary rock; disse- 

(in Permian sedimentary minated deposits 

basins) 

- occurrences of unproven value 



Number 


Location 


Type of deposit 


. . .Zinc. 

Most economic concentrations are as zinc sulphide or oxide, 
from ores of zinc, zinc/lead, or copper/zinc/lead. The zinc is believed 
to be concentrated by residual solutions from magmatic processes at depth, 
or by brines from basin sediments. "The (zinc) provinces in the Central 
Plains already have provided large production, and retain a high poten- 
tial for extensive future production .. .A potentially large but as yet 
unevaluated future resource (here) may be the zinc sulphides that fill 
fractures in coal beds in the Illinois and lowa/Kansas /Missouri (e.g. 
Forest City Basin) coal basins." The prominent zinc province that trends 
northeast, mentioned previously, lies between coal basins (localities 
Z3 and Z9 in Figure 45) . 

Zl northern Wisconsin stratabound; massive sulphides in volcano- 

genic rocks 

Z2 Wisconsin/Illinois/Iowa stratabound replacement and pore filling 
("Upper Mississippi in platform carbonates and breccia 

Valley district") 

- includes an active zinc mine 

Z3 Illinois stratabound or stratiform; in sedimentary 

(in coal basin) rocks 

- no known occurrences 

Z4 Illinois/Kentucky vein or pore filling, from hydrothermal 

("Hicks dome" district) solutions 

- at uplifted domal structure at the north end of the 
Mississippi embayment 

Z5 Missouri same as for Z2 

(southeast lead district/ 

Viburnum trend) 

- a major producing area; has several active mines with 
zinc as the primary or coproduct metal 

Z6 Kansas /Okla. /Missouri stratabound replacement in carbonates, 

("Tri-State district") or supergene (oxidized) enrichment, or 

residual deposits 

- has formerly-active zinc mines 

Z7 northern Arkansas same as for Z2 

Z8 Missouri same as for Z2 

("central lead district") 

Z9 lowa/Mis sour i/Kansas/ stratabound, disseminated in sedimentary 

Oklahoma basin rocks 

- in coal basin; unassessed localities 
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Number Location 


Type of deposit 


■( 

i 

I 


Lead 

Most of the U.S. production of lead comes from Mississippi Valley- 
type deposits in the Central Plains region (tlis Ozark uplift in southeast 
Missouri, see Fig. 25; and through Missouri). These are stratiform de- 
posits in vein fillings and pores in breccia, and as bedded replacement 
in Paleozoic platform carbonates that lie along the margins of uplifted 
regions . 


LX 


L2 


L3 


L4 


Wiseonsin/Iowa/Illinois 
Cupper Mississippi 
Valley) 


fissure veins and replacement in sedimen- 
tary rocks, generally related to intru- 
sive bodies; or stratiform or stratabound 
disseminations and pore fillings formed 
diagenetically in carbonates 


- includes an active mine with lead as a byproduct 


Missouri 

(central district) 

southeast Missouri 
(Viburnum trend) 


fissure veins and replacements in sedi- 
mentary rocks 


stratiform or stratabound disseminations 
and pore fillings in carbonates 

- includes several active mines and inactive ones; is the 
largest lead-mining producing district in the U.S. 

SE Kansas/SW Missouri/ 

NE Oklahoma 
(Tri-State district) 

- includes an inactive mine 


same as for LI 


L5 

northern Arkansas 

same 

as 

for L2 

L6 

s . Illinois /Kentucky 
(Hicks dome) 

same 

as 

for L2 

L7 

central Tennessee 
(Nashville dome) 

same 

as 

for LI 


It might be noted that the most magnetically-anomalous region 
of the mafic lower-crustal layer in the central U.S, is a north-south 
elongate zone that underlies much of the state of Missouri (see Figure 
46 here, an excerpt from Figure 38) . It has an interpreted magnetization 
intensity of 5-6 amp/m, compared to the U.S. continent lower-crustal 
average of 3.5 amp/m and a typical range of values of 2. 5-4.0 amp/m in 
the central U.S. as a whole. These magnetization values are calculated 
(Schnetzler and Allenby, 1983) by considering the magnetic moment strength 




i 
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Figure 46 . Magnetization intensity in the lower-crustal 
mafic layer, for central midcontinent (excerpt from 
Figure 38; Schnetzler and Allenby, 1983). Units are 
amp/m. 


of the source crust (as obtained by inversion of satellite magnetic 
anomaly data) and the thickness of the lower-crustal layer. The latter 
is assumed to be the source of the long-wavelength satellite anomalies. 
The resulting map shows the lateral variation in magnetization of the 
lower crust. The central-Missouri value in excess of 6 amp/m is about 
as high (i.e. as magnetically intense) as any lower-crustal rock that 
has yet been modelled for the earth (see Table of values, page 80 here). 
This major anomalously-r:agnetic zone is perhaps related to an ultrabasic 
intrusive complex or to a continental fragment of particularly favorable 
composition or metamorphic grade. In any event, it apparently underlies 
the region of the central midcontinent which has the greatest number. 


:o 



variety, and richness of metallic mineral deposits (see Figure 45) . 

This is central and southern Missouri, and extending slightly south 
into Arkansas. This area includes the following economic mineral de- 
posits and provinces: 

— Zinc Z5-Z8: includes the "Lead district/Viburnum trend", a. 

major zinc-producing area, in the southeast, and the "Tri-State 
district" in the southwest, and the "Central district" 

— Lead L2-L4: includes the Viburnum trend (the largest lead-pro- 

ducing district in the U.S.), and Tri-State district, and Central 
district 

-- Copper C4: in the Lead district, and has one of the two copper 

deposits being mined in the midcontinent at present 

One can, at this stage, only speculate on a possible genetic 
relationship between the deep-seated highly-magnetic mafic rock, as iden- 
tified by inversion of satellite magnetic data, and the relative abun- 
dance of metallic ore deposit localities and occurrences found in the 
overlying crustal rocks in this region. 

In addition to the above metallic mineral deposits, there are 
also economic iron-ore resources in the region (Hutchison, 1983) . These 
include the sedimentary iron formations (e.g. Biwabik-Gunflint) in nor- 
thern Minnesota and upper Michigan, and the Kiruna-type massive iron 
ore in southeast Missouri. The latter is at Pea Ridge, about 40 km 
north of the Viburnum trend’s lead deposits. The iron ore, primarily 
magnetite and of Precambrian age, is believed to be intrusive and gene- 
tically related to magmatic activity from below. The rhyolitic St. Fran- 
cois Mountain area at Pea Ridge has ten known iron occurrences and three 
major ore deposits. 
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MAGNETIC PROPERTIES FOR MAGNETIC MODELLING 


Magnetic modelling of crustal bodies to match observed anomalies , 
or inversion of anomalies to their causative sources, will be aided by 
knowledge of what plausible physical properties of rock types would be. 

The influence of crustal structural and petrologic factors has 
been discussed (pages 78-85) . These include: 

i) the thickness of the "magnetic lithosphere" that is magnetizable 

ii) the relative proportions of the crust that are (sialic) upper- 
crust, and (more mafic) lower-crust 

iii) relative contributions of induced and remanent components of 
magnetization 

iv) relative metamorphic grade of the crystalline basement rocks 

v) lateral geologic ihhomogeneity of crustal blocks or terranes 
(having, for example, different thicknesses, composition, age, 
structural texture, physical and magnetic properties) . 

In addition to these "bulk" characteristics, however, we should 
consider the magnetic mineralogy and associated magnetic properties of 
the rock types in the crust and upper mantle. This is because it is the 
magnetic minerals — and their relative abundance, composition of impuri- 
ties, chemical conditions (e.g. oxidizing, reducing, hydrating) during 
initial crystallization or later metamorphism, thermal history and cur- 
rent temperature, grain size — that determine the magnetization. For 
example, it is the Curie temperature (dependent on chemical composition 
of the magnetic mineral grains) which determines, in combination with the 
local geothermal gradient, the magnetic "bottom" of crustal sources. 

Magnetic mineralogy 

The most important minerals for rock magnetism are in the compo- 
sitional range bounded by the solid-solution series, ulvospinel-magnetite 
(titanomagnetite series), Fe 2 TiO^-FegO^, and ilmenite-hematite (ilmeno- 
hematite series), FeTiOg-FegO^. This is shown in Figure 47 . By far 
the dominant magnetic mineral in importance for magnetization is magne- 
tite, Fe^O^. The magnetic mineralogy, and associated magnetic and phy- 
sical properties, are given on the pages to follow (excerpt from Carmichael, 
1982; pages 246-255). 
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Figure 47 . Magnetic mineral compositions in 

ternary diagram FeO-TiO^-Fe.CU . Solid-solution 
series "I" is titanomagnetiref "II" is ilmeno- 
hematite. (from Carmichael, 1982) 
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MAGNETIC MINERALOGY, CRYSTALLINE AND MAGNETIC 

PROPERTIES 

The magnetization of rocks is retained in certain magnetic minerals. These minerals 
and their physicochemical state control the intensity and stability over time of the 
remanent and induced magnetization. Proper geological interpretation, as for mag- 
netic prospecting or paleomagnetism, of the magnetic character of rock samples de- 
pends on knowledge of the mineralogic, structural, magnetic, and mechanical proper- 
ties of the minerals present. 

Representative values for important properties are given here. The values are “typi- 
cal”. There is sometimes considerable discrepancy between published values; this is 
often due to the differing conditions and composition of the “natural” minerals and 
rocks studied. 


Important Minerals in Rock Magnetism 

The major minerals having magnetization, or being of interest in studies of magnetic 
minerology are 


Iron oxides 

titan □ magnetite scries: ulvospinel-magnctitc 

xFc 2 Ti0 4 .(I - x)FciO* 


ilmcnohematite series: ilmenitc-hcmatitc 

yFeTiOj.<l - y)FciOj 


maghemite 

martitc 

geothite 


lepidocrocite 

akaganeitc 

Sulfides 

pyrrhotitc series: 


yFe 2 Oj 

oFciOj 

aFcOOH (the most common of the natural 
hydrous ferric oxides, ne. T 
"limonite” FejOj.HiO) 

yFcOOH 

/JFeOOH 

troiliie-pyrrhotite 

yFeS.U-yJFc.-gS 


pyrite FeSi 

marcasite FeS 2 

mackinawite FeS 

Carbonates 

siderite FcCOj 

magnesite MgCOj 


Iron Oxides 

The most important and common rock-forming magnetic minerals for the study of 
rock magnetism and paleomagnetism are the iron oxides. They can be represented in 
the ternary diagram FeO-TiOi-Fe^Oj shown in Figure 9. There are three main solid 
solution series:. 


1 . Titanomagnetite series — cubic structure (inverse spinel) 

ulvospinel • magnetite 
xFe, 1 + Ti” 04 .(I-x)Fe 3 *(Fe 1 *Fe 3 *) 0 4 

2, Ilmenohematite series — hexagonal/rhombohedral structure 

ilmenite - hematite 
yFe 2 *Ti 4 *Q s .( l-y)Fei J *Oj 


^ ^ .- 


origiwm. PJ*ew 
OF POOR QUALITY 


RUTILE 
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FIGURE 9* Compositions in ternary diagram Fe0-Ti0j-Fe 2 0,. Solid-sal mion series are i — litanomag- 
netite. It — ilmenohemattte. III — pseud obrookite. Dotted lines represent direction of oxidation. 


[ 

\ 


. - i 




3 . Pseudobrookite scries — orthorhombic structure 

zFe 3 34 Ti 44 Os-(l -z)Fe 24 TI 3 44 O s 

Possible impurity phases are MnO, MgO, A1*0 3 , and V1O3. The solid solution region 
of greatest interest in nature is shaded in Figure 9. The dotted lines are lines of constant 
Fc:Ti ratio and represent trend of oxidation/reduction. The oxidation increases to the 
right. 

Ti0 2 has three polymorphs. It is called rutile with a tetragonal structure {c axis = 
0.6 A) or anatase if tetragonal (c axis = 1.8 A), and brookite if orthorhombic. 

The diagram is drawn on the basis of molecular ratios. Thus, 


FeTiOj 

FeO • TiOj 

r 

Fej TiO, = 

- FeO * — TiO, 
3 3 1 

i 

4, 

Fe 3 O. = 

FeO - Fe s 0 3 

f 

t 


Figure JO shows crystallographic structure and type of magnetism in the ternary 
diagram. In the ilmenohematite series, the magnetism varies with the composition as 

y~I antiferromagnetic {Fe and Ti both occupying all cation layers equally) 

0,45 v 0.95 Ferrimagnetic (ordered state of Fe and Ti t with Ti ions occupying every 

second cation layer perpendicular to the “c" axis) 

0 y 0.45 Antiferromagnetic (parasitic (spin-canted) antiferromagnetism; Fe and Ti 

occupying all cation layers, t.e.. disordered state) 

For 0,45 < y < 0.60, synthetic specimens show self-reversal of magnetization. 

Summarized in Figure II are some magnetic properties of the minerals. The Curie 
temperatures for the end-members of the series are shown in degrees centigrade. They 
vary uniformly from one end-member to the other. The increase in lattice parameter 
“a” is shown by arrows, likewise for saturation magnetization J„ and the constants 
K for magnetocrystalline anisotropy and A for magnetostriction. An auxiliary series is 
magnetite — hausmannite, Cl-x)FeAp 4 -xMn 3 0 4 . For 0 < x < 0.6, it is cubic structure. 
For 0.6 < x < 1 it is tetragonal. 
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FIGURE 10. Structure and magnetism (at room temperature) in ternary diagram. 



FIGURE 11. Magnetic properties in the ternary diagram. Curie temperatures in °C, “a" is cell dimension. 
Arrows indicate direction of increasing parameters. 


The relationships between the main iron oxides are shown in Figure 12. The temper- 
atures given are approximate. Values vary widely in the literature, depending on the 
sample condition and experimental conditions such as ambient atmosphere (air, vac- 
uum, etc.). The structural conversion from maghemite to hematite is pressure-depend- 
ent. The goethite to hematite conversion could occur during consolidation of sedi- 
ments, as 2FeOOh Fe,Oj + H 2 0. The magnetite to hematite oxidation could occur 
during initial cooling of igneous rock, at temperatures of about 600 to 1000°C. The 
magnetite to maghemite conversion could occur as low-temperature oxidation, as in 
late cooling or weathering. 

The conditions for precipitation of different iron minerals are shown in Figure 13. 
The stability fields are outlined by the parameters Eh (volts) representing oxidizing 
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FIGURE 12. Conversion of iron-oxide magnetic minerals. Temperatures are approximate. 




FIGURE 13. Conditions for precipitation of different iron minerals. (From Gass. 1., « al. t Understanding 

the Earth, MIT Press. Cambridge, M A, 1971 . 179. Article by Watson. E. K.) \ 

) 

potential, and pH (concentration of H ions) representing acidity/alkalinity. The sta- j 

bility fields also depend on the concentration of components (e.g., Fe) in solution. j 

ABCD Is the range of naturally occurring solutions. The upper hatched area is the 
stability field of orFeiOa. j 


The main minerals will now be described briefly. These oxides and oxyhydrates can 
be described in terms of stacked “almost-close-packed” layers of oxygen atoms with 
iron atoms behavine as interstitials. These irons have an atomic radius such that thev 
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will not quite fit into the available sites of the oxygen lattice without distorting it. The 
distances between the oxygen planes in all the structures are about the same, ranging 
from 2.3 to 2.5 A. The stacking sequences are characteristic of close-packed lattices, 

i.e., cither the -ab-ab-ab- (hexagonal close packed) type, or the -abc-abc- (face centered 
cubic) type. 

In the following, Z = number of formula units per unit cell, T* = temperature of 
magnetocrystalline transition, T c = Curie temperature, T* = Neel temperature, E = 
Young’s modulus of elasticity, and G = shear modulus of elasticity (modulus of rig- 
idity). 

Values are at room temperature unless otherwise noted. 

Magnetite 

1. Chemical formula: Fej0 4 ; Z = 8; percent by weight: Fe 2 *(24.1), Fe 3 *(48.3), 
0 2 (27.6). 

2. Structure: cubic, spinel of inverse type (see Figure 14). Stacking of hexagonal 
close packed oxygen planes with an -abc-abc- sequence along the <1 1 1> direc- 
tion. The general spinel structure is A* 2 B 2 °0 4 ’ 2 , where A occupies a tetrahedral 
site and B octahedral sites. In the unit cell, there are 32 oxygens, 64 available A 
sites, and 32 available B sites. The "normal” spinel has 8 A’s in tetrahedral sites 
and 16 B’s in octahedral sites. The "inverse” spinel has 8 B’s in tetrahedral sites, 
and 8 A’s and 8 B’s in octahedral sites. Cations occupy layers between two adja- 
cent oxygen planes; successive layers have cations either in the octahedral (B) 
sites between individual oxygens, or one third in the octahedral and two thirds 
in the tetrahedral sites. The inverse spinel arrangement for magnetite is 
Fe^Fe^Fe 3 *) 0 4 2 ’. The occupied tetrahedral sites form a diamond-type lattice. 
The octahedral sites have sixfold coordination, and the tetrahedral fourfold. The 
oxygen atoms form a face-centered-cubic lattice. The structure converts from 
inverse spinel to orthorhombic below about -155°C. The transition is abrupt for 
synthetic magnetite crystals, but can be spread over about 10° for natural crys- 
tals. The presence of impurities lowers the transition temperature. At high pres- 
sures (about 225 to 240 kbar), Fe 3 0 4 undergoes a phase transition to a monoclinic 
phase with density about 6.4 gm/cm 3 . Cell dimensions (cubic phase): = 8.394 
A; distance between (111} oxygen planes ~ 2.9 A; distance between Fe and O, 
octahedral site — 2.06 A, tetrahedral site — 1.87 A; angle between O-Fe-O, oc- 
tahedral site — 88.1°, 90°, 91.9°, tetrahedral site — 109.5°. 

3. Mineral characteristics: crystals are generally octahedral form {111}, but may oc- 
cur as cubic {100} or dodecahedral {110}; twinning and parting on {1 1 1}; hardness 
about 6 on Moh’s scale; density = 5.18 gm/cm 3 . 

4. Magnetic properties: ferrimagnetic at room temperature; saturation magnetiza- 
tion J, ~ 98 emu/gm (at 0°K); ~ 92 emu/gm (at room temperature) = 480 emu/ 
cm 3 ; critical single-domain size about 0.1-1 pm; magnetic anisotropy — magne- 
tocrystalline transition at T* = -140°C; 


T|c < T < T c 

-155° < T < T 

T < —155° 

K x ~ -1.35 X 
K a - -0.48 X 


-- easy axis <11 1> 
hard axis <100> 
k -- easy axis <100> 
hard axis <11 1> 

— easy axis orthorhombic 
10* ergs/cm 3 
10 s ergs/cm 3 


K, negative 
K t positive 


T* transition suppressed if grain size too small (less than about 0.1 pm); Aioo ~ 
-20 x 10’* cm/cm; A Il0 ~ 60 x 10’*; A tI1 ~ 78 x 10’*; As ~ 40 x 10 ’*; T c = 575°C; 
domain (Bloch) walls thickness about 500 to 1500 A and energy about 1 erg/cm 3 . 
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• OCTAHEDRAL (B) SITE 
® TETRAHEDRAL (A) SITE 
o OXYGEN IONS 




FIGURE 14. Spinel structure of magnetite. Top cube is 1/8 of unit cell. 


5. Mechanical and other properties: E ~ 1.6 x 10 4 kg/cm 2 : G ~ 0.5 * 10 4 kg/cm 2 ; 
Poisson’s ratio o *v» 0.3; slip plane {111}, slip direction <110> for dislocations; 
ultimate compressive strength (crystal) ~ 1000 kg/cm 2 ; electrical resistivity q ~ 
lO^-lO 2 ohm-meter (places it in the range of semiconductors, 10' 5 to 10* ohm- 
meters. Other ferrites are generally much less conductive.) 

Ulvospind 

1 . Chemical formula: Fe 2 Ti0 4 ; 

percent by weight: Fe 2 *(50), Ti 4 *(21.4), O a ’(28.6) 

2. Structure: cubic, spinel of inverse type (see Figure 14). Fe 2 * and Ti 4 * replace 2Fe 2 * 
in magnetite. Cell dimension: a ~ 8.5 A 

3. Mineral characteristics: occurs as intergrowth on { 100} faces of magnetite; density 
about 4.8g/cm\ 

4. Magnetic properties: paramagnetic at room temperature: T« ~ -150°C; weakly 
ferrimagnetic below T e (should be antiferromagnetic, in theory). 

Maghemite 

1 . Chemical formula: yFe 2 0 3 

2. Structure: a polymorph of hematite (crFe 2 0 3 ); cubic, defect spinel of inverse type 

(see Figure 14). Isostructural with magnetite except 1/9 of cations missing. These 
vacancies are distributed either randomly throughout octahedral and tetrahedral 
sites, or in octahedral sites alone, as in 4 Fe 2 3 *0 3 3[Fe 2 *0.Fe S /3 , *.Di' 3 0 3 ] 

where □ represents a vacancy. Cell dimension: a~8.35 A; transition to aFc 2 0 3 
at T~250 to 500°C; distance^ between Fe and O, octahedral site — 2.05 A, tetra- 
hedral site — 1.86 A 


( 


i 
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3. Mineral characteristics: may form as oxidation product of magnetite (3Fe 2 * 
2Fe J *), or converted from lepidocrocite (see Figure 12). The tetrahedral sites are 
less available for oxidation than the octahedral sites, since the former are cova- 
lently bonded and the latter ionically bonded; Density = 4.88 g/cm\ 

4. Magnetic properties: ferrimagnetic; J, - 83.5 emu/g (407 emu/cm J ); T c ~ 675 
to 750°C. (depends on structural state). 

Wustite 

1 . Chemical formula: FeO 

2. Structure: cubic, NaCl type; stacking of hexagonal close-packed oxygen planes 
with -abc-abc- sequence along the <1 1 1> direction. Lattice is like two interpene- 
trating f.c.c. lattices, one of oxygen and the other of irons. Defect structure, 
deficient m Fe, i.e., F *Q. 83H)*95°; cell dimension: a~4.30 A; distance between 
Fe and O — 2.15 A; angle between O-Fe-O — 90°. 

Hematite 

1. Chemical formula: crFe 2 0 3 ; percent by weight: Fe 2 *{70), O 2_ (30); Z = 6 for hex- 
agonal unit cell; Z = 2 for rhombohedral unit cell. 

2. Structure: hexagonal system, trigonal subsystem, rhombohedral class (see Figure 
15). Stacking of slightly distorted hexagonal ^dose-packed oxygen planes with 
-ab-ab-ab- sequence, along the [0001] direction, corundum-type structure. 
Groups of three oxygen atoms form a common face of two neighboring octa- 
hedra with a face parallel to {0001}. The stoichiometric number of iron atoms 
occupy octahedral interstices (sixfold coordination) between oxygen planes. They 
form FeiOj groups in the form of a trigonal dipyramid. Two thirds of the avail- 
able cation positions are filled, with each Fe being surrounded by six oxygen in 
a near-octahedron. Cell dimensions (refer to Figure 15): hexagonal unit cell — a 
= 5.035 A; b = 13.75 A; rhombohedral unit cell — a rk - 5.426 A; a - 55° 16'; 
distance between oxygen planes about 2.16 A; distance between Fe and O — r T 

- r 2 = r 3 = 2.08 A, d, = d 2 = d 3 = 1.95 A; angle between O-Fe-O (angle 

defined by bonds listed) — r,r 2 = r 2 r 3 = r,r 3 = 77°; = r 2 d 3 = r 3 di = 86.1°; 

r 3 d 2 = r 2 d, = r,d 3 = 90.8°; d»d 2 = d 2 d 3 =* d,d 3 - 102.5°. 

3. Mineral characteristics: crystals have form of positive rhombohedron {1011}, 
negative rhombohedron {0ll2}, or pinacoid {0001}; parting on {0001} or rhom- 
bohedral plane {0112} due to twinning. Twinning on {0001} as penetration twins 
and on {01 12} usually lamellar. Hardness 5-6; density = 5.27 gm/cm\ 

4. Magnetic properties: each sheet of Fe is ferromagnetic, but sheets are coupled 
antiferromagnetically, i.e., two antiferromagnetic sublattices. There is a system- 
atic deviation from oppositely-directed spin configuration, resulting in weak 
spin-canted (anti)ferromagnetism. Magnetic anisotropy: magnetocrystalline tran- 
sition T k at about -23°C {^13°C for synthetic material); -23°<T<T C — atomic 
moments in (0001) plane. Weak ferromagnetism. Fe atoms on trigonal axis of 
rhombohedral unit cell have spins parallel in groups of two, pointing in direction 
of “a” axis of hexagonal unit cell. Easy axis is <101 0>, hard axis is [0001 1. 
T<-23° — atomic moments point along trigonal axis. Perfectly anti ferromag- 
netic (no net magnetization). Successive spins directed oppositely. Easy axis is 
[0001]. Transition suppressed for grain size sufficiently small. T c <T<725°C(T*) 

— anti ferromagnetic. Has a hard isotropic stable magnetization which endures 
until T^ (Neel temperature). Probably associated with unbalanced spins due to 
lattice defects or impurities, and is found in natural, not synthetic, material. J* 
~ 0.45 emu/gm (2.4 emu/cm J ); decreases with grain size; T tf (weak ferromagne- 
tism disappears) = 675°C; T N (antiferromagnetism disappears) = 725°C. Criti- 
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cal single domain size, about 10 to 100 pm; Am'W.3 x 10-* cm/cm (i.e., normal 
to (0001)); A,! 2 ~ 8 x 10-* cm/cm (i.e., in (0001) plane). Both A’s become negative 
below region of magnetocrystalline transition. 

5. Mechanical properties: E ~ 2 x JO* kg/cm 2 ; G ^ 0.8 x io* kg/cm 1 ; Poisson’s 
ratio a ~ 0.27. 

Ilmenite 

1. Chemical formula; FeTiOj; percent by weight: Fe 2 *(36.8), Ti**(31.6), 0*'(31.6); 
Z = 6. 

2. Structure: hexagonal system, trigonal subsystem, rhombohedral class. Layers of 
Ti separated by Fe layers of alternating magnetic polarity. Cell dimensions: hex- 
agonal unit cell — a = 5.08 A; c *= 14.13 A; rhombohedral unit cell — a r * = 
5.52 A; a = 54°51\ 

3. Mineral characteristics: twinning on {0001}, parting on {Q00l},{0112}; hardness 
5-6; density = 4.78 gm/cm*; may occur as lamellae on octahedral {1 1 1} of Fe»0«. 

4. Magnetic properties: paramagnetic at room temperature; T W ~-210°C; T<T *, 
antiferromagnetic. 

Martite 

1 . Chemical formula: a-FejOj 

2. Structure: polymorphous after magnetite. Has cubic spinel outer form but rhom- 
bohedral internal structure. 

Iron Oxyhydroxides 

(Minerals of the family of hydrous ferric oxides, generally designated as Limonite) 

Goethite 

1 . Chemical formula: orFeOOH, or 2HFe0 3 ; Z = 4. 

2. Structure: orthorhombic, corresponds to hematite. Stacking of hexagonal close- 
packed oxygen planes with -ab-ab-ab- sequence along [001], Iron atoms occupy 
only the octahedral positions. Converts to hematite (see Figure 15). Cell dimen- 
sions: a = 4.60 A; b = 9.95 A; c = 3.02 A. 

4. Magnetic properties: antiferromagnetic at room temperature, but may have some 
stable remanence due to spin-canting or unbalanced spins. Fe atoms alternate 
direction of spin along axis of magnetization, [001]. T»~120°C. 

Lepidocrocile 

1 . Chemical formula: yFeOOH, or 2(FeO.OH); Z = 4. 

2. Structure: orthorhombic. Stacking of oxygen-hydrogen planes with an -abc-abc- 
abc- sequence along <05 1> direction of the orthorhombic unit cell (corresponds 
to <1 1 1> direction of a distorted cubic arrangement). Hexagonal packing of ox- 
ygen atoms in a sheet is not regular. Iron atoms occupy only octahedral sites. 
Each H atom is associated with an O atom, forming a discrete hydroxyl group. 
Cell dimensions: a = 3.06 A; b = 12.4 A; c = 3.87 A. 

3. Mineral characteristics: {010} cleavage. 

Akagancite 

1. Chemical formula: /JFeOOH; Z = 8. 

2. Structure: tetragonal. 

4. Magnetic properties: T* from -160°C to 20°C. 
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FIGURE 15. Hematite structure. Top — rhombohedron in hexagonal system: Left — schematic of Fe-O 
configuration; Right — rhombohedron of FeiOi units. 


Sulfides 

The main magnetic sulfides of interest to rock magnetism are in the troilite-pyrrho- 
tite series, troilite-pyrrhotite, yFeS.(l - y)Fe,- f S. Troilite is hexagonal in structure, 
witn Fc «nd S atom layers alternating. It has a niccolite (NiAs) structure (see Figure 
16). It is antiferromagnetic, with a Neel temperature of 320°C. The tetragonal form 
of FeS is termed mackinawite. For x = Vi % or pyrite (FeS 2 ), the structure is cubic and 
the crystal is paramagnetic at room temperature. 

Pyrrhotite 

1 . Chemical formula: Fei-,S, where 0<x<0. 125. 

Z = 2 

2. Structure: hexagonal system, rhombohedral class. In Fei-.S, x = 0 - hexagonal; 
0<x<0.07 - mixed structure; 0.07<x<0.1 - hexagonal; 0.1<x<.125 - monoclinic. 
When hexagonal, it has a defect NiAs structure (see Figure 16). It is related to 
FeS in that Fe 2 * replaced Fe 2 *, leaving some cation vacancies. Cell dimensions 
(forx = 0.115); a = 3.446 A; c = 5.848 A. 

3. Mineral characteristics: crystals usually tabular on (0001); hardness 3.5 to 4.5; 
density ~ 4.6 gm/cm\ 

4. Magnetic properties: 0< x <0.09 - antiferromagnetic; 0.09< x <0.14 - ferrimag- 
netic; J, ~ 18 emu/gm (83 emu/cm J ) for monoclinic phase; this maximum J, 
occurs for FcrS#, or Fe 0 .. 7 S S (i.e., x = 0.125); T c ^ 320°C; hard axis of magnet- 
ization [0001]. 

Carbonates 

Siderite 

1. Chemical formula: FeC0 3 ; 

Z = 6. 
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FIGURE 16. Hexagonal (NiAs, niccolite) structure for pyrrhotite. 
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FIGURE 17. Variation of magnetism with temperature, for magnetic minerals and ferromagnetics. 


2. Structure: hexagonal system, trigonal subsystem, rhombohedral class. Fe atoms 
are at points of a face-centered rhombohedron, COj groups midway between 
Fe’s. COj groups are trigonal, planar parallel to [0001], Cell dimensions: a = 
4.72 A; c = 15.46 A. 

3. Magnetic properties: antiferromagnetic; T* ~ -230°C; magnetization along 

[ 0001 ]. 

4. Mineral characteristics: crystals commonly rhombohedral (1011}, sometimes 
{0ll2) t {0221], {4041], or tabular on {0001}; cleavage {10ll}; hardness 4; density 
= 3.96gm/cm\ 

Figure 17 shows the variation of magnetic properties with temperature for magnetic 
minerals and ferromagnetics. 
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Values of Curie temperatures and magnetization, and their varia- 
tion with mineral composition, follow (as "Table 7" and "Figures 18, 19, 
and 20", on excerpt pages 256-260) . This is then followed by magnetic 
properties of rocks — remanent magnetization, susceptibility, and Koenigs- 
berger ratios— in "Table 16" (pages 272-279, Carmichael 1982). 

Experimental study of the variation of magnetic susceptibility, 
with pressures of upper-crustal significance (up to 2.7 kbars) on magne- 
tite and rocks, is given in Nulman, Shapiro, Maksimovskikh, Ivanov, Kim 
and Carmichael (1978) . Additional information on the more chemical as- 
pects of magnetic mineralogy is given in Haggerty (1979) . 

An uncertainty with regard to magnetic properties in situ at 
depth, as in the lower crust, is the likely Curie-temperature there. 

This will depend oil the compositional state of the titanomagnetite . For 
example, titanium-rich titanomagnetite has a low T c ("Fig. 11", p. 122). 
However, if there are oxidizing conditions (high oxygen fugacity) and 
high temperatures, as in the lower crust or in regional thermal meta- 
morphism at depth, then there could be exsolution within the titanomag- 
netite grains to form ulvospinel and magnetite, and oxidation toward il- 
menite. The ulvospinel would be non-magnetic, the ilmenite (of, say, 
intermediate composition) would have a low T c and small magnetization, 
but the magnetite would have a high (575°C.) and strong magnetization. 

An additional possibility for a contributor to the magnetization 
of the lower crust, and perhaps even the upper mantle, could arise from 
the serpentinization of deep mafic rocks (Haggerty, 1979) . The suggestion 
is that the olivines and sulphides would undergo reduction, yielding al- 
loys of iron/nickel/cobalt/copper. These would be appreciably magnetic, 
and would also have much higher Curie-temperatures (600-1 100° C.) than the 
usual iron oxide minerals. Their ability to retain magnetization would 
thus extend to greater depths. 
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Table 7 

J s AND T c /T n OF MAGNETIC MINERALS, ROCKS, AND METALS 




Density 

TcjTn 




Material 

Composition 

Notes (gm/cm 3 ) 

(°C) 

(emu/cm 3 ) 

(emu/gm) 

Ref. 



a b 

C 




Minerals 







Magnetite 

Fe,0 4 

Ferrimagnetic 







(pure) 5.18 
(natural) 5 -5 A 

578 

480 d 

92.3{24°C) 

10 





S10 d 

98.2{0°K) 

10 





471 


10 




580 



63 




575 

480 





0.6% Ti 


410-430 


16 



0.06% Ti 


480-500 


16 

Titanomagnctite, x * 0*55 

210 



37 

Ulvospxnel 

FejTi0 4 

antiferro- 4.8 

magnetic 

-153 



63 

Maghemite 

•yFcjOj 

Ferrimagnetic 4.88 

675 e 

407 d 

83.5(24°C) 

10 





417 


10 





545-675 e 


63 





750 


73 

Hematite 

aFe*O a 

Spin-canted (anti) 
ferromagnetic or 
ferrimagnetic 







5.26 

675 

2.6 

0.5(24°C) 

10 



(Synthetic) 

670 

2.1 

0.39(24° C) 

10 




680 

2.6 

0.5 

63 




680 

2.2 

0.42 

8 



•Tefcct^parasitic 

aniiferromagnetic f 

725 




llmenohematite. 

5?* 0.75 







max. ferrimagnetism 
in series 

-10 


19.5 

34 

Hmenite 

FeTIO, 

Antiferro* 4J2 

-205 


-0 

74 



magnetic 

-216 



63 

Wustite 

FeO 

Antiferro- 

-87 



74 



magnetic 

-83 



10 

Goethite 

aFeOOH 

Antiferro- 4.27 

magnetic 

120 



63 

Akaganeite 

pFeOOH 

Antiferro- - 

magnetic 

196 to 23 



63 

Magnesia- 

MgFc 1 0 4 

Ferrimagnetic® 4.18 



24.5 

63 

ferrite 



310 

1 LG 


10 





no 

23 

74 




440 

no 

25 

4 





140 

29(0° K) 

4 




tjp. .. *$■ . *-’4 . 
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Table 7 (continued) 

J s AND T c /T n OF MAGNETIC MINERALS, ROCKS, AND METALS 





Density 

T C iT N 




Material 

Composition 

Notes 

(gm/cm 3 ) 

<°c> 

(emu/cm 3 ) 

(emu/gm) 

Ref* 



a 

b 

c 




Jacobsite 

MnFe 2 0 4 

Ferrimagnetic 8 

4.95 

300 

416 

84 

63 






400 

81 

4 






408 


10 






560 

I12fO°K) 

4 

Chromite 

FeCr a 0 4 

Ferrimagnetic 8 

4.5-5. 1 

— 1 85 



10. 63 

Trevorite 

NiFe 2 0 4 

Ferrimagnetic 8 

5,35 

585 

270 

51 

63,74 






300 

56(0°K) 

4 






267 


10 

Franklinite 

ZnFe 2 0 4 

Ferrimagnetic 8 

5. 1 -5.3 

-258 



74 





-264 



10 

Fayalite 

Fe 2 Si0 4 

Antiferro- 

magnetic 


-147 



63 

Pyroxene 

FcSiOj 

Antiferro- 

magnetic 


-233 



63 

Pyrrhotite 

F C ,. X S 

Ferrimagnetic 

4.6 

300- 

325 62 

13.5 d 

10 





320 

90 d 

19.5 max. 
at Fe-jS, 

63 





300 



8 

Trailite 

FeS 

Antiferro- 

4.83 

320 



74 



magnetic 


340 



10 

Pyrclusite 

MnO a 

Anti ferro- 
magnetic 

5.06 

-189 



10,63 

Chromium 

CrO a 

Ferrimagnetic 


119 

515 


74 

dioxide 








Siderite 

FeC0 3 

Antiferro- 

magnetic 

3.96 

-233 



63 

Rhodochro- 

MnCo 3 

Antiferro- 

3.70 

-243 



63 

site 


magnetic 






Metals 








Iron 

Fe 

Ferromagnetic 

7.87 

770 

1714 

218 

2 






1760 

224{0°K) 

2 

Nickel 

Ni 

Ferromagnetic 

8.9 

358 



38 






485 

54 

2,38 






510 

57fO°K) 

2,38 





360 



4 





370 



82 

Cobalt 

Co 

Ferromagnetic 

8,85 

1120 

1420 

160 

2 






1445 

163{0°K) 

2 

High permeability metals: 







4% Si-Fe 




690 



74 

45 Permalloy 




440 



74 

Mumetal 




400 



74 

Supermalloy 




400 



74 
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Table 7 (continued) 

J s AND T c /T n OF MAGNETIC MINERALS, ROCKS, AND METALS 





's 


Density 

Material Composition Notes (gm/cm J ) 

a b 

Tc>T[sf 

(°C) 

c 

(emu/cm 3 ) (emu/gm) 

Ref. 

2V Permendure 


980 


74 

Rocks 





Sandstones 

Red; Britain; 12 sites, 
Triassic and Devonian 


(0.2— 1.2) 
X 10^ 2 

75,20 

Basalt 

Titanomagnetite x = 0.6 

200-400 


37 

Magnetite ore 

Lodestone (Fe^O*) 


50-80 

66 

Lunar rocks 





Soils 

8 samples 


G.9-1.5 

23 

Breccia 

16 samples 


0.05-2 

28 


Anorthosite breccia 

765 

0.145 

44 

Anorthosite 

2 samples 


0.7- 2*6 

28 

Basalt 


760 

0.2-2.2 

44 


10 samples 


0.1-2 

28 


3 samples 


Q.1-3[Q‘K) 

28 

Gabbro 

2 samples 


Q.2-0.7 

28 

Igneous rocks 


760-790 


28 

Breccia and tines 


745-790 


28 

Meteorites 





4 Iron meteorites 

(ie., Ni-Fe) 

GQQ-780 


32 


3 Magnetic state is at temperatures below T c or T N ; would be paramagnetic above* 

^ Used to convert emu/cm 3 to emu/gm or vice versa; density values from Handbook of Materials Sci- 
ence, VoL I, Lynch, C. T. t Ed., CRC Press, Boca Raton, Fla., 1974, p. 235, and Vol. III, p. 184. 
c T c is Curie temperature, for fcrrimagnetic or Ferromagnetic materials; T N is Neel temperature for 
antiferromagrierics. 
d Calculated, from other J s units. 
e May have been converted to hematite. 

f Present in natural specimens, due to effect of lattice imperfections. 

* Spinel structure, ic.. A* 3 B^ 3 0 4 


Variation of J,, T e , and Cell Dimension With Chemical Composition 
The titanomagnetite series Is ulvospinel-magnetite, xFejTi0 4 .(l“x)Fe 3 0 4 . Figure 18 
shows the variation of Curie temperature (T c ) and cell dimension (cubic lattice param- 
eter, i.c*, unit cell size, “a”) with ea x” in the above series. 

The relationships for titanomagnetite can be represented by 

cell parameter, a - 8.395 + 0.l35x A 
Le.,a “ 8.395 A for magnetite (x = 0) 
a * 8.53 for ulvospinel (x = 1) 

Curie temperature, T c = 575 - 725x + 600x (14 - xHI - x) °C. 

Le., T c = 575 "C for magnetite 
= -150 for ulvcspinel 
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FIGURE 18. Variation of T, and cell parameter “a” in titanomagnetite series. 1 




FIGURE 19. Variation of Curie temperature (T«) and cell parameter "a” in FeO-TiOi-FeiOj ternary dia- 
gram. (A) for T„ in °C, (B) for “a'\ in A. 


Figure 19 shows the variation of T, and “a” in the compositional range of interest 
in the Fe0-Ti0 2 -Fe 2 0 3 ternary diagram. 

The variation of saturation magnetization, 3 „ for the titanomagnetite series at room 
temperature is given by J. — 92(1 - x) — 42x(l - x) eimi/gm, i.e., 3 . = 92 emu/gm 
for magnetite, J, = 0 for ulvospinel. 

Table 8 shows a typical experimental determination of the variation of i, with com- 
position in the titanomagnetite series. 
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Fojpa F0T1O3 

Mole fraction of FeTiO s 


FIGURE 20* Variation of T, and rhombohedral cell parameter 
“a**” in ilmenohcmaiite series** 


Tables 

VARIATION OF J, WITH 
COMPOSITION IN THE 
TITANOMAGNETITE AND 
ILMENOHEMATITE SERIES* 2 34 


Titinomagnetite Iimenohematitc 

xFe,Ti0 4 , (l-x)FejO« yFcTiO*.(l--y)FeiO* 



J.(emu/ 


J, (emu/ 

X 

cm) 

y 

gm) 

0 

93 

0 

0.5* 

0.04 

90 

0.25 

0.6* 

0.10 

82— 

0.45 

0.6* 

0.18 

73 

0.5 

2.9* 

0.31 

59 

0.6 

9*2* 

0.56 

29 

0.7 

17.4* 

0.68 

15 

0.75 

19.5* 



0.8 

18.2* 4 



0.9 

5.3* 4 



0.95 

l.fi* 4 



1.0 

0* 4 


* Spin-canted (anti) ferromagnetic. 

* Ferrimagnetic. 

* Antiferromagnetic. 

4 For x » 0.75, T* is below room temperature, so the 
material is paramagnetic at room temperature. 

The ilmenohematite series is ilmenite-hematite, yFeTiOj. (l-y)Fe t O». 

Figure 20 shows the variation of T £ and ceil parameter a,* with “y” in the series. 
Relationships can be represented by Curie temperature T c = 680 - 890y °C, i.e., T c 
= 680 °C for hematite (y = 0), T c = -210 for iimenite (y - i); volume of unit cell, 
V =3 100.4 + 3.8y A J , i.e., V = 100.4 A J for hematite, V = 104.2 for iimenite. The 
variation of J, with composition “y” is shown in Table 8. 
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Table 16 

MAGNETIC PROPERTIES (J n , K, Q) OF ROCKS 



•In 

k 




(in 10- 4 

(in 10"* 



Material 

Notes emu/cm 3 ) 

cmu/cm 3 ) 

Qn 

Ref. 


a 


b 


Sedimentary rocks 
Soils 

Typical 

1-100 


29 


Texas Coast 

5 



Marine sediments 

Sandy clay, Texas 10-150 
continental slope 

15-35 

-5 


Silty shale 

Ventura basin, 5-40 

Calif. 

20-Z20 

-5 


SUtstone 

Pre Cambrian, Britain; 
4 samples 


0.02-2 

6 

Clays 


20 


65 

Shale 

137 samples 

5-1478 
ave. 52 


25 

Sandstone 

230 samples 

0-1665 
ave, 20-30 


25 


Rcdbeds,Precam- 


1.6-6 

6 


brian, 9 samples 

100 


26 


Redbeds, U.S., 82 2-20 

0,4-40 

1-3 

45 


sites 

3-76 


87 


Redbeds, Wyoming, 4-29 
Trlassic 

2-13 

ave. 4.4 

21.22 


Britain, 12 sites G.2-3.3 

10-28 


75,20 


Redbeds, typical 0.5-50 



21,22 

Limestone 

66 samples 

2-280 


25 



ave. 23 
0-5 


26 

Dolomite 

66 samples 

8 


25 

Coal 


2 


65 

Typical sedimentary 1-100 

3-300 

0.02-10 


rocks* average 





Igneous rocks 





Typical igneous 
racks, average 

100-40,000 

50-5000 

1-40 


Granite 

Pluton, Yosemite 100-800 
Calif. 

1000-4000 

0.3-1 

58, 62, 24 



to 


26 


97 samples, Okla. 1000-180,000 

l 

28 

62, 24 


41 samples 

280-2000 J 




30-2700 


25 


Minnesota, 31 samples 

0-4000 


42 


<1.4% Fe 3 0 4 

ave. 470 




Without Fe 3 0 4 

1-5 


58 


Intrusives, Japan 


0.1 -0.5 

58 

Acidic intrusives 

58 samples 

3-6527 




ave. 647 
30-60 
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MAGNETIC PROPERTIES (J n , K, Q) OF ROCKS 



J n 

k 




(in 10-* 

(in 10-* 



Material 

Notes emu/cm 3 ) 

emu/cm 3 ) 

Qn 

Ref. 


a 


b 



Minnesota, 17 samples 

350 


42 

Granodiorite 

Nevada 


oa-o.2 

13,58 

Diorite 


200 


26 

Dolerite (diabase/dikes) 





sills* England, 5 
samples 


2-3.5 

6 


dikes, India, 23 

55-1100 


53 


samples 

ave. 337 



Diabase 

Typical 1900-4000 

1500-230Q 

80-1000 

2-3.5 

25 


Minnesota, 19 samples 

800-12000 


42 


<3.4% Fe 3 0 4 

ave. 2600 




dikes, Precambrian 

100-20000 

0.2-4 

62 

Gabbro 

Minnesota, 37 samples, 

80-6100 


42 


<0.9% Fe 3 0 4 

ave. 1000 




Minnesota 

2000 

1-8 

26 


Sweden 


9.5 

58 



70-2400 


25 

Intrusives 

Sudbury basin, 1000-60,000 

20-5000 

0.1-20 

62 


Ontario 

Precambrian, basic 

2000-9000 

1-2 

76 

Basic 

73 samples 

44-9711 
ave. 2596 




Precambrian, India, 
5 samples 

3675-430Q 


53 

Basait/diabase 

Minn., 64 samples 

2500 


42 

Volcanics 

Montana, Eocene, 1 1 ,000 
455 samples 

700 

~30 

62 


rapidly-cooled 


30-50 


Basalt 

Australia, Cenozoic, 2100 
127 samples 

900 

~5 

58, 62 


Minn., 37 samples 

20-8400 


42 


<2.5% Fe 3 0 4 

ave. 2950 

3000-8000 

40-9600 


58 


Iceland,Tertiary, 


6 

6 


70 samples 

NSW Australia, 2000-30,000 

Tertiary, 7 sites 



6 


India, Deccan traps, 

1000-6000 


53 


60 samples 

W. GreenIand,Tertiary 

ave. 2300 

1-39 

58 



2000 


26 


Seamounts, N.Pacific 


8-57 

71 


Seafloor, EM~7 Mo- 


15-105 

71 


hole, NE Pacific 


ave. 40 



Seafloor, mid-Atlantic 

24-2900 

1-160 

71 


ridge * 


ave. 48 
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Table 16 (continued) 



MAGNETIC PROPERTIES (J n 

, K, Q) OF ROCKS 




Jn 

k 




(in 10** 

(in lO" 4 



Material 

Notes emu/cm 3 ) 

emu/cm 3 ) 

% 

Ref. 


a 


b 



Seafloor, depth 5000-8000 

1-6 meters 

300-600 

25-45 


(and for comparison:) 


Qt c 



Titanomagnetite grains. 




x = 0.6, size 1-2 Mm 


35 

37 


10 Mm 


17 



20 Mm 


10 



Magnetite powder, 2% by wt 
dispersion in epoxy, grain 
size 1-2 Mm 


20 

37 


4—8 Mm 


5.4 



37-75 M m 


1.5 



75-150 Mm 

Titanomagnetite powder, 
x = 0.55,2% by wt disper- 
sion in epoxy, grain size 


1.3 



1—2 Mm 


39 

37 


5-15 Mm 


23 



75-150 Mm 


1 



Magnetite powder Jj d 

dispersed, grain (emu/cm 3 of 

(per cm 3 

q t 



size: magnetite) 

of magnetite) 



1.5 Mm 0.55 

0.19 

7,2 

11,51 


6 0.15 

0.19 

2.0 



19 0.041 

0.19 

0.54 



21 0.032 

0.21 

0.38 



58 0.046 

0.22 

0.52 



88 0.041 

0.21 

0.49 



120 0.044 

0.24 

0.46 


Metamorphic rocks J n 


% 


Metasediments 

Precam brian 

20-200 


76 

Granite/gneiss 

Precam brian, India, 

30-100 


53 


12 samples 

ave. 59 



Gneiss 


0-240 


25 

Slate 

Minn., 26 samples. 

0-100 


42 


<.2% Fe,0 4 

ave. 50 


76 

Greenstone 

Precambrian, 8 samples 

10-60 



Minn., 15 samples. 

40-880 


42 


<.2% Fe,G 4 

ave. 100 



Basic metaigneous 



76 


Precambrian 

200-4000 

0.5-2 

Serpentinite 

61 samples 

250-6000 
0-5824 
ave. 349 


87 

25 

25 

Peridotite 


12,500 




5000 


26 
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MAGNETIC PROPERTIES (J n , K, Q) OF ROCKS 



J n 

k 



(in 10^ 

(in 10 * 


Material 

Notes emu/cm 3 ) 

emu/cm 3 ) Q n 

Ref. 


a 

b 


Ores 




Magnetite ore 

Sweden 

1-10 

70 


Sweden, 3 1-63% Fe s 0 4 

240,000- 

490,000 

67 


Sweden,86-95% FCj0 4 

1,000,000- 

1,120,000 

67 


76% by wt Fe 3 0 4 

350,000 

47 


90% Fcj 0 4 , grain size 
0*2 x *2 mm 

400,000 

47 


Lodestone iron ore, 5,000,000- 
7 samples 35,000,000 


66 


Lodestone ore, 13,100,000 

Arkansas 

323,000 94 

10 


Lodestone ore, 11,100,000 

311,000 80 

10 


Japan 

Biwabik & Soudan, 

50,000- 

42 


U.S.; 15 “26% Fe 3 0 4 
6 samples 200,000- 

120,000 

66 


800,000 



Hematite ore 

Sweden 

330-800 

10 


Pre Cambrian 

60-750 

76 

Chromite ore 

FeCr 3 G 4 ; 27- 
58% Fe 

600-100,000 

50 

Hausmannite ore 

Mn 5 0 4 ; Sweden 

130 

67 

Pyrite ore 

Sweden 

420 

50 


Sweden 

8-400 

50 

Pyrrhotite ore 

Sweden 

60 

50 

Lunar rocks 

(inl0“*emu/gm) 

(inlO^emu/gm) 



typical 0. 1 -1000 

500-2000 0.001 -l e 


Soils 

7 samples 

1100-3500 

28 

Breccia 

15 samples 

50-3300 

28 

Anorthosite 

2 samples 

200-900 

28 


Breccia 

400 x 10" 6 emu/cm 3 

44 

Gabbro 


50 x lO^emu/gm 

28 

Basalt 

8 samples 

50-700 x 10^emu/gm 

28 



100-300 x 10^emu/cm 3 

44 

Meteorites 

(in cgs) 

(in cgs) 


Meteorites 

60 samples (48 0.05 -.3 

irons, i.e.,Ni-Fe; 

12 stony-irons) 

0*2-4 

32 


a J n is natural remanent magnetization (NRM). 

b Q n is Koenigsberger ratio for NRM, i.e + , Q n - where H £ is earth’s field (about 0*5 oe). 

c Q T is Koenigsberger ratio for TRM, i.e., Q T = ^XRM^^E w b firc JxRM* 5 TRM acquired in Earth s field* 
d Lab TRM, in H = 0.4 oe. 
e * 5 nH = 0Joe. * 
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Table 17 

SUSCEPTIBILITIES OF ROCK TYPES, CALCULATED FROM THEIR 
MAGNETITE AND ILMENITE CONTENT 25 55 


Magnetite Content and Susceptibility / cgs units 
Minimum Maximum Average Ilmemte, average 


Material 


k*10* 

V* 

kxlQ* 


kx 10* 


kx 10* 

Quartz porphyries 

0.0 

0 

1.4 

4,200 

0.82 

2,500 

0.3 

410 

Rhyolites 

0.2 

600 

1.9 

5/700 

1. 00 

3,000 

0.45 

610 

Granites 

0,2 

600 

1.9 

5/700 

0.90 

2,700 

0.7 

1000 

Trachyte-syenites 

0.0 

0 

4.6 

14,000 

2.04 

6.100 

0,7 

1000 

Eruptive nephelites 

0.0 

0 

4.9 

15,000 

1.51 

4,530 

1.24 

1700 

Abyssal nephelites 

0.0 

0 

6.6 

20,000 

2.71 

8,100 

0.85 

1100 

Pyroxcnites 

0.9 

3000 

8.4 

25,000 

3.51 

10,500 

0.40 

5400 

Gabbros 

0.9 

3000 

3.9 

12,000 

2.40 

7,200 

1.76 

2400 

Monzonitc-latites 

1.4 

4200 

5.6 

17,000 

3.58 

10,700 

1.60 

2200 

Leuctte rocks 

0.0 

0 

7.4 

22,000 

3.27 

9,800 

1.94 

2600 

Dacitc-quartz-dioritc 

1.6 

4600 

8.0 

24,000 

3.48 

10,400 

1.94 

2600 

Andesites 

2.6 

7800 

5.8 

17,000 

4.50 

13,500 

1.16 

1600 

Diorites 

1.2 

3600 

7.4 

22,000 

3.45 

10,400 

2.44 

4200 

Peridotites 

1.6 

4800 

7.2 

22,000 

4.60 

13,800 

1.31 

1800 

Basalts 

2.3 

6900 

8.6 

26,000 

4.76 

14,300 

1.91 

2600 

Diabases 

2.3 

6900 

6.3 

19,000 

4,35 

13,100 

2.70 

3600 


“ Using - 0.30 emu/cm , ; = 0.137 cmu/cm J . 




FIGURE 23. Relation of susceptibility to magnetite conum for 
some igneous and ore rocks (Preeambrian, of Minnesota)/ 1 
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Magnetite V 


FIGURE 24. Relation of susceptibility to magnetite content for 
some mctamorphic rocks (Adirondacks.) (Reference 29, from data of 
Reference 15). 


Magnetite V 


FIGURE 25. Susceptibility and magnetite content or rocks and ores. 
(Reference 7; data from Reference IS for tines, Reference 67 for stip- 
pled area.) 





Table 18 

DISTRIBUTION OF MEASURED SUSCEPTIBILITIES FOR 
MAJOR ROCK TYPES 7 ” 



Number 

Percent of samples with k (in 10"* emu/cm J ) 

Rock type 

of 

samples 

<100 

100—1000 

1000—4000 

>4000 

Basic extrusive 

97 

5 

29 

47 

19 

Basic intrusive 

53 

24 

27 

28 

21 

Granite etc. (he*. acidic 

74 

60 

23 

16 

1 

igneous) 

Metamorphic (gneiss. 

45 

71 

22 

7 

0 

schist, slate) 

Sedimentary 

48 

73 

19 

4 

4 


Note: Where the terms indicate: basic = mafic; high in iron/magnesium silicates, i.e., the 
denser and darker-colored ferromagnesian minerals; acidic = siliceous; high in 
quartz; extrusive = formed by cooling after extruding onto the land surface, or 
seafloor; typically finer-grained; intrusive * piutonie; formed by cooling at depth; 
typically coarser-grained* 
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° (I945J 

K m v ] WERNER 
lf|r formula 


r+K M fr-v> 


I + 3X M [1-VJ 


k m v 


0 


(I9AU 

PU2ICHA 

formula 


l+f irK M tl-V| 


FIGURE 27. Uariation of susceptibility with magneiite content for 
magnetite-rich ores (Sweden) with various formulas. Using K„, “ 1.5 
emu/cm J .“ * T a— with N = 4n/3. 


Tabic 19 

RELATIVE PERMEABILITY OF MINERALS AND 
ROCKS" 


Materia! 


p, (relative permeability) 

Minerals 

Quartz (diamagnetic) 


0.999985 

Calcite (diamagnetic) 


0.999987 

Rutile (paramagnetic) 


1.0000035 

Pyrite 


1.0015 

Hematite 


1.053 

Ifmssltc 


1.55 

Pyrrhotite 


2.55 

Magnetite (ferrimagnetic) 

5.0 

Rocks 

V* magnetite 



0 

^1.0 

i 

r 0.2 

1.006 

Granites < 

0.5 

1.017 

1 

J 

i 1.0 

1.04 

Basalts 1 

2.0 

1.08 

« 

3.0 

1- 12 

1 

' 5.0 

1.18 

Iron ore j 

10.0 

13 4 

\ 

20.0 

1.56 
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Depth to Curie- temperature isotherm 


Depths to the Curie temperature (T c ) isotherm, below which 
remanent and induced magnetization is lost, were given here (pages 78- 
79) for typical continental (ancient shields) and oceanic crustal ter— 
ranes. These calculations, consistent with interpreted magnetic bot- 
toms from modelling, assumed the dominant magnetic mineral was titano- 
magnetite and typical geothermal gradients are as given in Figure 48 . 

The gradients (A, B, and C) are given for titanomagnetite of compo- 
sition x = 0 (i.e. pure magnetite), 0.1, and 0.3, respectively. The 
variation of Curie temperature with composition in the titanomagnetite 
series is shown in "Figure 18" of the preceding excerpt pages. The 
T c inc reases slightly with hydrostatic pressure and therefore with depth, 
at the rate of 1. 9-2.0 °C/kbar for magnetite (Carmichael, 1982). The 
geothermal gradients in Figure 48 are based on data from Blackwell, 
Wyllie, and Roy (Carmichael, 1977) . 

Magnetite (x = 0) has a T^ of 575°C. For titanomagnetite, values 
are 5I0°C. for x = 0.1, 450°C. for x = 0.2, and 410°C. for x = 0.3. 
Typical geothermal gradients are 25-30°C/km for a global average near 
the Earth's surface; and averaged over the upper 40 km of the Earth, 
representative ones are about 17°C/km for oceanic crust and about 10°C/km 
for continental shields. 


If the Moho physicochemical boundary between crust and mantle is also 

a lower bound on magnetization (because of the minerals in the upper 

mantle being non-magnetic) , then the magnetic bottom of sources could 

be shallower than the T -isotherm where there is an appropriate condi- 

c 

tion of thinner crust or lower geothermal gradient. Interpretation of 
an essentially non-magnetic upper mantle comes from petrologic and geo- 
chemical study (Wasilewski et al, 1979), but there is work supporting 
a contrary view (Haggerty, 1979) . 


: ' * £ 
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SUMMARY 


Our objective has been to use the satellite magnetic anomaly 
field observed by Magsat, with its characteristically long-wavelength 
information, to help interpret gross crustal structure, composition, 
properties of the upper and lower crust, major geologic terranes (e.g. 
accreted continental fragments) , tectonic development , and any rela- 
tion to economic resource provinces such as mineral deposits and oil 
and gas basins. This was to assess the utility of such large-scale 
geophysical surveying for geological studies in "stable interior"-type 
geological areas. Analysis and interpretation was done for the study 
area of the central midcontinent of the U.S. The satellite data set 
was examined in conjunction with other geophysical and geological data. 

The anomaly data set has been derived by processing the ori- 
ginal NASA satellite data, by cleaning, correcting, and correlation, 
and then reducing to a l°xl° latitude/ longitude grid at a common al- 
titude of 400 km, applying wavelength filtering, and then reducing 
the data "to the magnetic pole" to aid in interpretation. 

Satellite anomalies are expressions, primarily, of broad and 
deep-seated magnetic properties and gross geologic crustal structure. 

The magnitude (amplitude) and pattern of anomalies depend on such 
considerations as: 

i) the thickness of the "magnetic lithosphere" that is magneti- 
zable. That is, the crustal zone that has as its lower bound (and 
magnetic "bottom") the Curie- temperature of the magnetic minerals, be- 
low which the rock is non-magnetizable. An additional lower bound 
may be the Moho physicochemical boundary, below which the rock-forming 
minerals may be nonmagnetic because of compositional factors. 

ii) the relative partitioning of the earth's crust into a sialic 
upper-crustal "layer", and a more mafic (more ferromagnetian/ iron-rich 
minerals) lower-crustal "layer". The latter layer is believed to be 
much more magnetic (due to higher susceptibility, and thermally-activated 
viscous remanent magnetization) , and may play a dominant role in pro- 
ducing the long-wavelength crustal anomalies. Of consequences are 
both differences in lower-crustal layer thickness, and lateral varia- 
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tious in the rock magnetization intensity. 

iii) magnetic mineral content of rocks. It is the amount of mag- 
netite (in the titanomagnetite mineral series) with properties deter- 
mined by its chemical origin and thermal history, which is most im- 
portant in producing magnetization. 

iv) relative metamorphic grade of large terranes of crustal rock, 
particularly for deeper basement rock. Granulite-grade metamorphism 
seems to provide the most favorable condition for larger magnetization. 

v) boundaries between, or major intrusions in, crustal "blocks" 
of different thickness, rock type (composition, age), physical prop- 
erties, and structural style. These borders can yield distinguishable 
anomaly patterns, such as magnetic lineaments. 

vi) the orientation of the local Earth's magnetic field. This 
determines the direction of the induced (and deep-seated viscous) mag- 
netization, and the resulting anomaly pattern will vary in appearance 
with latitude unless the date are "reduced-to-the-pole". 

vii) satellite magnetic mapping presents anomalies which have char- 
acterise:’ c wavelengths greater than 250 km. 

tie major anomalous features in the central Midcontinent are: 

i) a large magnetic high over east-central Kentucky /Tennessee, 
at the intersection of major magnetic/geologic lineaments. The ano- 
maly is attributed to a mafic basement rock complex in the mid- and 
lower-crust, which is anomalously ultrabasic and magnetic, and possibly 
also to some localized crustal thickening. 

ii) a magnetic anomaly high over NW Texas /Oklahoma, and then ex- 
tending in an arcuate moderate high up to the northwest through Mis- 
souri to Lake Michigan. The former may be due to crustal thickening, 
with the latter trend to the northeast being due also to a different 
geologic terrane (the Precambrian Mazatzal /Granite-rhyolite belts) 
and which has an anomalously more-magnetic lower crust (especially 
under Missouri) . 

iii) a magnetic low over the Mississippi Embayment/aulacogen, which 
extends from the Gulf Coast up to the New Madrid (southeast Missouri) 
seismic zone. This is believed due to crustal thinning, associated 
with the rift attempt in late Precambrian/early Paleozoic time, com- 


bined with anomalously less-magnetic lower crust. 

iv) a magnetic anomaly low over South Dakota, for which there 
is some evidence for contributions from localized crustal thinning, 
an elevated Curie- temperature isotherm, and a magnetization deficiency 
in the lower crust. 

v) a magnetic low over west-central Georgia, probably due also to 
an anomalously less-magnetic lower crust. 

vi) a tendency for the satellite magnetic anomaly pattern to fol- 
low and reflect the southern "rifted continental margin" formed in late 
Precambrian/early Paleozoic time with the opening of the proto-Atlantic. 

vii) a lack of distinguishable magnetic expression for the Central 
North American Rift system/Midcontinent Geophysical Anomaly, a long, 
deep, and major crustal structure and intrusion developed in late-Pre- 
cambrian (Keweenawan) time. This is apparently due to the opposing 
contributions of reversed and normal polarity remanence in the upper- 
crustal mafic rock seam, resulting in a reduced net magnetization, and 
to the north-south alignment of this relatively narrow structural fea- 
ture. This seam, a basalt-dominated sequence of extrusives and in- 
trusives, is interpreted as having its lower portion with reversed 
remanence, and its (younger) upper portion with normal remanence. In 
support of this, there is known to be a reverse-to-normal transition 
of the geomagnetic field in Keweenawan time. 

We have also compiled data on magnetic properties of minerals 
responsible for magnetization in the Earth's crust. Such information 
will help in subsequent modelling and interpreting crustal sources 
for anomalies mapped by satellite or other long-wavelength surveys. 

It will provide guidance for selecting properties to deduce crustal 
composition, structure, and thermal regime (as depth to the Curie- 
temperature isotherm) . 

Magsat anomaly data is a promising new aid for reconnaissance 
exploration of major geologic and tectonic features of the Earth's 
crust. It may be helpful in identifying and interpreting plate tec- 
tonic boundaries, now and in the past as at ancient collisional zones 
where continental blocks have been sutured together. If there is a 
causal relationship between such tectonism and the processes of con- 
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centration and emplacement of ore deposits, then there is the prospect 
for assistance with exploration for economic mineral resources. It 
might be noted that the (interpreted) most highly-magnetically-anomalous 
region of the mafic lower-crustal layer in the central U.S., an elon- 
gate zone underlying Missouri, underlies the region of the central 
Midcontinent which has the greatest number, variety, and richness of 
metallic mineral deposits (lead, copper, zinc, iron) . 
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ANALYSIS AND USE OF "MAGSAT" SATELLITE MAGNETIC DATA TO 

HELP INTERPRET CRUSTAL CHARACTER OF U.S. CENTRAL MIDCONTINENT 

R. A. BLACK (Iowa Geological Survey, Iowa City, Iowa USA) 

R. S. CARMICHAEL (Department of Geology, University of Iowa, 

Iowa City, Iowa USA 52242) 

NASA's MAGSAT satellite measured magnetic fields from October 1979 until June 
1980. The processed magnetic data yield long-wavelength anomalies that arise 
from crustal and upper-mantle sources, interpretation of the anomalies, done 
in conjunction with correlative data such as gravity and aeromagnetics, heat 
flow and geothermal gradients, and known nearsurface geology, can lead to bet- 
ter understanding of major deep-seated geologic structures and crustal compo- 
sition. Such work has application to resource exploration, as well as to geo- 
tectonics. As part of the NASA project to investigate and use the MAGSA T 
data, we are developing analysis techniques to help interpret the structure 
and character of the lithosphere in central North America. The region in- 
cludes the "Midcontinent Gravity Anomaly" paleorift zone (1200 km long, 80 
km wide, intruding 1.1 billion years ago) and the New Madrid rift/seismic 
zone, both of which are of paleotectonic and neotectonic interest. Our pre- 
liminary analysis of the initial MAGSAT data, combined with correlative geol- 
ogical and geophysical data, shows the utility of the satellite data for re- 
gional crustal and basement study. Work in progress supported in large part 
by NASA contract NAS5-26425. 
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Use of MAGSAT satellite data to interpret crustal geology, structure, 
and geophysical properties of the U.S. central Midcontinent 

R. S. Carmichael, R. A. Black, and R. A. Hoppin 

Department of Geology, University of Iowa 

Iowa City, Iowa U.S.A. 

NASA’s MAGSAT satellite measured magnetic fields from Oct. 1979 until June 
1980. The processed magnetic data yield long-wavelength anomalies that arise 
from crustal and upper mantle sources. Interpretation of the anomalies, done 
in conjunction with correlative data such as gravity and aeromagnetics, heat 
flow and geothermal gradients, and known nearsurface geology, can lead to 
better understanding of major deep-seated geologic structures and crustal 
composition. Such work has application to resource exploration, as well as 
geotectonics. As part of the NASA project to investigate and use the MAGSAT 
data, we are developing analysis techniques to help interpret the structure 
and character of the lithosphere in central North America. The region in- 
cludes the "Midcontinent Gravity Anomaly” paleorift zone (1200 km long, 80 km 
wide, intruding 1.1 billion years ago) and the New Madrid rift/seismic zone, 
both of which are of paleotectonic and neo tectonic interest. Our preliminary 
analysis of the initial MAGSAT data, combined with correlative geological and 
geophysical data, shows the utility of the satellite data for regional crustal 
and basement study. 
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Abstract for Ann. Midwest Meeting of American Geophysical Union 
Minneapolis Minn., Sept, 1981 (from EOS Trans. A.G.U., p. 615, 
August 1982) 

GEOPHYSICAL INTERPRETATION OF THE GEOLOGY OF THE 
CENTRAL SEGMENT OF THE MID CONTINENT GEOPHYSICAL 
ANOMALY 

R, R. Anderson 

R. A. Black (both at: Iowa Geological Survey, 

Iowa City Iowa 52242) 

The midcontinent Geophysical Anomaly (MGA) is a 
pronounced linear gravity high with flanking gra- 
vity. lows stretching from northern Lake Superior 
to central Kansas. The anomaly is created by a 
central basalt-dominated horst, bounded on the 
east and west by deep, elastic-filled basins. It 
is divided into three distinct sections, separated 
by apparent transform faults into northern, central 
and southern segments. The central segment, the ma- 
jority of which lies in Iowa, displays measured 
Bouguer gravity anomaly values ranging from greater 
than +67 mgals to less than -113 mgals. In this 
study the Bouguer gravity anomaly and aeromagnetic 
field of 19 profiles across the central segment of 
the MGA have been modeled. These models provide es- 
timates of the thickness of the basalt-dominated vol- 
canics of the central horst (locally greater than 
12,000 m), location of the high-angle faults which 
bound the horst, and determination of the geometry 
of the elastic-filled marginal basins which reach 
maximum depths in excess of 10,000 ra. This is the 
most detailed examination of this portion of the 
MGA yet undertaken. 
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Magsat Session 


Magsat Mission: Crustal Studies 


Results of the Magsat Mission MAGI 

M, Settle, NASA, and R. Langel, 

Goddard Space Flight Center 

The Magsat satellite mission obtained worldwide measurements 
of magnetic field strength with an accuracy and spatial resolution 
that surpassed all previous global surveys* In addition, Magsat was 
the first satellite to measure directional characteristics of the 
earth's magnetic ‘field on a global basis* Magsat carried a scalar 
magnetometer capable of measuring total field strength to an ac- 
curacy of +3 gammas* and a vector magnetometer capable of 
measuring component field strength in thn*a orthogonal directions 
to an accuracy of ±6 gammas. TI 

altitudes of 350-560 km over a ^se Magsat Satelli 

products are improved spherical * Data Interpret Cru 

main (core) field* and crustal an Resource Potential c 

wavelength variations in the crust * Robert S, Carmichael at 
are available through the Natio The a^essmem of futi 

Goddard Space Flight Center, G KKuthcmial eneruv . is b< 


l\ Studies Large-Scale Crustal Magnetization Models MAG3 

Derived from Satellite Magnetic Anomaly Data 

M, A, May hew, STS , Inc. 

An equivalent layer magnetization model for the U.S. and en- 
MAG1 virons derived from Magsat total field anomaly data is presented. 

Information contained in such a model is essentially lateral varia- 
tion. to a certain limit of resolution, in the vertical integral of 
magnetization from the earth's surface to the Curie isotherm. 

. . . . , Resolution for the desenbed technique and Massat data is about 

y and spatial resolution ... . ^ e , . , 

: . „ f w 200 km dipole spacing. Apparent maenetizatton vanations are due 

In addition. Maesat was „ - . . - . . . . * . t . . . 

. , . . * * . to Cune isotherm depth variations* regional vanations m crustal 

I characteristics of the . _ *1 . . . 

. - . , t magnetic properties, or to discrete mafic plutontc complexes. In 

Magsat earned a scalar „ r . r _ 

- * , the western U.S*. the first effect appears to dominate. Excellent 

field strength to an ac- . . , 4 . . 

, . _ correlations between apparent magnetization and laree-scaie tec- 

jgnetometer capable of ^ inces> incIudi the CoIorado lateau> Grande 

orthogonal d.mctmns ... £ c l_ Midcontinent rift , lhe 

Use of Magsat Satellite Magnetic Anomaly MAGS achita belt, and the 

Data to Interpret Crustal Character and 
Resource Potential of the U.S* Midcontinent 
Robert S, Carmichael and Ro\s Black, Urtiv. aflat ui 

The a>sessniem of future resources of petn»le mi. minerals, and MAG4 

geuthcmml energy, is best served by an expl*-rati»*n strategy that 


achita belt, and the 


MAG4 


Data Selection Techniques i 
of Magsat Data Over Austra 

B. David Johnson and C. N. G. 
Macquarie Uniw, Australia 
The Magsat data require critic 
selF-consistent data set suitable 
quent interpretation. The interact 
scribed in this paper, involve the 
oriented data base and a color gra 
employed to select profiles that Ii 


minimum altitude, and satisfy va* 

initial scan selects profiles which have their perigee (lowest alti- 
tude of orbit) within the survey region. Subsequent scans through 
the data are made to select profiles that contain data below a given 
altitude. Each profile selected by these automatic search criteria is 
displayed for visual validation. Various interactive procedures are 
available to remove data spikes, trim the ends of the profile, and 
detrend the data values. This corrected profile data may then be 
compared with other profiles already in the data base and then 
finally stored in the data base if selected. The large degree of 
redundancy in the Magsat data enables the rejection of noi^v or 
bad profiles and the detection of time-varying effects in the data. 
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combines geophysical exploration with geologic interpretation ot 

the structure, composition, and thermal and tectonic development * en ) s as ° n S wavt - 
of die earth's upper crust. A promising tool for reconnaissance tn • 

exploration of the crust and its energy and resource potential is the ; orfaita drift while retain- 
use of satellite sensing and measurements. NASA's Magsat sate!- oved a ^st-fit first-order 
lite measured in 1980 the earth's long-wavelength magnetic field tiie e “ ect a tltu e 
anomalies arising from major crustal sources. Acquired at altitudes road sca * anoma bes. we 
of 350-560 km, the data provide global coverage with resolution 11111011 altitude by using a 
notpreviously achievedby satellite sensing. Magsat anumalies have nlation ‘ o interpret these 
a magnitude ramie up to 26 gammas, and a magnetic-source reso- te rm s - W<J ev _ we a mo ™ 
lution of about 150-200 km. of magnet. anon ,n add.- 

*nal matnx ot reciangular 
4 ». — vu w i^tcr.m mu TuiiunMu of the magnetic suscepti- 

bility, which eventually vanishes at the Curie depth. Both of these 
new methods are verified to be effective for the simulated anom- 
alies, The results from the actual data must be compatible with 
such surface data as Magnet U.S. aeromagnetic data or a crustal 
model based on the seismic profiles* 
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be effective up to the number of block* 35 1 using a simulated data 

set for a simple theoretical anomaly. 
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Use of Magsat Satellite Magnetic Anomaly MAGS 
Data to Interpret Crustal Character and 
Resource Potential of the U.S. Midcontinent 

Robert 5 . Carmichael and Ross Black , Uniw of Iowa Summary 

The assessment of future resources of petroleum, minerals, and 
geothermal energy, is best served by an exploration strategy that 
combines geophysical exploration with geologic interpretation of 
the structure, composition, and thermal and tectonic development 
of the earth's upper crust. A promising tool for reconnaissance 
exploration of the crust and its energy and resource potential is the 
use of satellite sensing and measurements, NASA's Magsat satel- 
lite measured in 1980 the earth's long-wavelength magnetic field 
anomalies arising from major crustal sources. Acquired at altitudes 
of 350-560 km, the data provide global coverage with resolution 
not previously achieved by satellite sensing, Magsat anomalies have 
a magnitude range up to 26 gammas, and a magnetic-source reso- 
lution oF about 150-200 km. 

The satellite data, used in conjunction with correlative geo- 
physical and geological data such as aeromagnetics, gravity, 
Landsat imagery, and near-surface geology can lead to better 
understanding of major deep-seated geologic structure and condi- 
tions and the associated potential for the localization of energy and 
mineral resources. With the central U.S. as a study area, we have 
been investigating the use and applicability of the satellite mag- 
netic data in interpreting the midcontinent's crustal structure and 
composition. The anomaly data sec has been derived by processing 
the original satellite data, reducing to a t* x 1° latitude/longitude 
grid at a common altitude datum of 400 km, applying wavelength 
filtering, and then reducing the data ”to the pole” to aid in inter- 
pretation. The major anomalous features are a magnetic high over 
Kentucky, a high over northern Texas and extending to the north- 
east towards Lake Michigan, and magnetic lows over the Missis- 
sippi embavment, and South Dakota. These features are inter- 
preted in terms of variation in crustal thickness, depth of Curie- 
temperature isotherm, and petrologic variation between different 
geologic provinces. Supporting information is provided by gravity 
anomalies, and crustal thickness as determined From regional 
seismic data. Abstract 

Exploration for future supplies of energy fuels (oil. gas), min- 
eral resources, and geothermal energy, will be best pursued by 
combining geophysical prospecting with geologic interpretation of 
the structure, petrologic character, leetnme history. and thermal 
and chemical condition nt the earth's crua A promising tool tor 
aiding in reconnaissance e^pW.ifum ot cnergj and po*\ 
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inces, and potential localization of economic objectives, is the use 
of satellite sensing and measurements. 

NASA's Magsat satellite measured in 1980 the earth's magnetic 
field anomalies arising from major crustal sources. This was the 
first satellite to measure component- field n X T Y f Z) data to 
yield relatively low-altitude global magnetic jiiomuly maps. These 
maps have a resolution not previously achieved b> satellite mag- 
netic-field surveying. The satellite operated at altitudes of 350-560 
km, mapped anomalies having a magnitude range up to 26 gam- 
mas {nanoteslasj, and yielded a magnetic source resolution of 
about 150-200 km. 

The satellite magnetic data, and resulting maps, constitute a 
new aid for interpreting the crust of the earth and as a reconnais- 
sance component in a geophysical exploration strategy. They can 
be integrated with correlative geophysical and geologic data sets 
such as aeromagnetics, gravity coverage, Landsat remote sensing 
imagery, and near-surface geology. This can lead to a better 
understanding of crustal geologic structure and composition, and 
the associated energy and mineral resource potential. The long- 
wavelength magnetic anomalies measured at satellite elevations 
are expected to reflect differences in crustal thickness, the petro- 
logic character and properties of major geologic provinces, and 
variations in geothermal gradient and thus the depth of the Curie- 
temperature isotherm. 

With the central U.S, as a study area, we have been investigat- 
ing the applicability of satellite magnetic data in interpreting the 
midcontinem’s crustal geologic structure and composition. 
NAS A/ Goddard Space Flight Center has provided Magsat data as 
measured along the satellite tracks, with some preliminary correc- 
tions applied. We then did computer processing and analysis of the 
data set, with steps including: (l) removal of spurious data points, 
arising from the data acquisition process rather than from terres- 


trial sources; (2) statistical smoothing, comparison, and correla- 
tion of individual data tracks, to reduce the effect of geomagnetic 
temporal and transient disturbances; (3) reduction of data by 
weighted averaging to a grid with 1° X l a latitude/ longitude spac- 
ing, with altitudes interpolated and weighied to a common datum 
elevation of 400 km (from the orbit range of 350-560 km}; (4) 
wavelength filtering, to remove (or retain) anomaly features of 
desired wavelengths; and (5) reduction of the anomaly map "to the 
magnetic pole’* to compensate for variation of magnetic inclina- 
tion with latitude and aid in interpreting the anomalies with respect 
to causative crustal features. 

Figure 1 shows the Magsat anomaly map we extracted from the 
data for the central U.S. Here, as for subsequent figures, it is for 
magnitude (i.e., scalar} data, as calculated from the X. Y , and 
Z-component vector data as measured by the satellite. The actual 
area of processed data is about 50 percent larger than shown, and 
this belt- way around the study area was used to prevent “edge 
effects’* in the processing to follow. 

Figure 2 shows the previous map, now filtered with a low-cut 
wavelength filter which removes wavelengths less than about 400 
km. This map has similarities to maps of some previous satellite 
studies over the U.S., but is more detailed and revealing because 
of (1) the care taken in removing bad data points and tracks, and 
comparing profiles to optimize data reliability, and (2) treating the 
data on a finer 1° x 1° grid rather than the customary but coarser 2° 
x 2° grid spacing often used for global-scale maps. For compar- 
ison, Figure 3 shows the Magsat anomaly map as derived by 
NASA (March 1981) from a preliminary global data set. There is 
less detail in the latter, and it is a prime purpose of our study to 
assess, for prospective resource exploration and crustal interpreta- 
tion, how far the satellite magnetics resolution can be extended by- 
data processing and analysis. 
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s Fig. 3, NASA’s preliminary Magsat scalar anomaly map. Data 
from below 400 km altitude, on 2° x 2° blocks. Contour interval 
of 2 gammas. Zero-level is arbitrary. 


Several majtTr anomalous features are opened or the pruce^ed 
scalar magnitude map (Black. 1981 ) of Figure 4. They are: ilia 
bulls eye magnetic high over Kentucky, as noted on previous 
POGO satellite data and interpreted by others as due to a mafic 
basement rock complex and intrusion in the lower crust: (2) an 
arcuate magnetic high extending from its maximum over nonhem 
Texas up to the northeast and the Lake Michigan area, probably 
originating from a zone of greater crustal thickness developed in 
late Precambrian time; (3) a magnetic low over the Mississippi 
embayment/rift, as noted by others and having an origin asso- 
ciated with a combination of crustal thinning accompanying failed 
(paleorifring, petrologic character of the low er crust, and elevated 
Curie-temperature isotherm; and f4) a magnetic low over South 
Dakota. 

Supporting correlative information for this study of satellite 
magnetic anomaly data includes the distribution of basement rock 
provinces, gravity anomaly data, and crustal thickness as deter- 
mined from regional seismic data. For example, in the midconti- 
nent region the areas of thicker crust fi.e., thickness over 45 km) 
are those with the magnetic highs on the satellite map — north 
Texas/Oklahoma, and Kentucky /Tennessee. 

This work is supported in large part by a NASA/ Goddard Space 
Flight Center contract. 


Figure 4 shows the data of Figure 2, now reduced to the mag- 
netic pole. The alteration is not dramatic, because of the relatively 
high magnetic latitude {about 45-65° N.) of the study area with 
respect to the north magnetic pole. The major magnetic highs have 
been shifted north — the north Texas one up to the Oklahoma 
border, and the Tennessee one up into Kentucky. 
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Fra. 4. 1° x 1° scalar data at 400 km altitude, high-pass filtered and 
reduced to the magnetic pole. Contour interval of 2 gummas. 
Zero-level is arbitrary. Plot by Black (1981). 
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CRC Handbook of Materials Science is e current, readily accessible 
iwide to the physical properties of solid state and structural 
Materials. Interdisciplinary in approach and contant, it covert a 
The CRC road variety of types of materials, including materials of present 
all of the ommercial importance plus new biomedical, composite, and laser 
thesenew Materials. In view of the importance of having critically evaluated 
on Cosm roperty data available on materials in order to solve modern 
Symbols i*robiems. this handbook is published to fill the gap between 
and a Sola Merit ical data collections carrying all the information available for 
i single material class and general reference works, with only 
Updated i mited property and classification data on materials. 

Selected ' 

7 x 10. 198J tilizing a largely tabular format for easy reference and compare- 
iihty of various properties, this reference ia of uaa to those seeking 
Mformation on materials outside their specialized field of compe- 
• nee as well as for experts in the field seeking detailed property 
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GENERAL PROPERTIES 

•>0 pp.. 7 x 10. 1974, ISBN-0 87819-231-X. 

• ECTlON 1 - THE ELEMENTS. 

ECTlON 2— ELEMENTAL PROPERTIES. 

ECTlON 3- MISCELLANEOUS TABLES OF PHYSICAL PROP- 
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• Current Data 

• Tabular Format 

Reliable, comprehensive data on various properties of rocks, min- 
erals, and other related materials is compiled in an organized 
manner. Format is largely tabular for easy use in comparison and 
referencing. Chapters are written by recognized experts working 
in education, industry, and government. The information is of 
particular value to those in geology, geophysics, materials engi- 
neering, and allied fields. 

VOLUME I: Mineral Composition: Chemical composition and prop- 
erties of rocks, mineral/* and crystals, pore fluids, ores, coal, 
petroleum. Earth's crust, meteorites. Electrical Properties: Resistivity 
and dielectric constants of minerals, dry and wet rocks, sedimentary 
rock sequences. Earth's interior. Spectroscopic Properties: Absorp- 
tion/transmission, reflectance and emission, and spectral charac- 
teristics of minerals and rocks in the visible and infrared range. 
416 pp., 7 x 10, 1982. ISBN-0-8493-0226-9. 

Catalog no. 0226ZN, $69.50 
Outside U.S., $80.00 

VOLUME II: Seismic Velocities: Compressions! and shear wave 
velocities for rocks, minerals, marine sediments. Earth's crust, ice; 
variation with fluid saturation, pressure and temperature. 

Properties : Properties of magnetic minerals and rocks, and their 
variation with different parameters. Engineering Properties: Factors, 
tests, and properties relating to rock appraisal, characterization, 
and assessment of hardness, strength, and deformation. 360 pp., 
7x 10. 1982 ISBN-0-8493-0227-7. 


Catalog no. 0227ZN, $64.00 
Outside U.S., $74.00 
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Abstract for Ann. Meeting of Geological Society of 
America, Indianapolis Indiana, Oct. 1983 


EARLY PROTEROZOIC DEVELOPMENT OF THE SOUTHERN 
ARCHEAN BOUNDARY OF THE SUPERIOR PROVINCE IN THE 
LAKE SUPERIOR REGION. 

ANDERSON* Raymond R. , and BLACK, Ross A., Iowa Geological Survey, 

123 N. Capitol St., Iowa City, IA 52242. 

The earliest Proterozoic rocks in the western Great Lakes area repre- 
sent two major regimes. The oldest regime (2100-1900 Ma) is character- 
ized by sequences of supracrustal rocks, dominantly conglomerates, 
quartzites, arkoses, iron formations, and carbonates, and coeval mafic 
, and intermediate volcanics. The second regime (1840-1800 Ma) includes 

i deformation and metamorphism of these earlier rocks and emplacement of 

a suite of dominantly calc-alkaline volcanic rocks and plutons and 
deposition and metamorphism of a variety of sediments including con- 
glomerates, pelites, graywackes, tuffs, and carbonates. 

These early rock suites are described by most authors as the product 
of a single orogenic event, the Penokean Orogeny. In this study, it is 
proposed that these events represent two unrelated tectonic events. 

The first, here called the Lake Superior Taphrogen, is characterized by 
extensional tectonics, graben development, platformal and basinal sedi- 
mentation, and mafic-dominated vulcanism culminating in rifting and 
ocean development on the southern margin of the Superior Province. The 
second, the Penokean Orogeny, included the closing of the ocean, sub- 
sequent collisional tectonics, and the development of a calc-alkaline 
volcanic suite. 

The rifted continental margin and subsequent suture zone is repre- 
sented by the Niagara Fault Zone in Michigan and northern Wisconsin and 
its southwest-trending geophysical extension, the Storm Lake Trend, 
which crosses Minnesota, Iowa, and extends into Nebraska. This paper 
will address the evidence for the location of this zone and propose a 
scenario for the chronology, stratigraphy, and structural development 
of the area in Early Proterozoic time. 



